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Sur les Conditions d’Auto-Oscillation des Systemes Vibrants 


Par YVES ROCARD 


Ecole Normale Supérieure, Paris 


A Discourse delivered 5 May 1948, on the occasion of the presentation of the Third Holweck 
Prize and Medal to Professor Y. Rocard 


MS. received 5 Fuly 1948 


RESUME. Le professeur Rocard a rassemblé en une doctrine plus compréhensible les 
conditions dans lesquelles les systémes vibrants peuvent devenir auto-oscillants. On a 
depuis longtemps des moyens mathématiques trés bien adaptés pour déterminer ces con- 
ditions. [I restait a les analyser et A en comprendre la nature. On peut sur quelques 
exemples voir la complexité du sujet ; l’auteur cite le roulement des bogies, des loco- 
motives (la conicité des roues engendre |’auto-oscillation sous forme de “ lacet’’’), des 
auto-oscillations du genre des coups de bélier entretenus dans les conduites, enfin les 
vibrations d’ailes d’avion, ou “ wing flutter ’’. 

Ij faut distinguer dans ces systémes ceux qui résultent du couplage d’un certain nombre 
de circuits oscillants ou pendulaires ayant chacun leur période propre (le “‘ wing flutter ”’ 
est de cette sorte) et ceux qui résultent du couplage d’éléments n’ayant pas de période 
propre, mais seulement des constantes de temps: les télécommandes sont de cette sorte, 
et aussi les systemes d’équations un peu arbitraires que |’on peut construire dans les théories 
économétriques. On peut établir alors les résultats suivants : 

1. Si un systeme n’a qu’une fréquence propre, les forces qui peuvent le rendre auto- 
oscillant ont forcément le caractére d’une résistance négative. 

2. Si un systeme résulte du couplage de plusieurs “‘ degrés de liberté ” ayant chacun 
leur fréquence propre, tout couplage passif (par self, capacité, inertie, etc.) écarte les 
‘fréquences propres composantes les unes des autres pour donner les fréquences prcopres. 
du systeme, mais il existe des couplages réactifs qui peuvent rapprocher les dites fréquences 
propres, sans l’introduction d’aucune résistance positive ou négative. On démontre alors 
qu’un tel systeme devient auto-oscillant au moment oti deux fréquences propres se 
confondent. 

On peut imaginer dans les problemes de “ wing flutter”? des cas trés simples qui 
illustrent ce théoreme, et on peut en donner des modeéles électriques également trés simples. 

3. Dans un systeme tel que ci-dessus (§ 2) mettons maintenant des forces de résistence 
(proportionnelles aux vitesses) dans l’espoir d’amortir les auto-oscillations. On démontre 
ce résultat surprenant que, tant que les fréquences propres ne sont pas profondément 
modifiées, ces résistances ne font que faciliter l’auto-oscillation, qui prend naissance alors 
que deux fréquences propres se rapprochent mais ne sont pas encore confcndues. 

Ces résultats font d’ailleurs perdre a la théorie des oscillations de relaxation beaucoup 
de son importance philosophique, le mécanisme de la résistance négative n’apparaissant 
plus comme obligatoire dans l’auto-oscillation. nek 

4. Dans le cas de systemes tout a fait généraux, ot on ne peut distinguer des fréquences 
propres composantes, il n’est naturellement pas possible de se livrer A une analyse aussi 
pénétrante, mais on démontre encore que si l’on attaque le systeme par une petite perturbz- 
tion sinusoidale dont on fait croitre la fréquence, toutes les conditions de stabilité connues 
(Routh, Leonhard, Nyquist) expriment une vérité tres simple: quand, 4 impédance 
croissante avec la fréquence, on passe par une résonance (la réactance s’annulant) le terme 
de résistance trouvé doit étre positif. 
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f dois tout d’abord exprimer mes remerciements a la Physical Society de 
Londres pour la haute distinction dont elle m’honore aujour@’hui. J'ai 
eu la rare chance, pendant la derniére guerre de pouvoir utiliser pour la 
défense de la Liberté mes propres possibilités scientifiques. Cette action 

comportait en elle-méme sa récompense, et pour celle qui vient s’y ajouter aujourd’ 
hui, je ne puis qu’en reporter le mérite sur ceux de nos collégues qui ont perdu la 
vie du fait de ’ennemi. Outre notre grand Holweck, je citerai spécialement 
MM. H. Abraham, E. Bloch, G. Bruhat. Tous trois ont dirigé le laboratoire de 
Physique de l’Ecole Normale Supérieure, dont je suis maintenant chargé, et 
cette Ecole s’honore aussi d’avoir vu beaucoup de ses éleves combattre pour la 
Liberté, aux cdtés de ses meilleurs défenseurs, les Anglais. Elle sait de quelles 
traditions elle est la gardienne, et elle ne faillirait pas devant des taches nouvelles 
qui se présenteraient. 

Pendant la guerre, l’activité scientifique dans les pays occupés était pro- 
fondément affectée par des circonstances en général facheuses, notamment 
une certaine asphyxie intellectuelle qui venait du manque de bibliographie en 
langue anglaise. Cependant certains savants, dont l’existence antérieure n’était 
peut-étre pas assez repli¢e dans la méditation, ont trouvé dans les loisirs forcés de 
Europe Nouvelle l’occasion de mener a bien certaines taches de longue haleine. 
Ainsi, plusieurs de mes illustres maitres ont séjourné quelque trois semaines dans 
la prison de Fresnes. Ils ont assuré en étre sorti les uns et les autres avec la 
solution enfin trouvée de problemes mathématiques qui leur avaient paru jusque- 
la insolubles .. . 

Bien que n’ayant pas été soumis a un traitement aussi efficace, c’est un peu de 
cette maniere que mes propres travaux ont évolué, ce qui me fournit tout naturelle- 
ment le sujet de cette ‘‘ Address”’. 

Dés 1934 ,j’ai eu a m’occuper‘de questions concernant la stabilité des systemes 
mécaniques, et des conditions dans lesquelles ils deviennent auto-oscillants. 
Nous conviendrons donc aujourd’hui de rassembler sous le méme vocable 
d@ “oscillateur”” des objets aussi différents que 

— une locomotive qui a du lacet 
— une lampe triode qui oscille 
— une flute ou un tuyau d’orgue mis en mouvement par un jet d’air 
— le wing-flutter des avions 

etc... et nous chercherons ensemble quelles sont les causes qui provoquent 
l auto-oscillation. 

La locomotive est l’engin le plus difficile: on ne comprend méme pas comment 
elle roule. Le Professeur Andrade est, je crois, l’auteur d’un petit livre sur les 
locomotives, il saura me redresser si mon exposé contient des erreurs :—Les roues 
des locomotives sont coniques, le rayon de roulement varie donc avec l’écart 
latéral de lessieu sur la voie. Posons donc sur une voie un essieu avec ses deux 
roues, un peu de travers (figure 1), et poussons en avant, on voit sans trop de peine 
que l’essieu roule en avancant d’un mouvement sinusoidal entretenu pur: 
loscillation est pendulaire, mais sa cause est purement cinématique. Sicependant 
l'on pose sur une voie un bogie avec ses deux essieux et ses quatre roues, les deux 
sinusoides qu’ont envie de décrire en roulant les deux essieux sont incompatibles 
avec les liaisons (les deux essieux doivent rester paralléles, etc...). 
te Un simple bogie ne peut donc méme pas “‘rouler”. Et cependant, il a bien 
Vair de rouler. Quel est donc son mouvement ?— On sait depuis Reynolds que 
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des déformations des rails et des roues autour de aire de contact finissent par 
rendre le mouvement possible moyennant l’intervention de certaines forces 
proportionnelles aux dites déformations. Ces forces ne dérivent pas d’un 
potentiel, elles peuvent soit amortir les oscillations latérales propres, soit les 
auto-entretenir: c’est généralement ce qui se passe au-dessus d’une certaine 
vitesse critique, et j’ai autrefois rempli deux petits volumes de calculs de telles 
vitesses critiques de locomotives. Je calculais trés bien mais je comprenais trés 
mal, lorsque les malheurs de la guerre obligérent la généralité des Francais a 
s'intéresser aux bicyclettes. 

Or, on y trouve l’exemple le plus simple possible d’un oscillateur ot la cause de 
Tauto-oscillation réside dans les liaisons non holondmes du roulement. Soit en 
effet le petit chariot représenté sur la figure 2, et largement simplifié: il est articulé 
en O, et toute la masse m est rassemblée en O. Un systeme de ressorts tend a 
mettre le bras J, dans la direction du bras l,. La masse m tend alors 4 osciller, 


Figure 1. Figure 2. 


en déplacement latéral, comme un pendule. Cette oscillation fait changer 
Porientation des roues. Si alors on suppose que la machine se meut en avant 
{V>0) ou en arriére, les liaisons du roulement réagissent sur le mouvement 
latéral de m, lequel est régi finalement par |’équation: 


A ROR IO ar ISL g aN wiry eae 
wt (e-7) "5 1 nti) 7H are Cae ice (1) 


x tant l’abscisse latérale de m; K est le coefficient du couple Ky exercé par le 
ressort lorsque les bras /,, J, font entre eux un angle i. 

Cette équation est trés instructive: (i) Pour V =0 elle donne un mouvement 
pendulaire simple, sans amortissement. (ii) Sil’engin roule(V 40) on aun terme 
’amortissement, si le centre de gravité est en arriére du milieu quand on se tourne 
vers l’avant (donc, /,<J, si V>0 et inversement /,>/, si V <0). (iii) Dans le cas 
inverse, mettons /, </, et V <0 par exemple, on a un frottement négatif procurant 
un auto-entretien. (iv) Si la vitesse grandit encore on a, lorsque V?/l,/, ’emporte 
sur (K/m)(1//,+1/l,) un nouveau type dinstabilité, avec écart exponentiel sans 
auto-oscillation. Si en particulier /,=/,, l'amortissement positif ou négatif dis- 
parait, et c’est seulement ce type d’instabilité qui intervient. Nous le retrou- 


verons dans d’autres problémes. 
20=2 
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On voit donc que la locomotive était un exemple inutilement complique, et 
que des liaisons de roulement pur suffisent pour nous donner de beaux exemples 
d oscillateurs. 

Je ne m’étendrais pas davantage sur les oscillateurs a un seul degré de liberte, 
dont l’exemple le plus connu est celui de la triode oscillatrice, car en en étudiant 
d’autres, (une aile d’avion vibrant en torsion seulement, une anche dans un tuyau 
d’orgue) nous retrouverions toujours ce résultat. 

Un systeme a un degré de liberté devient auto-oscillant par une résistance négative. 

Si lon aborde maintenant les oscillateurs 4 plusieurs degrés de liberté les 
conclusions deviennent toutes différentes, et nous découvrons de nouvelles causes 
d’auto-oscillation. Quoique la théorie soit trés générale, nous parlerons seulement 
pour deux degrés de liberté. Commengons par examiner les moyens mathé- 
matiques a notre disposition. En mettant en équation ces problemes, nous 
trouverons généralement une équation différentielle portant sur une des variables, 
y, il sera suffisant de la considérer comme linéaire au moment out l’on sera a la 
limite de stabilité. Elle s’écrira donc symboliquement 


Aop" +a,p""1 abe be See Gh, == (); Se. 0 eels (2) 


p étant le symbole d/dt, et devenant aussi un paramétre algébrique lorsque l’on 
pose y=e”", et le systeme correspondant sera stable si l’équation (2) n’a que des 
racines a parties réelles négatives. 

Les conditions d’une telle stabilité s'expriment par divers moyens, trés connus, 
et en principe équivalents : le critérium de Routh (ou les déterminants de Hurwitz) ; 
le critérium de Nyquist. 

Mais on fait ainsi des calculs aveugles, sans analyser les mécanismes physiques 
a la base. 

Soit donc un systeme a deux degrés de liberté: sur la figure 3, on a représenté 
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Figure 3. Figure 4. 


deux circuits L, C et un organe de couplage A. L’ étude de ce systéme nous méne 
a une €quation dop* + a,p? + ap” + a3p +a,=0, et les conditions de stabilité sont, 
d’aprés Routh, les suivantes: 


(@,/@y) — (a@3/a,)>0 et en outre, a, az, a,>0, 


ay/ Ap 


—_——— > 0 ainsi que a,, a.>0. 
@y[Ay — A3/a; q vee 


(a3/a,) — 


Mio: (3) 


Elles sont bien classiques, mais on les écrit rarement sous cette forme. Nous 
avons voulu au contraire dans les équations (3) mettre en évidence que les termes 
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ampairs 3, @, Mintervenaient que par leur rapport. Or ces termes impairs 
viennent uniquement des résistances qui peuvent figurer dans les circuits, et si l’on 
convient de faire tendre les résistances vers zéro en leur gardant leurs valeurs 
relatives, les équations (3) montrent que les conditions de stabilité ne changent pas. 
D’ou la question suivante: wn systeme totalement denué de resistances, positives 
ou négatives, peut-il étre auto-oscillant ? 
Cela revient 4 étudier les racines de Péquation: 


yp -Osp" +a, =U) ee ee (4) 


Sur la figure 4, nous avons représen té & la fois le plan complexe de la variable p?, 
et le plan complexe de la variable p. 

x et B sont les deux racines en p? (deux nombres réels négatifs) lorsque 
4,” —4a)a,>0. Mais ces deux racines passent en «’, ’ lorsque a,” —4a a, <0. 
Elles sont alors imaginaires conjuguées. Extrayant leurs racines carrées pour 
obtenir p;, Ps, Ps, Ps, les quatre racines de l’équation (4) on en trouve deux (p,, po) 
a partie réelle positive, en divisant par deux l’angle argument de «’ par exemple. 
Dans ce cas, |’équation (4) décrit un mouvement instable. 

Or, le cas limite étant obtenu pour a,?—4aya,=0, lorsque l’équation (4) avait 
ses deux périodes propres confondues. Nous avons donc obtenu la proposition 
suivante : 

La propriété est visiblement générale, on peut la démontrer pour un nombre 
quelconque de degrés de liberté, et méme pour une infinité (oscillations avec 
propagation, lignes électriques, tuyaux sonores...) et énoncer: 

Dans un systeme a nombre quelconque de degrés de liberté, dépourvu de résistances 
positives ou negatives, Vinstabilite nait avec la confusion de deux fréquences propres. 

Si l’on cherche des exemples de tels systémes, on trouve que tous les couplages 
passifs écartent les fréquences propres, donc ne peuvent donner naissance a 
Pinstabilité. C’était bien nécessaire pour la conservation de I’énergie. 

Ainsi, sur la figure 5, les deux pulsations propres sont (dans le cas simple ot 


Figure 5. Figure 6. 


L,=L,, C; =C,) sans couplage (M=0) w, =w2=1//(LC), avec couplage (M0) 
w,=1/\/{((L4+M)C}; w.=1/./{(L—M)C}. Ces couplages passifs sont d’ail- 
leurs caractérisés par la symétrie des coefficients dans les deux équations. 

Un systeme mécanique équivalent au schéma de la figure 5 est fourni par une 
aile d’avion qui vibre en flexion-torsion (figure 6). 

Laxe elastique E est le lieu des points ot l’on peut mettre des poids sur I’aile 
sans latordre. L’angle «, décrit la flexion de l’aile. L’angle «, décrit la torsion. 

Si Je centre de gravité G n’est pas sur l’axe élastique on concoit sans peine 
qu’en empoignant l’aile en A et en la levant brusquement, la force d’inertie 
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appliquée en G va l’obliger a se tordre d’un certain angle «2: flexion et torsion: 
sont donc couplées et, avec les variables choisies, le couplage est par inertie. On 
obtiendra donc finalement deux équations du type: 


2 a ax5 7 
(1 as +K;) al ee 


1 dt? dt? 
ee (5) 


12 dt 


Par exemple, pour un avion de chasse, on peut avoir: 


Cee ee ) 
{tae castes IK =(). 
If + (22 TE + Ky) a, | 


Sans couplage (1;,=0) Aprés couplage (14.0) 
Fréquence propre de flexion 13 cycles/seconde 10 cycles/seconde 
Fréquence propre de torsion 24 cycles/seconde 30 cycles/seconde 


Faisons maintenant voler l’ avion qui posséde cette aile, a une certaine vitesse V, 
elle va devenir un oscillateur, et d’ailleurs, cesser presque immédiatement: c’est la 
le ‘wing flutter’’, et voici comment on I’explique. 

En vol, des forces aérodynamiques s’exercent, ces forces sont modifiées par la 
flexion et latorsion. Les nouvelles forces issues de la flexion sont négligeables (et 
ont d’ailleurs le caractére d’un amortissement), celles issues de la torsion sont 
proportionnelles a « (angle qui vient augmenter l’incidence) et au carré de la 
vitesse. Elles peuvent étre appliquées au voisinage de l’axe élastique, auquel 
cas elles ne réagissent pas sur la torsion elle-méme, mais, en tout cas, elles exercent 
un moment de flexion considérable (force F sur la figure 6) de sorte que la premiere 
des équations (5) se récrit ainsi (8, coefficient) 


da do2 5 
(1.5 de +Ky) ont In Ge dE =pV2 Op, rw wwe (6) 


donnant naissance au déterminant en p: 
Tp? +k, Typ? — BV? 
Tyop* yp? 0 Ky 


On trouve sans peine que si le centre de gravité est en arriere de l’axe élastique 
(1,2 <0) il existe une vitesse critique, on s’en explique l’origine ainsi: 

La seconde équation (5) montre que, aux fréquences w basses (K,>I,w?), a, 
est en opposition de phase avec %(I4,<0, d?/dt? = —w*). Dans l’équation (6), le 
terme —J,,w*x, vient donc diminuer le rappel élastique K,«, mais le terme 
—PV*ax. va ’augmenter, la fréquence propre de flexion va donc augmenter sans 
cesse avec la vitesse et se rapprocher de la fréquence de torsion: lavion meurt 
quand elles se confondent. 

On trouve facilement I’ equivalent électrique dun tel systeme: les équations (5) 
conviennent a la figure 5 si «, et « représentent 
les charges des condensateurs C,, Cy. On réalise 
le cas de léquation (6) par le schéma de la 
figure 7, ot un amplificateur A attaqué par une 
tension en «, la transforme en une force électro- 
motrice réinjectée dans le circuit de gauche (sous 
la forme d’un tres fort courant passant dans une 
trés petite résistance). Le gain de A est ici 
Panalogue du parametre V? pour I’avion. 
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La question de l’autc-oscillation par confusion des fréquences propres est 
donc bien claire, elle exige d’ailleurs des forces qui ne dépendent pas d’un potentiel, 
mais qui n’ont pas a étre en phase avec des vitesses, la résistance négative n’est pas 
le seul mécanisme qui améne & l’auto-oscillation, comme on avait tendance A le 
croire aprés les beaux travaux de M. van der Pol, continuant Sir Edward Appleton, 
sur les oscillations de relaxation. 

Mais les systémes réels sont beaucoup plus compliqués que ceux que nous 
venons de décrire, et il est nécessaire d’ex:miner l’effet des résistances sur la 
stabilité. 

Introduisant alors des résistances dans un schéma tel que figure 7 nous 
trouvons ce tres curieux théoréme: 4 fréquences propres constantes, |’ zntroduction 
de resistances positives ou négatives ne fait qu’ augmenter latendance al’ auto-oscillation: 

Sil s’agit des circuits de la figure 7, l’amplificateur accrochera pour une valeur 
plus faible du gain; s’il.s’agit de Vaile de la figure 6 elle cassera pour une vitesse 
plus faible si sa structure est amortie. 

Pour le voir, reprenons l’équation du quatriéme degré, en p, et les conditions 
(3). 

Posons @3/a,=d,/a), une premiére condition de stabilité est 0 <= ae 
quant a la seconde, elle s écrit 
a> A4/ Ap TRS | a; 


; >0 ou +> ens Meshes 7. 
Qy — (@3/@y)\ 1—p) ay il) do (7) 


Or 1/(1 — 2) est constamment supérieur 4 4 quand p varie de 0 4 1, donc la limite 
@instabilité a)? =4a9a4/u(1 —) est atteinte pour des valeurs des paramétres telles 
que a,” >4aoa,, c’est-a-dire avant que les deux fréquences propres, en train de se 
rapprocher, ne soient confondues. 

On trouve sans grande peine des démonstrations tout a fait générales de cette 
propriété, pour un nombre quelconque de degrés de liberté. Cela revient a 
énoncer le théoreme suivant: 

Sit un polynome en x*, a coefficients tous positifs, a une racine double negative 
(donc quatre racines imaginaires pures pour x, valant + jw), et si on lud ajoute un 
polynome quelconque ne contenant que des termes de degre impairs en x, le polynome 
resultant aura au moins une racine a partie reelle positive. 

On trouve aussi un théoreme semblable pour les équations fonctionnelles 
(et non plus algébriques) auxquelles donne lieu I’ étude de la stabilité des systemes 
a propagation. Il y aurait beaucoup d’exemples a développer ot lon verrait 
comment la théo1ie précedente éclaire les conditions quiregnent a l’auto-oscillation, 
mais le cadre de cette conférence générale ne nous permet guere de nous étendre. 

Consacrons plutSt les instants qui nous restent a déduire encore quelques 
conséquences nouvelles concernant le sens des conditions qui expriment la 
stabilité des systemes. 

Soit H(,) ?équation différentielle ou fonctionnelle qui exprime l’évolution de la 
variable y caractéristique du systeme H,,).y=0. 5Silon pose y=e?"”, on obtient 


R+jX=0, wees 


R et X devenant deux fonctions de w; |’équation (8) en général nest pas 
satisfaite, mais si exceptionnellement on peut trouver une valeur de w telle que 
R(w) = X(w) =0 c’est que les conditions concernant la limite de stabilité sont juste 


satisfaites. 
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Soit en particulier un pendule amorti 
m(d2y/dt®) + f(dy/dt) + Ky =0, 
R=K —mo?%, A= fo, 


la limite de stabilité exigerait f=0, et si f était fonction de la fréquence, il faudrait 
plus particuliérement f= 0 pour w =a, tel que K—mw,2=0. On remarque aussi 
que en faisant varier w, lorsqu’on paste la resonance K —mw*=0, on trouve une 
resistance fw positive. 

Il est A remarquer que l’expression Rijw, ici K/jw +jmw, est la reactance du 
systéme, et qu’elle est une fonction croissante de la fréquence. (Il y a un théoreme 
qui l’exige.) 

Or, il y a de fortes raisons de penser qu’un systtme quelconque, pour les 
fréquences voisines de sa fréquence de résonance est assimilable a un pendule avec 
une certaine masse, une certaine élasticité et un certain frottement (positif ou 
négatif): 

En effet, au lieu de laisser le systéme au repos, attaquons-le par un 
signal ae/"’, mesurons y en amplitude et en phase, de cette mesure résultent 
R(w) et X(w). Attachons-nous a une valeur w, de w telle que R(w),=0, 
qui donne X(w,)=X,=f,w, essayons les fréquences voisines w,—Aw donne 
R(w, — Aw) = — Aw(dR/dw) (mesuré). Essayons de poser R= K, —m,w* pour » 
voisin de w, on aura 

(Ds kee 
(2) —2m,,Aw = —Aw(dR/dw), (connu), 
d’ou mw, = —4(dR/dw), et K, = —4w,(dR/dw),. 


On trouve donc K, et m,. Sices deux termes ont des valeurs positives, i] est clair 
que, physiquement, si le systéme est stable, c’est que le frottement f, doit étre 
positif, donc X, doit étre trouvé positif. Or K, et m, sont positifs si (dR/dw), est 
négatif, donc la condition de stabilité que nous trouvons est la suivante: 


L’expression X(dR/dw) doit étre négative pour toute valeur w,, ws, .... de w 
qui annule R comme R ne dépend que de w?, on peut, plus simplement encore dire: 
AAI C03") 20 (9) 


En termes physiques, on dirait plutot: la condition de stabilité est qu’on 
trouve un frottement X positif a toute résonance de réactance R=0 parcouru dans 
le sens des inrpedances croissantes (sinon, on trouve une masse négative, une élasticité 
négative, et la stabilité veut alors un frottement négatif, tous les signes de l’équation 
étant alors changés). 

Les remarques qui préc?dent, extrémement simples, nous donnent alors la clé 
des conditions mathématiques déja connues, mais dont on comprend mieux 
maintenant la nature. 

Ces conditions connues, pour la stabilité, sont (i) le critérium de Nyquist, 
(ii) Pentrelacement des racines de R et de X. Il a déja été démontré que ces 
conditions sont équivalentes. Or notre condition (9) est équivalente a l’entre- 
lacement : En effet a chaque racine successive w,, w. de R, dR/dw? change de signe, 
donc, dans le cas de la stabilité, X(w,), X(w») etc. .. . sont alternativement de 
signes contraires, donc les racines de X et de R sont entrelacées. 

Comme exemple de problemes de stabilité traités suivant ces méthodes, nous 
pourrions citer beaucoup de problémes de télécommande, mais sans doute la 
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plupart des personnes qui m’écoutent prendront-elles plus d’intérét 4 la solution 
de que!ques problémes concernant la stabilité des automobiles. 
J’en citerai rapidement deux: 


1. Stabilité de route dune voiture qui roule, direction bloquee. 

La figure 8 représente l’avant de cette voiture, vu de dessus, roulant normale- 
ment. Les aires de contact des bandages avec le sol sont hachurées. La figure 9 
représente ce qui se passe si une force latérale agit sur la voiture, les roues se 
déforment, sans que l’aire de contact change de place, mais elles vont du fait de 
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Figure 8. Figure 9. 


Figure 10. Figure 11. 


< 


leur déformation rouler dans une direction oblique. L’angle « est proportionnel 
ala force F. 

La prise en considération de cet effet permet la mise en équation du “ roule- 
ment’’ compte tenu des oscillations du rectangle (1) représentant la voiture par 
rapport au rectangle (2) des aires de contact, ce rectangle (2) étant lui-méme 
susceptible de “rouler”’ dans une direction oblique suivant les forces qui tirent 
sur lui, et qui viennent du déplacement relatif du rectangle (1) (figure 10). 

La résolution de ce probleme m’a conduit aux résultats suivants . 

(a) Si les quatre pneus sont également gocnflés, la voiture ne présente pas de 
“‘vitesse critique d’instabilité”’ si, le centre de gravité tombe entre le 
milieu et un certain point un peu sur I’avant. 

(b) Si le centre de gravité est, pour simplifier, pris au milieu de la voiture, 
celle-ci ne présente pas de vitesse critique d’instabilité si les pneus avant 
sont un peu moins gonflés que les pneus arriére, dans une proportion 
d’ailleurs limitée. 
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Telles sont les meilleures conditions de stabilité d’une automobile. J’espere 
que la prospérité bientét retrouvée dans nos deux pays nous permettra a tous de les 
expérimenter.... 

2. Le Shimmvy des voitures comportant un essteu avant. 

Ce vieux phénoméne, qu’on tend a oublier, était bien génant en 1928 environ. 
En 1948, il reste intéressant, parce que c’est un cas pur d’oscillateur selon nos idées, 
avec l’auto-oscillation par confusion de deux fréquences propres. Le figure 11 
montre un essieu avant, avec deux roues, suspendu sous une voiture. Les masses. 
et les ressorts de la suspension, des pneus, donnent une période propre d’oscillation 
pour l’angle ¢. Mais les roues roulent vers |’avant, et vont créer un effet gyro- 
scopique. Si par un mouvement ¢, la roue de gauche s’éleve, le plan des deux 
roues change d’orientation, l’effet gyroscopique le fait tourner de 4, la voiture 
roule en cercle, de part les liaisons du roulement. I) nait alors une force centrifuge F 
qui pousse la masse MW vers l’extérieur et tend a donner une variation contraire a 
dé; voila le couplage dissymétrique qui rend le systtme auto-oscillant: Une 
rotation 0 produit un couple sur ¢ (en plus de l’effet gyroscopique avec son couple 
en d0/dt agissant aussi sur 4), tandis qu’une rotation ¢ ne produit sur ¢ que le 
couple gyroscopique. 

Le nouveau couple introduit est proportionnel a V?, on voit sans peine quwil 
confond les fréquences propres en 0, ¢ a une certaine vitesse critique V., dont on 
obtient l’expression. 

Si Pon réduit le Shimmy a ces remarques tres simples en le schématisant un 
peu, on y trouve une seule cause dinstabilité, et on a un schéma tout a fait calqué 
sur celui de Vaile d’avion de la figure 6. 


CONCLUSION 

Sur les exemples développés au cours de cette lecture, nous espérons avoir 
montré qu'il était possible de se livrer a une véritable analyse des causes des. 
auto-oscillations, au lieu de se borner a une discussion mathématique rigoureuse 
moins aveugle de leur stabilité. On trouvera une doctrine compléte exposée a 
ce subjet dans l’ouvrage ‘“‘ Dynamique générale des vibrations’’ paru en France 
en 1943 et dont une seconde édition plus étendue s’imprime en ce moment. 
Une traduction anglaise est en cours, sans qu’on puisse encore préciser sa date de 
parution. 

Je tiens en terminant a adresser un vif remerciement a la Physical Society de 
Londres pour m’avoir donné cette occasion tres flatteuse pour moi d’exposer ces. 
quelques travaux. 


The Growth of Crystals 
By.G. =P; THOMSON 


Imperial College, South Kensington, London 
Thirty-second Guthrie Lecture, delivered on 4 June 1948 


N choosing a subject for the Guthrie Lecture it seems important to consider 
the use to which such a lecture can be put. There is so much published now 
on physics that any addition to the literature needs justification. I have no 

great discovery to report, but there is a field in which considerable work has been 
done at Imperial College both in the Physics Department and in that of our 
President *, and which I think it may be profitable to bring to your notice because 
it seems very promising, and it is one which has not perhaps received due attention 
in this country—the growth of crystals. 

While the last thirty years have seen an enormous increase in knowledge of how 
the atoms are arranged inside crystals, there has not been a corresponding advance 
in understanding of the kinetics of the process of growth, or of the causes which 
determine outward form. 

Certain thermodynamical arguments have long—perhaps too long—dominated 
this study, and though I hope to show that they are liable to misuse, it is best to 
start from them. 

In general a crystal must start from small beginnings, and the same sort of 
argument which governs the growth of water drops in a mist or a Wilson chamber 
applies here. Free energy resides in the surface separating two phases and retards 
the formation of a second phase inside the first. On dimensional grounds a 
surface effect must be more and more important compared with a volume effect the 
smaller the volume. Hence for very small drops the effect of surface tension is 
greatly to increase the vapour pressure with which the drop is in equilibrium, 
according to the well known formula ln p,/p,, =2va/rkT =8, where o is the surface 
tension, r the radius of the drop, v the volume occupied per molecule of the liquid, 
and 6 is the supersaturation measured as the logarithm of the ratio of the vapour 
pressure in equilibrium with the drop to that in equilibrium with a flat surface. 
The drop is of course unstable. _ Increase of radius makes the vapour oversatur- 
ated with respect to it so that it grows more. 

If a crystal is started without a previous nucleus we should expect a similar 
expression toapply. Itis true that Gibbs (1906) pointed out—“..... We cannot 
imagine a body of the internal structure and external form of a crystal to grow or 
dissolve by an entirely continuous process, or by a process in the same sense 
continuous as condensation or evaporation between a fluid and a gas, or corres- 
ponding process between an amorphous solid anda fluid. ‘The process is rather to 
be regarded as periodic.”’ Kossel (1934) has raised a similar objection. Yet by 
considering the growth of a complete layer over the surface of a small crystal we 
can derive an expression similar in form to the above with the modification that 
the linear parameter corresponding to 7 can only increase by discrete steps corre- 
sponding to the addition of layers one cell thick all over the crystal, while the 

* Professor G. I. Finch. 


+ I am indebted to Mr. Humphreys-Owen for calling my attention to this and the following 


juotation. 
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quantity corresponding to o may have to be averaged for different faces if the crystal 


includes more than one “‘form”’. 
Stranski (1928) has shown that for a cube of side 27, 


1 PE GR IS 
nb = eal; + an) 
where the second term is due to edge energy. 

We may digress here to note how, from the thermodynamic point of view, the 
equilibrium shape is determined. As P. Curie long ago pointed out, the crystal 
truly in equilibrium with solution or vapour will take a shape such that the free 
energy is a minimum for a given volume. ‘Take some point P, inside the crystal. 
Let p, etc. be the perpendiculars from P on the faces of the crystal of area F’ etc. 
Then the volume V is the sum of the pyramids with P as vertex and the faces as 


n nv 

bases : V=42p,F,, but the free energy A = Xo,F, where o, etc. are the surface 
1 1 

tensions. Now let one set of faces grow at the expense of another. ‘The changein 


volume can be written either as SFdp or 1>(Fdp +pdF). If we keep the volume 
1 1 


nh 
constant both are zero, and hence XpdF=0. But the condition that A is constant 


1 
n 


gives LodF =0 (neglecting any variation of o with F), so that faces for which it is 
1 


possible to vary F independently must be so situated that poco. ‘This implies 
that the point is taken midway between each pair of faces with equal c, and when so 
taken p/o is constant for all possible forms. This result is due to Wulff (1901). 

But all this has its limitations. When we consider quantities, we find that for a 
crystal of a few cubic centimetres the difference in the equilibrium concentration 
between one form and another, e.g. (100) and (110) for rocksalt (Brandes 1927) is 
only of the order of lin 108. (Notice that, like all surface energy effects, it vanishes 
_as the crystal grows indefinitely). This is too small to make any appreciable 
difference except perhaps for the slowest possible growth such as might occur in 
geological time. In any laboratory process small fluctuations of temperature or of 
concentration due to diffusion will far outweigh this. 

What matters is the rate of growth, and this may well be governed by some 
initial step. Gibbs (1906) early recognized this : ‘Since the difficulty in the 
formation of a new layer is at or near the commencement of the formation, the 
necessary value of the Potential may be independent of the area of the side except 
when the side is very small. ‘The value which is necessary for the growth of a 
crystal will however be different for different kinds of surfaces, and will probably 
generally be greatest for the surfaces for which the Surface Tension is least.... 
Asa final result, a large crystal will generally be bounded by those surfaces on which 


the deposit of new matter takes place least readily. If one kind of surface cannot _ | 
form a closed figure, a crystal will be bounded by two or three kinds of surfaces. | 


determined by the same condition. The kinds of surfaces thus determined will 
probably be generally those for which the Surface Tension has the least values. 


But the relative development of different kinds of sides will not be such as to make_ | 


the sum of the Surface Energies a minimum.” Gibbs’s conception has been _ 
-developed by Kossel (1928) and Stranski (1943) and their pupils in the following | 
way. Consider a simple cube lattice built up of cubes of matter and suppose for | 
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simplicity that the only forces are those between the faces of the little cubes that 
come into actual contact. 

Suppose the crystal has grown till its side is n units and that it is a cube. Then 
of the 6n° faces all will be in contact except the 6v? on the surface. If ¢ is the 
binding energy for each pair of faces, the average binding energy per unit is 
(3° — 3n?)6/n =34(1—1/n). When the 
crystal is infinite this is 346, and a unit 
which requires this work to detach it will 
just be in equilibrium with a solution or 
vapour saturated for a large crystal. 

Now build a new layer on one face. 
The first unit, say at a corner, only has one 
face in contact with the crystal—work 4 . 
The next along an edge has work 2¢, there 
are 2(n—1) of these or 2(m—1) if the 
layer only covers m? of the units, m<xn. 
The remainder all have the full work 3¢. 
This corresponds to what Kossel calls the 
“repeatable step”. Once we have got 
over the initial difficulty, the crystal will 
grow or diminish according as the solution 
is under- or over-saturated for this ‘“re- Figure 1. 
peatable step’. But what is the measure 
of the initial difficulty? The total separation energy for the whole plane is. 
b+4(m—1)b6+(m—1)236=3dm2—2dm. For a line added to an incomplete 
plane the separation energy is 3¢6m—d¢. The defect is due to the fact that we 
have extended two faces of the cube by an amount r2m each, where 7 is the side of 
a unit cube, and added two new steps of equal area, and so produced surface energy 
4mr*a; so 4mr?o =2md or 20 =¢/r?, as is otherwise obvious. 

The defect in separation energy requires an extra concentration of energy in 
the form of a probability fluctuation. he amount required is 2md =4mrc, and 
increases in proportion to the side of the embryo produced. Small fluctuations 
are frequent, but unless the embryo produced is large enough to be in equilibrium 
with the degree of supersaturation present, it will relapse back into its elements and 
come to nothing. ‘The production of an embryo thus becomes more likely the 
greater the supersaturation, because the greater this is the smaller the embryo: 
needs to be. 

While the work for the “repeatable step” is the same for a small or large 
embryo the small embryo is unlikely to survive. It is, indeed, equally likely to gain 
or lose at each step, but its chance of growing really large is small, for it will probably 
be extinguished by chance fluctuations. ‘The case is rather analogous to a man 
gambling at Monte Carlo. If he goes on he is pretty sure to lose all his capital, even 
apart from the zero in favour of the bank, just because of probability fluctuations. 
The smaller energy of attachment of the first few pieces acts against the crystal 
something like the zero against the gambler. ‘The crystal can only be saved by the 
solution being supersaturated. ‘There appears to be no way of supersaturating 
Monte Carlo ! 

A similar explanation may apply to the growth of large crystals at the expense of 
small, as in the process of annealing. 
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In calculating the energy fluctuation necessary to produce an embryo which will 
be in equilibrium with a particular degree of supersaturation account has to be 
taken of the work done by the gaseous (or osmotic) pressure, and this reduces the 
fluctuation required from Fo to F'c/3 in the case of three dimensions (Volmer 1921) 
and from up to up/2 in the case of two dimensions (Volmer and Adhirkari 1925, 
1926), where wu is the perimeter in equilibrium with the supersaturation and p is 
the work of separation per cm? surface. 

Now the trouble is that for moderate supersaturations the energy fluctuation is 
extremely large. Thus Brandes (1927) calculates for rocksalt with 1% super- 
saturation a value of 2:3 x 10-*erg for the surface embryo. This is 6 x 10*RT 
at room temperature and growth should be negligible, while in fact a crystal of 
visible size would grow fairly fast at such asupersaturation. ‘lhe chance of energy 
‘concentration includes of course the Boltzmann factor, exp(— W/RT). 

Amsler (1942) has made experiments on the delay times for the growth of 
crystals from solution, detecting the onset of crystallization by measuring the 
change in the electrical conductivity. The delay time is subject to a probability 
fluctuation, but when the stirring is good the results are surprisingly consistent. 
‘The chance should go as exp (— C/85?) where C = 167M?No?/3p?R?T?, where M is 
molecular weight, p is the density, NV is Avogadro’s constant and 6 the super- 
saturation. ‘The variation with 6 is as expected, but the constant gives a value of o 
ot the order of 1 dyne/cm. Now a is the difference of the surface tensions of the 
solid and liquid if the contact angle is zero. The calculated value for rocksalt is 
107-5 dyne/cm. for a large crystal and that for the solution is 78-8 dyne/cm. giving 
.a difference of about 30 dynes which is cubed in the exponential so that the factor 
here is of the order 303. Amsler thinks that the change in surface tension of the 
-solid with size may account for the difference but it seems unlikely that it would so 
closely approach that of the liquid, especially as results with NaCl and KCl are 
-closely similar. 

In any case, the difficulty remains in the two-dimensional case which Amsler’s 
work merely confirms. 

The way out seems to lie in an early suggestion of Volmer (1921), to whom 
indeed much of the above theory is due. He supposed that the surface of a growing 
crystal is covered with a mobile supersaturated layer. He and Estermann (1921) 
brought forward evidence for the motion of mercury atoms over the surface of a 
growing crystal, showing that the material reaching a face directly from the vapour 
might be many times less than that which actually formed on it. He also (Volmer 
and Weber 1926) found that crystals of benzophenone appeared to lose and gain 
matter by flow over the surface, and measured the diffusion of a verv thin film of 
benzophenone over a glass plate. It should be mentioned perhaps that the 
transference by vapour pressure was appreciable in these latter experiments. 
Other instances of atoms migrating over surfaces are the effects observed by Stern, 
Estermann and Cockcroft with molecular beams, and the technically important 
effects studied by Langmuir and others in which alkali atoms migrate over tungsten. 

If now we suppose that a layer exists on a growing crystal which can diffuse 
rapidly over the surface, but which is much more concentrated than the vapour or 
solution rotind the crystal, we can explain the rapid formation of embryos in such a 
case. The fluctuations of energy needed to form an embryo of given size, and 
‘the size of the embryo needed to produce growth, are now quite different from 
what we have calculated above. The equilibrium condition is enormously 
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hastened, as by a catalyst. There is strong evidence for the existence of such a 
layer, not only in the work of Volmer, but in a line of experiments started by Berg 
(1938) and recently developed at Imperial College by Humphreys-Owen, to whom 
I am indebted for allowing me to give some of his results. The experiments 
consist in illuminating a crystal growing from solution between optical flats in such 
a way as to produce interference fringes. From these the variation of concen- 
tration from point to point of the solution can be determined, lines of equal 
concentration can be drawn and the flow normal to the force of the crystal deter- 
mined, since this by the law of diffusion is proportional to the normal concentration 
gradient. Berg and Humphreys-Owen agree in finding that the surface of the 
growing crystal is far from being a surface of constant concentration (see figure 2), 
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and also that the normal flow is not constant along a growing face, being most at the 
centre and least at the corners. Since the face in fact grows flat, there must be a 
flow of material from the centre outwards in a layer too thin to be observable, 
which means less than about 0:01mm. The distribution of flow is sufficiently 
often symmetrical about the middle point for one to regard this as the normal case, 
but there are frequent cases of exception duc presumably to greater ease of growth 
at one side cr the other (figure 3). Humpreys-Owen finds that the maximum 
flow in the layer is proportional to the inflow in the normal case, and is about 7% 
of it. 

At this point I should like to show you a film kindly lent me by Mr. Ferguson of 
Imperial Chemical Industries which shows some experiments on crystal growth : 
the film has been shown before but the results have not yet been published. You 
will notice that the growing face increases in thickness discontinuously. Steps 


figure refers, unlike that in figure 2, is growing nor- 
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are formed, often at a corner, which sweep across the crystal. Marcellin (1918) 
seems to have been the first to observe this phenomenon. Sometimes one step: 
will catch up on another and the whole move forward together. Isuggest that we 
see here the absorption of the mobile layer at the step, and that the speed with which 
the step advances is governed at least in part by the mobility of the layer. ‘Thus in 
figure 4 material from the old face of the crystal will be flowing towards the 
advancing step to build it up, while a new layer will be reforming on top of the step 
and flowing in the opposite direction. Both these layers will flow in a way 
determined by some law of diffusion, though it seems reasonable to suppose that 
their thickness is limited to some saturation value. If the step is shallow it will be 
able to advance more rapidly than if it is deep, for it will not have to wait for so much 
material to flow in from parts of the mobile layer ahead of it. This flow can only 
occur by a species of surface diffusion, depending on changes in the concentration 
of the mobile layer. It ought to be possible from Humphreys-Owen’s results to 
deduce a good deal about the properties of the layer, such as its coefficient of 
diffusion. Experiments of this kind should be combined with experiments of the 
Ferguson type under comparable conditions. 

Kowarski (1935), who was one of the first to observe the stepwise growth of 
crystals—in his case paratolindine from the vapour—observed a curious effect, 
which he attributes, I think rightly, to the Volmer layer. He found that, by 
adjustment of temperature, liquid drops could be formed on the crystal surface. 
In certain cases when the crystal was growing in area—the crystals were in the 
form of thin flakes—these drops were drawn to the advancing edge of the crystal, 
and especially to the foremost corner, where they were held as though by an elastic 
force. This seems to indicate a flow of the mobile layer to the edge of the face 
which it is prolonging, carrying the drop withit. All this fits well with the Volmer 
picture. Kowarski also studied the rate of surface growth in different directions 
and showed that it was inexplicable unless the vapour could move over the surface 
from where it first struck. He finds a maximum velocity for the layer of the order 
2.0r-3 «10-7 em)seéc: 

In considering crystal growth one’s interest naturally centres on the atomic 
processes, but there are other factors which may be more important in practical 
cases—for example diffusion. Ifacrystalis growing fast it will exhaust the material 
in its immediate neighbourhood, and any protrusion will be more favourably 
situated than the rest. We have now seen how up to a point a crystal can continue 
to grow flat faces though different parts of them are at different concentrations, 
but there must be a limit to this, and when it comes the crystal grows out to un- 
exploited regions. ‘This leads to dendritic growth. Even for a crystal of compact 
shape the effect is important and may govern the rate of growth. We can make a 
crude model (Valeton 1923, 1924) which brings out the main points, by supposing 
that there is a discontinuity of concentration C,’—C, at the surface of the crystal,. 
supposed a sphere of radius a, and that the flow of material into the crystal is 
R(C,’—C,) per‘cm?. ‘Then if D be the coefficient of diffusion (taking 1 cm? of 
solid crystallite as unit of material) it is easy to show that 


da 7 (CogalC. 
dt alD+1]k’ 


where C,, is the concentration of the distant solution. The two terms in the 
denominator can be regarded as series resistances, one due to diffusion in the bulk 
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of the liquid and the other due to the surface layer. If the resistance due to diffu- 
sion is the important one, da?/dt is constant, but if that due to the surface layer is 
predominant, then da/dt is constant. 

Humphreys-Owen finds that for his crystals the first was the case (his crystals 
were in two dimensions which gives a somewhat different law, but the difference is 
slight). It should be noticed, however, that the relative importance of the terms 
depends on a, but his crystals were fairly small, approximately 2 mm., and the law 
should hold a fortiori for larger ones. 

It is well known that the habit of crystals is influenced by the presence of 
colloids in the solution, which favour dendritic growth. They probably act by 
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Figure 5. Showing normal growth of a crystal for which the 
square of the characteristic dimension increases linearly 
with time. The zero of time is the moment when ob- 
servation started, not that when the crystal started growth. 


reducing the coefficient of diffusion of the salt and so increasing the difficulty of 
supply of the material, making this relatively more important than the resistance to 
embryo formation. 

There may, however, sometimes be another effect which, indeed, almost 
certainly occurs when some dyes are present in solution. ‘These are preferentially 
absorbed on certain faces, and so inhibit their rate of growth. Itis usually supposed 
that the faces finally produced are those which grow most slowly, the more quickly 
growing faces extinguishing themselves. ‘This is not self evident, as is seen by a 
comparison of figure 6(a) and (b). It is conceivable that in certain circumstances 
a face might extend as it grows as in (b). This may occur in certain cases of very 
rapid growth where the supply of material is the determining factor to the exclusion 
of considerations of energy. 
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Another way in which non-atomic considerations may modify crystal growth 
is the purely geometrical effect of surroundings. ‘Thus crystals which tend to 
form flat plates and are grown from vapour on a solid substratum tend to be orien- 
tated with the plates parallel to the substrate (Thomson and Cochrane 1938). 
If the vapour is directed as a molecular jet they tend to be orientated normal to the 
jet, and when both causes are present a compromise is reached. These effects may 
be explained by the greater ease of growth afforded to a plate which happens to 
start parallel to the solid surface, or in the second case by the greater mass inter- 
cepted by the face normal to the stream. In both cases we must suppose that the 
molecules that strike the surface form a mobile layer over 1t which moves so as to 
build up the edge of the discs, as in Volmer and Estermann’s experiment. 

Professor Finch has informed me of a similar effect when crystals are grown in 
fine scratches. 

It is pretty general experience that very thin films formed by spluttering or 
evaporation are, if not amorphous, at least micro-crystalline (Thomson and 
Cochrane 1938). The work of Andrade and Martindale (1935) on gold is a good 
example of this, as is also Lotmer’s (1945) on antimony. Increased deposition or 
heating leads to crystallization. Andrade and Martindale found that the first stage 
was to form spherulites. 

When growth occurs it seems to do so over the surface of the substrate, and the 
crystals show strong orientation with regard to this surface (Thomson and 
Cochrane 1938). It has been reported (Kowarski 1935) that one crystal can 
sometimes influence the orientation of another at a distance, which suggests a very 
thin mobile crystalline layer. ‘The evidence for this 1s not yet conclusive, but it 
would be intriguing if the Volmer layer turned out to have an ordered structure of 
some kind. Other cases of crystals increasing in size with time are probably 
examples of the tendency for the large to swallow up the small. 

So far we have considered a crystal growing from a seed of its own composition, 
or from none atall. Very interesting effects occur when one crystal is grown upon 
another. ‘The classical case is the growth of sodium nitrate on calcite. The two 
are isomorphous and the sodium nitrate grows from solution with its axes parallel 
to those of calcite. Many other cases are known, for growth both from solution 
and from melts, for electrolytically deposited crystals and for crystals grown from 
vapour. 

Isomorphism is not essential. Sodium nitrate will form orientated crystals 
from the melt on mica (Lassen 1934), and very large crystals can be grown in this 
way, but though mica is monoclinic, the top layer of atoms is pseudo-hexagonal 
and similar to the (111) face of soldium nitrate. 

The view is often expressed that orientated growth of this kind, or epitaxis as it 
is called, depends essentially on a near equality of spacing between atoms in the top 
layer of the substrate and those in the layer of the crystal that grows to join it. 
Certainly such an agreement exists in many of the most striking cases, but I would 
like to suggest that the importance of this has been exaggerated and that other 
considerations are equally important. In particular I think that the existence of 
steps and cracks in the mother crystal are often the determining factor. Atoms 
align themselves along these steps, which become a zone axis of the daughter 
crystal, some prominent face grows out over the surface of the mother crystal, 
and thus orientation is established. Certainly this is more easily done if there is a 
similarity in the atomic patterns of the two layers, but even this is not essential. 
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Andrade and Martindale have shown the importance of cracks in determining 
the siting of crystals of silver and gold grown on glass and diamonds. Strictly, 
this is not quite the same thing, since they were concerned with the way in which 
minute crystals scattered over a substrate surface grouped themselves while we 
are interested in the direction of the axes of the individual crystals. Still, the 
attraction of the cracks for the atoms that are trying to form up as a crystal is 
strikingly shown in their work. 

I should like to call your attention to three or four examples in which epitaxis 
has been observed for groups of crystals. 

Let us take first the classical work of Royer (1928) on the growth of alkali halides 
on mica. He found that among 25 alkali halides, both of the face- and body-centred 
types, 11 showed good orientation on mica. In all cases the (111) face of the 
halide coincided with the cleavage of the mica. The pattern on this face is of 
course one of equilateral triangles, and in all cases but one these were so turned asto 
fit the corresponding pattern of the mica. Ifthe difference in the dimensions of the 
pattern was less than 12° orientation occurred, if more it did not. In one case, 
KCl, orientation sometimes occurred with a 14-5°% discrepancy. But there is 
another type of orientation in which the triangles are set at right angles. The fit 
here is much better for KCl—almost perfect—but it applies only to Lof the atoms. 
It is the less common form of orientation of the two. Four other of the halides 
could adopt this form of orientation with less percentage discrepancy than in the 
one they actually chose, but this is not observed. ‘This goes to show that the 
specific nature of the planes and lines on which the fit occurs is more important 
than the actual numerical closeness of agreement. 

Another example is given by metals deposited by spluttering on rocksalt. 

_ They have been studied by Lassen and Briick (Brtick and Lassen 1935, Briick 1936) 
using the method of electron diffraction. ‘They find many examples where the 
orientation actually chosen 1s not that which gives the closest fit, as can be seen from 
table 1. Briick interpreted his results in terms of the mean distance of atom in 
the mother and daughter crystals averaged over a cell, which he found tended to be 
aminimum, but in view of the different electrical character of the Na and Cl atoms 
in the surface it seems doubtful if much weight can be attached to this. 

In all Briick’s cases important lines of atoms coincided in the two crystals, 
and this appears also from the experiments of Riidiger (1937), who studied the 
deposit of Au, Ag, Pt on calcite, mica and fluorite. 

It should however be pointed out that in Briick’s case, at least, there is no 
single direction on the mother crystal, which is always a zone axis for the daughter. 

When metals are deposited on metals there is a strong tendency to follow 
completely the structure of the underlying metal as shown in the work of Cochrane 
(1936) and Finch and Sun (1936), who in both cases deposited the metal electro- 
lytically and tested the result by electron diffraction. Many suitable metals have 
nearly the same structure and spacing, so on almost any theory one would expect 
them to fit in this way. Interesting exceptions are iron on palladium and gold, 
where Finch and Sun found cube planes in contact, but the face diagonal of the iron 
parallel to the cube edge of the substrate. ‘This gives a good fit. ‘The same 
authors found a peculiar effect of nickel on gold in which the same arrangement 
as that for iron occurs, though complete parallelism would actually give a better fit. 

Cochrane found that a substrate of copper forces nickel to adopt the copper 
spacing, which is 2:5%%, different, for very thin deposits. ‘The strain causes 
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frequent twinning of the nickel when the deposit grows thicker. Chromium om 
copper showed two orientations in both of which important zone axes are parallel 
in the two crystals. 

An example of the influence of faults on crystal growth is shown in a paper by 
Willems (1943) which studies the growth of urotropine on gypsum. ‘The small 
parallelograms of the urotropine appear lined up with their sides parallel to steps in 
the gypsum. 

While mere agreement in spacing is not conclusive, even when accompanied, 
as it must be, by similarity in arrangement of the atoms, there is no doubt from 
Royer’s work that it plays an important part. If we take a figure of 12% as giving 
the maximum permissible disagreement for halides of mica, we can estimate the 
size of the embryo crystal which is sufficient to ensure indefinite growth under the 
conditions of the experiment. If the spacings of halide and mica are 1 and1+«, 
then after 7 steps the atoms will be «m apart. We may suppose that while this is 


Table 3. Metallic deposits on copper. Cochrane (1936) 


% error | Any % error 
Metal Orientation w.r.t. better Orientation Wickets Remarks 
Cu fit ? Cu 
Takes Cu lattice) 
Ni Completely || — 24% No _ — constants for 
thin layer 
Cu Ditto ’ Nil No = ya an 
Zn Does not orientate — —= = == pa 
Cd Ditto == = = aS om 
: (111)’s |]. 
Ag Completely || +13-2% 1 Yes Rotated 30° w.r.t. | —2:0% Ne oe 
one another. 
(110) Cr || (111) Cu. } +12-7% Tae 
Cr Cube edge of Cr || Plane _ “a 1 poe 
face diag. of Cu — 81% or orientation 
Co Completely || — 16%] No ae a 


less than 3, the forces on the whole tend to favour growth compared with a random 
position, but when it exceeds this value then they are unfavourable. Hence 
an=%, which for «=0-12 givesn=4. The spacing in question being 5-15 a., we 
get 20 a. as the size of the embryo. 

One might expect the necessary size of the embryo would decrease with in- 
creasing supersaturation, and it would be interesting to try if any effect of this kind 
can be observed. 

Another kind of orientated growth is that which sometimes occurs when a 
crystal takes part in a chemical reaction. Many years ago I was able to show 
(Thomson 1930) that a single crystal of copper could oxidize to Cu,O in such a way 
that the axes of the two crystals, both of which are cubic, were parallel. Since 
then a number of other examples have been found. In studying these we must 
remember that only certain faces of the original crystal will be exposed, and that, 
except in the case of cleavage faces, these are not necessarily the apparent faces, 
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since, for example, etching may produce pits covering the whole or almost the 
whole of the apparent surface. It is only on exposed faces that orientated growth 
can occur. 

An examination of the observed instances allows certain general conclusions to 
be drawn. In general we have todo with a reaction in which new atoms are 
introduced, as when metals are oxidized, or in which some new atoms are sub- 
stituted for old ones, as in the growth of one halide of silver on another, or silver 
bromide on thallous bromide, recently studied by Schwab (1947). 

In most cases the atoms that remain keep their relative positions, though 
the distances between them necessarily often change. In some cases, e.g. growth 
of cuprous bromide and iodide on the (111) face of copper (Usmani 1941), the 
arrangement of atoms is similar in the two crystals but of different orientation. 
Here the copper atoms in both the metal and the halide forma lattice of equilateral 
triangles, but the two are turned with respect to one another through a right angle. 
The two planes appear always to have the same symmetry and general arrangement 
of atoms. 

There is a preference for a fit in actual spacing where this is possible. ‘Thus 
the rotation of the planes of the copper halides just mentioned greatly improves 
the fit, which is about 3% for the iodide and still better for the bromide, though it 
necessarily implies that not every atom has a partner in the new crystal. 

The requirements of fit are not too exacting. Thus Schwab (1947) found 
hexagonal AgI can grow on the (111) face of AgBr with a difference of 13°, and 
some of the oxides show even larger differences, 18°/, for copper and 28°, for 
barium (Burgers and Van Amstel 1936). Schwab has pointed out that the growth 
of AgBr on T1Br is probably helped by the fit normal to the plane of contact, 
which is (110) in both crystals, the large bromine atoms filling spaces left by the 
smaller metal atoms in the layer. The discrepancy in dimensions in this case is 
2/5: 

There is much yet to be found out in the study of orientated growth. It is 
probable that different considerations are dominant in different cases. Equality 
of spacing is important but is not the only thing that matters. Crystal growth, and 
especially the part played in it by the Volmer layer, is a fascinating subject, which 
may be expected to increase greatly our knowledge of the other surface pheno- 
mena. ‘lhe study of orientated crystals produced by chemical action is one of the 
most promising ways of finding out more about the mechanism of chemical change 
on which, after centuries of chemical study, science is still grossly ignorant. 
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“Kikuchi-Line’’ Patterns— 
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ABSTRACT. Seven methods are described of indexing electron-diffraction Kikuchi-line 
patterns from single crystals for the purpose of practical applications of such patterns. 
The crystal lattice form and orientation can be determined by method 1 (see Part I), which 
also gives the indices of the diffractions relative to the chosen unit cell. ‘The other six 
methods are in general only applicable in conjunction with a knowledge of the crystal lattice 
form and orientation, in which case any one of them, or more conveniently several in con- 
junction, can be used to index the lines. 


$i. INTRODUCINION 

HE only generally applicable direct method, independent of information 
| about the crystal other than that contained in the diffraction pattern and a 
knowledge of the diffraction conditions, is the method described by Wilman 
(1948)—to be referred to as Part I—involving a construction of the reciprocal 
Jattice. When a suitable unit cell in this reciprocal lattice has been chosen, the 
crystal lattice orientation and Laue symmetry can be defined, and the Laue indices 
of the Kikuchi-line diffractions follow simply as the coordinates of the correspond- 

ing reciprocal-lattice points, relative to the chosen reciprocal axes. 

‘The other methods described here may be useful, for example, when too small a 
range of diffracted beams is recorded, or is distinguishable, to apply the absolute 
method described in Part I, or when the crystal nature and structure are already 
known but it is required to explore the intensities of a large range of diffractions, or 
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to define the orientation more precisely than can be done from the spot pattern 
when the spots are too diffuse or too few to apply the methods previously developed 
by the author (Charlesby, Finch and Wilman 1939). 

One or other of these methods, or several in conjunction, can generally be 
applied successfully, so that the desired information about the relative orientations 
and relative reflecting powers of the different net planes can be derived and specified. 


§2. OUTLINE OF THE MAIN FEATURES OF THE METHODS, 
AND SOME ESSENTIAL COMMON CONSIDERATIONS 


(1) The general characteristics of the methods 


‘The seven methods listed below are of three types. The first two methods are 
characterized by construction of a projection in which each net plane is repre- 
sented by a point. It is only the first method, the reciprocal-lattice projection, 
which incorporates information on the net-plane spacings as well as on the direc- 
tions of the net-plane normals represented in the second, gnomonic projection, 
method. ‘The main feature of the second type of method is the derivation of the 
ratios of the net-plane or the Laue indices from the ratios of the intercepts cut off 
by the net-plane trace, i.e., the line-pair median, either on the three lines which are 
along the projections of the crystal axes on the plate (method 3), or (methods 4 and 
5) on two perpendicular axes on the plate, one of which is along a definite (/,,/,) 
net-plane trace on the plate. The third type of method uses either the zone 
relationships (method 6) or the closely related parabolic and circular envelope 
relationships (method 7). All the methods are developed for the general case of a 
triclinic crystal, and the expressions applicable to crystals of other systems can 

be derived by giving a, b, c, «, 8, y appropriate values. 

Methods 2, 4, 5 and 6 supply only the Miller indices of the reflecting net plane 
as a first step, and the order number of the Bragg reflection must then be found 
from the calculated plane spacing d and the measured line-pair separation D. 
Methods 1, 3 and 7, however, provide the Laue indices of the line diffractions 
directly, the measurements of D being used in combination with the other data 
which define the net-plane direction. 

The first method alone actually carries out a construction of the reciprocal 
lattice and thus an independent determination of the crystal lattice, assuming only a 
knowledge of the camera length L and the position of O,, the foot of the perpen- 
dicular from the crystal origin O to the photographic plate. The wavelength 
must also be known if an estimation of the axial lengths a, b, c is required instead 
of just the axial ratios. Method 2 can also be used to a limited extent for a 
determination of the lattice form and orientation, but only in rather special cases 
where the crystal setting is a favourable one. 

In ‘all the other methods the further knowledge is assumed of the lattice 
constants a, b, c, x, 8, y and of the crystal orientation relative to the photographic 
plate or the primary beam. This information, if not definitely known or obtained 
by the method of Part I, can sometimes be inferred from the relative positions of 
the main bands in the pattern, possibly together with information provided by the 
geometrical form of the spot pattern and, if necessary, in conjunction with an 
examination of other patterns obtained at different settings, or patterns of other 
types such as powder and rotation patterns. 

The various methods are mostly best suited to different spheres of application. 
Method 1 is the most fundamental, generally applicable and independent. 
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Method 2 is not usually so convenient as the other methods and is only readily 
usable when the polar direction of the projection is along or close to a crystal axis. 
or a principal lattice row. Method 4 (Finch, Quarrell and Wilman 1935) was. 
developed by the author to study the relationships between the Kikuchi lines and. 
their parabolic and circular envelopes in reflection patterns at different azimuths ;, 
and method 5, also evolved in this laboratory (Thirsk and Whitmore 1940) is a 
form of method 4 more directly expressed for application to particular main 
azimuths. The neater and simpler method 3 was developed by the author in 1939 
for easier practical use, especially when the crystal orientation is not very near to a 
main setting but is known relatively closely from the Laue spot positions as in the 
case of one of the anthracene patterns reproduced by Charlesby, Finch and Wilman 
(1939), where this method was first applied. Methods 6 and 7 are easy to apply 
and are often exceedingly useful, especially if only the intensified line of a pair is 
distinguishable ; and when desirable these methods can be systematically applied 
in the way shown below. 


(ii) Preliminary considerations applicable to all methods 


The following considerations should be noted in indexing a pattern, whatever 
the method used. 


(a) Conventional viewpoint to be used in indexing prints. 'Vhe indexing of the 
pattern can be carried out on either a negative or a positive print; the latter is 
preferable for clearness. A positive print should be considered as corresponding 
to the electron distribution in the recording plane seen from the side of this plane 
opposite to that on which the electrons fall and on which the crystal lies (Part I). 
A negative print (or the initial ‘‘negative”’ recording observed from the sensitive 
emulsion side) should be regarded as corresponding to the electron distribution in 


the recording plane, viewed from the side on which the electrons fall and on which 
the crystal lies. , 


(b) The Bragg-reflection positions at the band edges. ‘These have been shown 
in Part I to be at, or practically at, the points where the intensity is equal to the 
background intensity in the region considered. 


(c) The black and white line positions. The positions of the Bragg reflection 
loci and the intersections of the net planes with the plate are determined by the 
orientation of the crystal relative to the plate, and are best specified relative to the 
foot O, of the normal from the crystal origin O to the plate. The primary electron 
beam OF is not in general along OQ,, and of a Kikuchi-line pair the one farthest 
from OF is the intensified line and the one nearest OF the absorption line relative 
to the general background intensity. 


(d) Convention for the sign of the Laue indices of the line diffractions. The 
direction of a net plane is specified equally well by the Miller indices (hkl) or 
(—h, —k, —1), but ina Laue spot pattern the diffraction directions represented by 
the Laue indices hkl and —h, —k, —Jare different, and the sign of each index is 
fixed by the directions of this diffraction and of the primary beam relative to the 
corresponding crystallographic axis. In Kikuchi-line patterns, owing to the 
partial reciprocal relationships between the diffraction lines of a pair, there remains 
the question of which line of the pair is to be called the Aki diffraction and which the 
—h, —k, —l. Itis important to adopt a convention to specify the sign to be used 
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for each line, not merely for uniformity, but in order to be able to use the zone and 
envelope Laue-index relationships (methods 6 and 7) without errors arising from 
this cause. 

A convention which is easy to apply, and which is now recommended, is the 
following : Consider the two lines of a pair as the intersection with the plate of 
the two halves of a cone whose apex is at the crystal origin O, each line of the 
pair corresponding to one half of the cone on one side only of the net plane (also 
through O) whose trace on the plate is the median of the line pair. Then, having 
chosen the positive directions of the axes a, b, c, the positive (hkl) plane normal is 
that from the origin O to the plane which cuts off intercepts a/h, b/k, c/l on a, b, ¢ 
respectively, and the Kikuchi line which corresponds to the half-cone lying on the 
same side of O as this plane can be assigned the Laue indices AA/ or the appropriate 
multiples nh, nk, nl. Conversely, the other line of the pair has Laue indices —nh, 
—nk, —nland corresponds to the half-cone lying on the opposite side of O, where 
also is the plane whose interceptson a, b, c are —a/h, —b/k, —c/l respectively. 
A paper or card model comprising the ad, bc, and ca faces of the unit cell is helpful 
in visualizing orientations of a particular net plane when the crystal orientation 
relative to the photographic plate is known. 

The above convention is automatically followed if the indexing of the lines is 
done by method 1, but it has to be specifically applied in all the other methods of 
indexing. 


(e) Conversion of plane-trace and Kikuchi-line positions to those which would be 
obtained on a different recording plane. 'Uhe expressions derived below in methods 
3, 4, 5 for the positions of the lines and the line-pair medians or plane normals take 
more convenient simple forms in the case where the known direction which defines 
the normal to the photographic plate has relatively small integral] indices. If the 
plate normal OO, in any particular case does not approximately coincide, as it 1s. 
usually arranged to do, with such 
a lattice row, the position that a 
line or a band median would have 
on a plate normal to any desired 
lattice row [uvw] inclined at 6 to 
the actual plate normal can be 
constructed as follows (figure 1). 

Let O, be the (observed) inter- 
section of the lattice row [uvw] 
with the plate and O,M be normal 
to O,O,. Let QM be a line-pair 
median, i.e., a net-plane trace on 
the plate, produced to meet O,.M foe 
at Mand O,0, at Q: Then, as 
can be seen from the side elevation in figure 1, the line OQ which lies in the net-- 
plane under consideration is prolonged to meet at Q’ the plane which is normal to 
OO,, i.e., [uvw]. If this new recording plane is now rotated about the line O,M 
so as to become coincident with the actual plate, O,Q’ rotates to the position 
O,Q"; whence the position of Q” is as shown in the plan part of figure 1. In the 
triangle OO,Q, let angle O,OO, =5, angle 0,0Q =é, and O,Q=y, then 


y/sin € = L’/sin {180° —€ —(90° —8)}, 
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and Q” is given by 
Q"0,/Q0,=(y +Ay)/y =(L'|y) tan €=c088/ {1 —(y/L') sin 8}. 

The plane trace on the new recording plane is then the line Q’M, which 
becomes Q”M when rotated into the plane of the plate. It follows that any other 
ordinate RN is also changed to RN where R"N/RN = Q’O,/QOs. 

As a particular example it can be seen that if 8 is less than about 5° and 
L'~20cm. with RN <6cm., then Ay/y<2%, thus under these conditions, to 
a first approximation, the pattern can be indexed as if OO, was the normal to 
the plate. 

The new Kikuchi-line positions can be found, if necessary, by treating as 
plane appropriate small parts of the corresponding diffraction cones, and applying 
the above construction. 


§3. THE METHODS OF INDEXING THE PATTERNS 


(i) The reciprocal-lattice projection method 


This has been fully described and illustrated in Part I and will not be discussed 
further here. 


(ii) The gnomonic projection method 


A gnomonic projection can be constructed very easily from the positions of 
the band medians. If any net plane is shown in side view by the line OP in figure 2 


5 


Plane of 
Projection 


Figure 2. Gnomonic projection. Figure 3.  Stereographic projection. 


(Kx, and K, representing a Kikuchi-line pair from this plane), then the plane normal 

is OP’, O,P’=Lcot¢=L?/p; ¢ is the inclination of the net plane to OO,. The 
gnomonic projection on the recording plane comprises points such as P’ where the 
plane normals meet this plane, and it is thus simply obtained by marking the 
position of O, on the pattern and then marking a point corresponding to each band 
median, at a distance proportional to 1/p from O, ona line through O, normal to the 
band median, and on the opposite side of O, from the band median. The appli- 
cation to x-ray diffraction patterns has been illustrated by Wyckoff (1931), 
Bragg (1933), and others. _ 

Although stereographic projections are not easy to apply to the indexing of 
Kikuchi-line patterns, it may sometimes be desired to construct such a projection, 
and this can be done as follows. Representing a net plane in edge view by OP 
‘(figure 3) and its normal by OP’, meeting the plate at P’, its stereographic projection 
on a plane through O parallel to the plate is P’, of radius R. By joining the point 
where OP’ cuts the sphere to the opposite pole S$, the point P” is found as the 
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intersection of this line with the equatorial projection plane. It is easy to show 
that the length OP’, ie. p”, is 


p” = Reosd/(1 +sing) =R[{1 + (p/L)2}! +p/L}4. 
To construct a’ stereographic projection from a Kikuchi line pattern it is thus 


necessary to measure off from QO, (or a point displaced suitably outside the pattern 


area) distances p” in a direction opposite to that of the normal drawn from O, to 
each line-pair median. 


(ui) The method of intercepts on the three axial projections on the plate 


The Laue indices of the Kikuchi lines can be calculated directly from the 
bandwidth D and the intercepts which the band median cuts off on the projections 
a,, b,,, c, of the crystal axes a, b, c, drawn from the point O, where a known lattice 
row [wvw] meets the plate 
normally (or nearly so). In 
figure 4, let O be the origin 
in the crystal and OO,, of 
length LZ, be the normal from 


O to the plate. Let a, b, ce ‘ 
mectathe plate ratw A. B.C, ue iver if 
respectively, so that O,A, vi é / 


O,B, C,C are in the directions i 
of a,, b,, ¢, and let the (AAl) / 
plane through O cut O,A, Photographic Plate Trace of (kL) Plane 


O7B) O, eat xX YoaZ. distant 
x,y, zfromO,. Also let N be 
the unit vector along OO, and n that along the normal from O outwards to the 
plane whose intercepts on a, b, ¢ are a/h, b/k, c/l; the angles between N and 
a,b, care 0, 9,, 0, respectively. Then we have 


O0,A= 0A — O00, =(L/cos 6,,)(a/a) — LN. 
The equation of the line O,A relative to O as origin is thus 
r=LN+mL{(a/acos@,)—N}, 


where mis a parameter. Since the point X on this line lies in the (AA/) plane 
through O, its radius vector r must also satisfy r.n=0. 
Hence, with a. n=hd, 


IN .n+mL[{(a.n)/acos6,}-—N.n]=0, 


1.€; mt = COS 4,,/{cos Oy zi. (hd/a COS Ba) 5s 


Figure 4. 


where @,, is the angle between Nand n. We now have, directly, O;A=Ltan8,, 
whence O,X =x+m0,A =Ltan 6, cos @,/{cos 0, —(hd/acos 0,)} 
and h=(a/d) cos 8, . [cos 6, —(L/x) sin 4,], 


and there are similar expressions for & and 1. 

The values of d and cos6,, must be expressed in terms of the bandwidth D, 
the Laue indices AAl of the diffraction line, and the length p of the perpendicular 
from O, to the trace XYZ of the (hk/) plane (subtending 90°—4@,, at O), thus 


d=Aj(2sin0)~AL/Dsin#,; cos6,=p/(L?+p?)!; sind, = L/(L* +p”). 
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Hence h=[pD/(L? + p?)|(a/A)[cos 0, —(L/x) sin 64]. 
Similarly k=[pD/(L? + p?)](6/A)[cos 6, — (L/y) sin 44], 
l=[pD/(L? + p?)|(c/A)[cos 0, — (L/z) sin 6,]. 

In order to calculate h, k and J for a line pair it is thus necessary to calculate @,, 
9, and 6, and the directions of a,,, b,, ¢, for the known or estimated crystal setting 
(actually in most main settings these are almost obvious or very simply calculated), 
and to measure AL, «, y, 2, p and D and then use these equations. A need not be 
known if only h :k : lis accepted as sufficient to determine h,kand/. ‘The lengths 
x, Y, 8 are positive or negative according to whether they are measured from O; 
along the “positive” or the ‘“‘negative”’ direction of a,, 6, or c,. The signs of 
h, kand also depend on p and @,,._ If p is always positive and thus 6, always an 
acute angle, i.e. to the positive direction of the net-plane normal 7 always down 
from O towards the plate, then with the index convention stated in § 2 the indices 
h, k, l are assigned to the line diffraction (black or white) which lies on the side of 
the plane trace XYZ nearest to O,, and the line on the opposite side XYZ has the 
indices with all three signs reversed. 

When only one line of the pair is visible, so that p and D cannot be measured, the 
bandwidth D can be estimated roughly from the sharpness of the line, or much 
more exactly from its position in relation to other lines which form with the line 
concerned a parabolic envelope with the associated zone axis intersection at its _ 
focus, so that D/2 can be measured. 

If p=0, and thus x = y=z=0, the only finite terms remaining in the equations 
for h, k and /are those containing the factors p/x, p/yv and p/z which are the cosines 
of the angles between O,P and O,A, O,B and O,C respectively. 

Calculation of 6, 64, 0, when OO, is along [uvw]. If the angles 6,, 0, 6, are not 
directly known, but a lattice row [uvw] is known to be along OQ,, the values of 
cos @,,, sin @,, etc. can be calculated as follows for insertion in the above equations 
for h, k and 1. 

If T is the (positive) [wvz] lattice translation, 


cos 6,=T.a/Ta=(ua+vbcosy +wecos f)/T, 
cos 6, =T . b/ Tb =(uacos y + vb + wecos«)/T, 
cos 6,=T.¢/Tc=(uacos B+ vbcosa+we)/T, 
where T? = u?a? + vb? + wc? + 2uvab cos y + 2vwhe cos x + 2wucacos p. 

Since sin? 6, = 1—cos?6,, this gives 
sin? 6, =[v°b? sin? y + wc? sin? B + 2vwhe (cos « — cos B cos y)]/T?, 
sin? 6, = [wa" sin? y + wc? sin? « + 2uwac (cos B — cos y cos «)]/T?, 
sin? 6, =[w°a? sin? B + vb? sin? « + 2uvab (cos y — cos « cos )]/T?. 

Thus, for example, 

h=(pD/(L* + p*);(a/A)(1/T) . [(ua + vb cos y + we cos B) 
— (L/x){v°b? sin? y + wc? sin? B + 2vebe (cos « — cos B cos y)}3]. 

These expressions take much simpler forms in the more symmetrical crystal 

systems, thus for example for cubic crystals, 
h=[pD|(L2-+p®)\(aA)(u2-+ 02 + 2%)-!fu —(L/x)(02 + 0°)4, 
k= [pD|(L2+p?)[(a/A)(u2 +02 + 20%)![o —(L)y)(u? + 0°), 
1=[pD|(L? +p2)|(a)d)(u2 + & +202) [oo—(L]2)(u? +2°)!]. 
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Calculation of the angles between a p» b,, and ¢, and the normal to the trace of the 
(Ayk4l,) plane passing through [uvw], normal to the plate. When the directions of the 
projections a, Dic rare not conveniently determined by inspection from the 
known orientation or the geometry of the spot pattern, i.e. the Laue zones along 
the spot rows, they can be calculated from the known indices [uv] of OO, and 
those (h,k,/,) of a plane which passes through this zone axis and whose Kikuchi 
lines are identified in the pattern. 

Thus, the vector r,* =h,a* +h, b* + i,c* lies in the plate, normal to (h,k,/,), 
and since a, =a—(a.N)N 


») 
cos (a,,r,*) =(a,. r,*)/a,7,* =(d,/a,){a.r,* —(a. N)(N.r,*)} =h,d,/asin8 


cos (b,,r,*)=,d,/bsin6,, and C0s(C), fs )=1,a7/cs1n 0, 


a) 


where sin 6,, sin 6, and sin@, are given by the expression stated above, in terms of 
u,vandw. For cubic crystals, for example, cos (ap, 01*) =hy/(hy? + k,2 +1,2)* sin 8,, 
UC. 

It is desirable to check the relative positive directions of a,,, b,,, c,, carefully, 


preferably with the aid of a rough model, before proceeding to index the pattern. 


(iv) The method of intercepts on two perpendicular axes along and perpendicular to the 
trace of a plane (hyk,1,), in terms of the azimuth ¢ in this plane 


An analytical derivation of the equation of a net-plane trace or Kikuchi line 
has been given and illustrated previously (Finch, Quarrell and Wilman 1935), and 
it can also be obtained in a similar way by vector methods. ‘The results are as 
follows. 

First imagine the crystal unit cell to be mounted in a fixed, but at present 
unspecified, position relative to three orthogonal reference axes Ox OV Ole 
(A convenient setting is that shown in figure 5.) Now transfer the cell, together 


Z 
A 12,-(4X1) 
; plane normal 


Figure 5. Figure 6. 


with these axes, to the diffracting crystal, and orientate it so that the axes a, b,c 
of the unit cell are parallel to the corresponding axes of the crystal, whose origin is 
at O (figure 6). Let OO, be the normal to the fluorescent screen or photographic 
plate and assume the known crystal orientation to be defined by the indices 
(hyky/,) of a net-plane whose trace on the plate is the line O,'X, and by the azimuth 
angle ¢ between two planes normal to (/,4,/,), one also normal to the plate and thus 
containing OO, and OO,’, and the other containing OY, and intersecting (h,h,/,) 
in OY. 
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The orthogonal set of axes OX,, OY,, OZ, is further defined by setting OZ, 
along the positive normal to the net plane (/,4,/,), i.e. from O to the plane which cuts 
off intercepts a/h,, b/k,, ¢/1, on a, b, ¢ respectively. A further orthogonal set of 
axes OX,, OY;, OZ, is such that OZ, coincides with OZ,, and OX, is parallel to the 
plate and lies in the (#,4y/,) plane through O, i.e. OX,Y,Y,. Finally the orthogonal 
axes OX, OY, OZ are such that OX is along OX, and OY is along O,O; the angle 
YOY, is ws. 

On this basis, if the normals to the net planes (,,/,) and any other plane (A2/) 
have direction cosines A,, 44, v, and A, p, v relative to OX;Y3Zs5, the equation of 
the line where the plane (AA/) through O intersects the plate is, relative to O,;X 
and O,Z on the plate (figure 6), 

Cyz= —A,x+ BL, 


where 
ere ht A,=A,cos¢—B;sin¢, 
B, = B, cosy Cysin ys, B,=B,sin¢+B, cos 4, 
C=C, cosy — B, sin, C=C 
A, = (Avy — vAy)/(1 — 27), 
Bs = {pu — py (Ady + eye + 791) }/(1 — 2?) 
Cy=My + pop + ry, 
in which 


1=\24 24+ 7=APt+p?t+vP=A?+ B+ CP=A?+ B+ C,2=A,?4+B,? + Cs. 


If the origin of coordinates on the plate is transferred from O, to O,’, which is 
on the (A,k,/,) plane trace where it cuts O,Z, the (AR/) plane trace has the equation 
C,z= —A,x + B,L/cos %, and if is small this is approximately C,z= — A,x + B,L, 
which is the same as the equatign for the plane trace relative to O,XZ with ys put 
equal to zero. 

The equation of the mth order Kikuchi-line pair from the (AR) plane is, relative 
to O,X and O,Z on the plate, the hyperbola 


(Aye + Bry + Cys)? =(a2 +92 + L2)(md/2d), 


and when x? + y? is small compared with L? this approximates closely to the pair of 
parallel straight lines, Cyze= —A,x+{B,L—(nAL/2d)}. The equation relative to 
O,/ is obtained analogously from the above equation for the plane trace, and when 
y is small it is effectively Cz = —A,w+ {BL + (mAL/2d)}. 

If the crystal setting and lattice constants are known, these expressions can be 
evaluated, for example as shown below. Since the equation of the (ARI) plane 
trace involves h, k and / only linearly, the two intercepts of any line-pair median on 
O,X and O,Z (or one intercept and the slope) suffice to allow the ratios h :k :1 to 
be determined. ‘The angle ¢ is taken here as the angle through which the crystal 
must be rotated in a positive direction (right-handed screw) about the positive 
direction of the (/,,/,) plane normal, to bring the crystal from the setting where 
OY, is in the plane YOZ, into the actual crystal orientation. 

The axes OX;Y3Y3 can be chosen conveniently relative to a, b, c as shown in 
figure 4, such that OY; is along 6 and OX,is in the ad plane on the same side of b as 
is a. In this case, for the general triclinic crystal, 


A=[(h/a)+(Rcosy)/6]d/siny; e=kd/b; 
v= {(h/a)(cos « cos y — cos B) + (k/b)(cos B cos y— cos «) + (I/c) sin? y}abed/V sin y, 
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where d=(V/abc)/{X(h?/a?) sin? « + 2X (hk/ab)(cos « cos B — cos y)}! 


and V = abc(1 —cos? x — cos? B — cos? y +2 cos « cos B cos y)}. 


There are similar expressions for Aj, p14, v, in terms of hy, ky, L,. 

The angle ¢ is not in general known directly as part of the specification of the 
crystal orientation, but it can be calculated for example from the indices [w,v,w,] 
of the lattice row which is along OO,’ in the (h,k,/,) plane, since then we have 


cos = (cos 8)/44 = (ua cosy +v,b +. wc cos a)/ 47, 


where 6 is the angle between [u,v,w] and OY,, i.e. b, and T, is the [u,v,2,] lattice 
translation given by 7? =Xu,?a? +2Eu,v,ab cosy. 

Thus in any particular case, knowing the crystal lattice constants, the orientation 
defined by the indices (h,ky/,) and [u,7,z,] and the angle ¢ which can be measured 
from the pattern in terms of the position of the (4,k,/,) line-pair median relative to 
O,, the procedure to index the pattern is as follows. First write down the values 
of Aj, #4, ¥, and A, 4, v from the above expressions by inserting the known indices 
hy, ky, 1, and a, b, c,d, V,«, B, y, and calculate cos ¢ and sin ¢, as just described, from 
[w,2,w,]. ‘Then putting these values in the expressions for Ay, By, C, etc., the 
equation for the (Ak/) line median is obtained, involving A, k, J linearly. Finally 
the position of each line-pair median in the pattern can be compared with this 
equation and two relationships between h, k and / thereby established which give 
h :k <1; from this the simplest integral value of h, k and / consistent with the band- 
width D must be deduced, bearing in mind the accuracy of the measurements. 

This method, like the next, is unwieldy for crystals which have low symmetry, 
but it has been easily applied to the patterns from some cubic and hexagonal 
crystals (Finch, Quarrell and Wilman 1935, Finch and Wilman 1936, 1937 a, b). 
For cubic crystals, for example, the expressions take the simple general forms 


A=hdla,  p=kdja, valdja, = d=a/(#2? + +2)}, 
M=ahda, py=kd,la, v,=hdJa, dy=al(he+k2+1,), 
A, = (hl, —1h,)|(h? +h? +2)(hy? +1), 
B,=[k—ky(hh, + kk, +U,)]/[(2R2)(Zh2\(h2 41,2), 
Cy = (hh, + RR, +1l,)/[(Zh2)(Zh2)]3, 
cos $6 =[v,/(Xm,?)[Eh2/(h.? +1,2)]}. 


The general analytical expressions derived above are also useful in studying the 
positional relationships between systems of net-plane traces or Kikuchi lines, such 
as the line envelopes. 


(v) The method of intercepts on two perpendicular axes along and perpendicular to the 
trace of an (h,k,l,) plane in which [u,v,w,] 1s normal to the plate 


Thirsk and Whitmore’s (1940) expressions are given below, in a notation 
uniform with the above, and the equation of a line pair is added. 

Figure 7 shows the disposition of the reference axes OX, OY and of the known 
plane (/,k,/,) normal to the plate, with positive direction of the known lattice row 
[u,v 2] along the normal OO, from the crystal origin to the plate. The trace MN 


PROC. PHYS. SOC. LXI, 5 28 


426 H. Wilman 


of any plane (hkl) then has, relative to OX and OY, the equation Px + Ky + L=0, 
where 


u,a2+v,abcosy u,abcosy+v,b? u,accosB 


A +w,cacos B +w,becosa  +4,becosa+wyc? 
= é! 
V(mh+vk+w,l) h, ky l, 
h k l 


\(74,1,) 
(AkYnormal plane trace 


ee (nhl, plane 


(Akl 
ae nN 


x 


— 


L eye) 0 
Z(h, k,l) plane trace 

Figure 7. 

Note: vertical through O should read “‘ plane normal’’. 


R={T,d,(a?b?c?/V?)|(ujh+v ,k+w,l)} 
x [Z(Ah,/a?) sin? « + X{(hR, +h,k)/ab}(cos « cos B —cos y)], 
L=camera length, . 
V =volume of cell = abc(1 — cos? « — cos? 8 — cos? y +2 cos « cos B cos y)?, 
d,=(h,k,l,) spacing =(V/abc)/{X(h,?2/a”) sin? « +22(h,R,/ab)(cos« cos B — cos y)}, 
T,=[u,v1w,] lattice translation = (ua? + 24u,v,ab cos y)?. 


Thirsk and Whitmore give, as an example, the expressions for the case of a 
hexagonal crystal; in that for R, however, their vc)? in the numerator should read 
Uv Ca 


The equation of the mth order Kikuchi-line pair from (ARI) is 
Pr+ Ry t+ Lt [(nAL/2dyy)(l +(x? +y°)/L? PT |(Suyh) diag]: 


(vi) The zone relation method 


The zone relationships. If a net-plane (hl) is parallel to two lattice rows 


[w,22,] and [vz] whose lattice translations are T, and T,, then 7*,,, is perpendi- 
cular to 7, and T,; thus 


wht+ukt+wl=0; uh+v,k+wl=0, 


which gives A : Rk :1=(v 2 — V1) :(WyUy — Wey) : (UyVg—UgV,). Conversely, the 
zone axis [uvw] parallel to two net-planes (/,k,,) and (Ayk,l,) is defined by 
U:v :W=(Ryl,— Rol) 2 (yhg—Ighy) : (hyRe —hehky). 

‘This zone relation method seems to have been the sole method applied before 
1935 to index the line-pair medians in Kikuchi-line patterns, for example by 
Kikuchi (1928), Shinohara (1932), Kirchner (1932) and Shinohara and Matukawa 
(1933). Nevertheless, no general systematic method has hitherto been stated for 
indexing the various zone-axes corresponding to the points of intersection of 
line-pair medians on the plate, and the method has only been applied with relatively 
simple crystal orientations where an easy geometrical construction suffices to 


identify the main zone axes. General systematic methods of calculation are now 
outlined below. 
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Calculation of the position of the intersection of a zone-axis [uvw] with the plate. 
Any lattice row (uvw] along OP (figure 8), with lattice translation T, makes with 
the plate normal OO, along the reciprocal-lattice direction [Aokyly]*, i.e. r5*, 
an angle 7 given by cos 7 =(T . ry*)/T79* = (uh+vk+wl)d,/T. The radial distance 
of P from O, in the plate is then O,P = L tan n. If the direction O,X of the trace 


(u vu)* 
reciprocal 
lattice plane 


tO (2oKolo) 


Figure 8. Figure 9. 


of a known plane (/,k,/,) through OO, is identified, and PQ is the normal from P 
‘to O,X, then O,P makes with O,X an angle €, given by 


sin £, =QP/O,P =(OP . r,*)/7,*O,P = (wh, + vk, + wh)d,/Tsin n. 


‘The distance O,P and the angle & thus specify the point P where [uvw] meets the 
plate. In the above expressions, d; ~AL/D,, where D, is the measured hykyl, 
line-pair separation near O,. From the above equations it follows also that 
‘y, = QP =0,P sin & = L{(auh,)/Xuhy\d,/d,. 


Calculation of the indices [wvw] of a zone axis from the position of its intersection P 
with the plate. From the above equations, conversely, 


Uhy + vk + wly=(T/do) cosy; uh, + vky + wl, =(T/d,) sin & sin; 
but to determine [uz] uniquely it is necessary to measure also the angle €, made 
by O,P with the trace of another net-plane, for example (/,k/,), through O,, whence 
Uhy + UR, + wl, =(T/d2) sin & sin 7. 


mhen w:v:w= 
(cos 7)/do, ko, I hg, (cosn)/do, I hy, Ro, (cosy) /dy 
| (sin€,;siny)/d,, ky, 4): |Ay, (sing,sinn)/d,, 4]: |hy, ky, (sin€,sin n)[dy| . 
|(sin€,siny)/d2, ke, ly hs, (sin€,siny)/d,, 1 hy, ky, (sin €,sinn)/d, 
Alternatively, from the above expression for y, (and y,), we have the following 
two equations, if py=d,/dy and p,=d,/d): 
uthoy, —hyLpy} + v{RoV1 — ky Lpy} + {lov —1,L py} =0, 
UN V2 — aL pa} + V{RoV2 — ka L ps} + w{lyy2—1,Lp.} =0, 
which gives w:v:w as in the case of a pair of zone-relation equations given above. 
28-2 
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Calculation of [hpRolo|* of this direction is initially specified by [woo]. If OO, | 
is initially known to be the [uvpwo] lattice direction, then the indices [HoRalo| ot 
this direction in the reciprocal lattice, required in the above calculations, may be 
obtained readily by comparison of the two equations 

Tj? =i ae en cee 
T, =ua+tupb + wyc=(Tp.a)a* + (Ty. b)b* + (Ty. c)e*. 
Thus, hy =Ty. a=A(tly a+ U9) COS y + Wy COs f); 
ky = Ty. b=)(uy acosy + Ub + Wy COS &) ; 


1,=T,).c=c(uy acos B +p b cosa +2 Cc). 


(vii) The envelope-relation method 


All the net planes such as (hk/) which are parallel to the same zone axis [wvw] 
intersect the plate in lines which are along the corresponding Kikuchi-line pair 
medians, and they all pass through one point, where the line [wow] drawn through 
the crystal origin O meets the plate. The indices (hkl) of all planes of the zone 
satisfy the relation uh + vk + wl=0, but the Kikuchi lines of the group of cozonal 
planes fall into subgroups, the lines of each subgroup being tangential to a common 
parabolic envelope. 

The envelope curves themselves are not in general present in the pattern as 
true parabolic curves (black or white), though in the main subgroups where a large 
number of lines touch the envelope locus at frequent intervals, the polygonal 
nature of the apparent curves is only recognizable if the crystal is highly perfect 
and if a narrow well-focused monochromatic electron beam is used to obtain high 
resolution (Finch, Quarrell and Wilman 1935, Shinohara 1935, Finch and Wilman 
1937b). ‘There is usually noticeable a characteristic weakening and curving away 
of Kikuchi lines just before they touch the envelope, also contributing to the 
prominence and appearance of continuity of the main envelope loci. — 

As Shinohara (1935) has shown, the lines associated with any one parabola arise 
from net planes such as (hkl) belonging to a zone [uvw]; the reciprocal-lattice 
points associated with these planes not only lie on the (wvw)* plane through the 
origin O, fulfilling the condition wh + vk + wl =0, but also lie on a line, for example 
in an [/,k,/,]* direction, lying in the (uww)* plane parallel to the plate or nearly so, J 
but not passing through the origin. (If it did, the lattice points along it would only jf) 
correspond to the orders of diffraction from a single net-plane.) If this reciprocal-_ |] 
lattice line [h,k,/,]* is densely populated, the approximation to the parabolic |) 
envelope curve will be prominent. It can be shown that, if the zone axis is normal 
to the plate, or nearly so, the axis of the parabola is parallel to the perpendicular ON 
to this reciprocal-lattice line [h,k,/,]*, and is thus, in the plate, along the trace of the 
net-plane (/,k,/,) (figure 9). The Laue indices of the successive Kikuchi lines 
touching the parabola differ by the same integers h,k,/, as the indices of this net- 
plane, or by a multiple of them. 

Thus, if the indices of two or more lines touching a particular parabolic locus 
are known, those of the plane (,k,/,) along its axis are the differences of these two ]) 
sets or a sub-multiple of these differences ; the indices of all the other lines touching | : 
the locus are also easily found even if only one line of each pair is distinguishable, 
and can be checked by the position or slope of the line (see methods 3, 4, ©) 
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Conversely, if the plane (/,k,l,) whose trace is along the parabola axis has been 
indexed and the indices ki of any line touching the parabolic locus are known, 
then those of the other lines are of the form h+nh,, k+nk,, l+nl,. 

If the sign convention of §2 is adopted in assigning the Laue indices Ak/ to any 
Kikuchi line, black or white, the above reciprocal-lattice interpretation shows 
(figure 9) that Kikuchi lines touching a parabolic locus on one side of its axial net- 
plane progressively further from its vertex have indices obtained from those of the 
previous one, by adding the Laue indices (or a multiple) of the first-order Kikuchi 
line due to the axial plane of the parabola and lying on the same side of the net-plane 
trace, i.e. the parabola axis. 

As has been shown (Finch, Quarrell and Wilman 1935), the parabolic envelopes 
in turn fall into groups whose envelopes are a system of circles with centres at the 
points where main lattice rows, if produced, would meet the plate. Shinohara 
(1935) has also shown that the circular envelopes represent the limiting loci 
outside which lie all Kikuchi lines characterized by the same Laue index h’ relative 
to a lattice row meeting the plate at the centre of the circle, the lines all lying on the 
side of the net-plane traces, opposite to that on which the centre of the circle lies. 
The angular radii of such circles are given by cos 6=1—(h'A/T), i.e. 6 = (2h'A/T), 
where T is the lattice translation of the lattice row directed towards their centre 
and h’ is any integer; thus the circles have radii nearly proportional to V/h’, 
and have the same positionsas the h’th-order Laue zones of reinforcement from the 
lattice row corresponding to an imaginary primary beam parallel to the row. 

If the Laue indices of a Kikuchi line are Aki relative to the axes a, b, c, and the 
lattice row is [wvzw], the line belongs to the group associated with the h’th-order 
circle round the [uvzw] trace on the plate, where h’ is the Laue index of the diffraction 
direction Al relative to the lattice row [uvw], that is hu+kvu+lw=h'. Many of 
the lines of a group will not touch the circle nor even lie close to it, but they fall into 
subgroups as above, touching parabolic envelopes which do touch the circle; 
the connection between any particular line and one of these circles is thereby 
established, and if [www] and hf’ are known, the above linear relation between h, k 
and / is obtained. A line may be seen to belong to two or more non-coaxial 
parabolae or circles, and in this case several relations such as the above hold 
simultaneously and allow hkl to be determined without difficulty, on this basis 
alone, without reference to the bandwidth. 
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ABSTRACT. The object of the research was to examine the diffuse pattern obtained. 
when an electron beam is transmitted through a crystal of anthracene. A special camera 
adapted for low temperature work is described. It has been found that the diffuse pattern. 
has a fine structure which is more apparent in the case of thick crystals. ‘The changes 
in intensity of the diffuse pattern with temperature have been established experimentally 
and it is shown that the phenomenon of the diffuse scattering is analogous to that found. 
in X rays. 


§1. INTRODUCTION 

LECTRON diffraction transmission photographs from single crystals of 
anthracene show, besides a cross-grating pattern, a superimposed pattern. 
which appears to be that of diffuse bands and areas. These bands and 
areas are related to the orientation of the anthracene molecule within the crystal 
lattice and to the separation of carbon atoms in the anthracene molecule; this. 
led Charlesby, Finch and Wilman (1939) to the assumption that the diffuse 

areas are due to what they call “‘ molecular scattering’’. 5 
The following is abrief summary of their theory. The molecules are 
vibrating independently of each other in the crystal lattice and their vibrations. 
have no definite phase relationship with each other. Although built in the 
crystal lattice they are producing an effect similar to that of a gaseous mass of 
oriented molecules in the electron beam. _ In the consideration of the vibrations. 
the intramolecular heat motion is neglected. Molecules are thus vibrating 
as rigid units about their mean positions in the crystal lattice. As the anthracene 
unit cell contains two molecules in different orientations the authors proceed to. 
construct qualitatively a molecular scattering pattern such as would be given 
by two independent molecules in these orientations. This diffuse pattern is. 
to be superimposed on the cross-grating pattern of the crystal. It is claimed 
that the calculated diffuse pattern agrees well with the observed positions of 
diffuse areas and bands and gives the correct intensity distribution. Furthermore, 


* Part of Thesis submitted to London University for Ph.D. degree examination; 
{ Now at Imperial Chemical Industries, Billingham. 
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they have derived an expression for the dependence on temperature of the 
intensity of the molecular pattern, on the assumption that the thermal displace- 
ment of any molecule from its mean position is given by the Maxwell-Boltzmann 
law. 

No confirmation for either the theory or the experiments of Charlesby, 
Finch and Wilman has been found in recent x-ray work. In this field another 
temperature-dependent phenomenon was found. It was examined experimentally 
by Laval, Preston, Lonsdale and Smith, and others, and it received theoretical 
treatment by Faxen, Waller, Born, Zachariasen et al. (see review by Lonsdale 
1943). ‘The experimental researches of Lonsdale and Smith (1941) were 
carried out on organic as well as inorganic crystals. No X-ray investigations 
have been reported on anthracene, but benzil, sorbic acid, hexamethyl-benzene, 
urea nitrate and urea oxalate have been studied. A very brief outline of the work 
of Lonsdale and Smith will be given here. 

A simple theory of diffraction of monochromatic x-rays by a crystal lattice 
predicts that with most crystal orientations (where the Ewald sphere of reflection 
misses all the reciprocal lattice points) no trace of any pattern should be found on 
photographs, no matter how long the crystal is exposed to x-rays. Experimentally, 
however, it was found that, with long exposures or very powerful x-ray sources, 
some definite pattern of diffuse spots of various shapes, and streaks of different 
types (sometimes connecting the diffuse spots), could be seen on the photographic 
plate. The effect was specially clear when highly monochromatized radiation 
was used. ‘The intensity of the diffuse pattern was greatly diminished when the 
specimen was cooled with liquid air. The temperature change was reversible. 

The diffuse pattern varied considerably from substance to substance and 
showed very interesting features for some chain and layer structures—for chain 
structures a characteristic sharp streak and for layer structures a large diffuse 
spot. ‘The diffuse streaks connecting diffuse spots are found to lie along or nearly 
along some reciprocal lattice axes. Generally the intensity of the diffuse spot 
is largest near the position of a Bragg spot, and although the intensity diminishes 
with crystal rotation from this position, diffuse reflections could sometimes 
be observed when the crystal was turned as much as 30°. The intensity of a 
diffuse spot is related to the intensity of the corresponding Bragg reflection; 
thus no diffuse spots were observed corresponding to reflections forbidden by 
structure-factor considerations. 

The phenomenon of diffuse scattering is accounted for by the Faxen-Waller 
theory which not only predicts its existence but also gives the formulae for the 
shape of the diffuse spots and accounts for the dependence of their intensity on 
temperature. ‘The Faxen-Waller theory assumes that atoms are vibrating 
thermally in the crystal lattice and that the vibrations are harmonic. As the 
maximum frequency of atomic vibrations (~101 sec.) is much lower than 
that of x-rays (10'* sec.-') the atoms or molecules in a crystal will appear 
stationary as far as the passage of x-rays is concerned but slightly displaced from 
their mean positions. It is possible to consider the atomic displacement as due 
to superimposed, standing, elastic waves in the crystal which can be of acoustical 
or optical type. Although the complete theory must take into account both 
acoustical and optical branches of the vibrational spectrum, nevertheless the 
acoustical, which corresponds to larger amplitudes, is of primary importance. 
The theory shows that a single stationary wave of this type would produce two 


A32 F. Skrebowski 


extra points very near each main reciprocal lattice point. ‘The shorter the wave- 
length the larger the distance of these extra spots from the main spot. The 
effect of thermal motion can be represented by a cloud of such extra points, 
giving extra reflecting power (or diffuse reflection) to every region surrounding 
a reciprocal lattice point. The shape of this extra reflecting cloud is dependent, 
among other factors, on the elastic constants’ (Lonsdale 1943). So far, the 
Faxen-Waller theory has not been worked out in detail for organic crystals but 
it explains qualitatively the peculiar diffuse reflections from chain-like molecules 
and layer structures and also the diffuse streaks extending between two spots. 

The results predicted by the molecular scattering theory of Charlesby, Finch 
and Wilman, and those predicted by Faxen and Waller are essentially different. 
While, according to the molecular scattering theory, it should only be possible 
to observe relatively large diffuse areas or bands without any associated fine 
structure, the Faxen-Waller theory predicts diffuse spots and streaks. Later, 
Charlesby and Wilman (1942), on re-examining their plates, found evidence 
of diffuse streaks as well as diffuse areas on the electron diffraction pattern 
of anthracene. There is also an attempt by Charlesby (1942) to reconcile 
the molecular scattering theory with that of Faxen-Waller, on theoretical 
. grounds. 

The aim of the present research is to study the diffuse areas observed by 
Charlesby, Finch and Wilman, and to see whether these can be interpreted in 
terms of the Faxen-Waller theory which has been used successfully in the x-ray 
work. 

§2. APPARATUS 

A camera of the type designed by Thomson and Fraser (1930) was modified 
so as to accommodate four photographic plates. This design of camera uses 
a cold-cathode discharge as electron source and a fine beam is obtained by pin- 
hole collimating, no magnetic focusing being employed. The accelerating 
voltage of the electrons was usually in the region of 28 kv. The distance from 
specimen to photographic plate was 42 cm. 

As the experiments were to be carried out at low temperatures, the main 
cone was redesigned to enable the specimen to be cooled down to the temperature 
of liquid air (figure 1). A large brass cone B is bored out. Inside it a tube T, 
about 3 cm. in diameter and 25 cm. in length, made of German silver, is soldered 
at its upper end to a flexible bellows C, which in turn is soldered to the main cone. 
The position of the tube can be adjusted by three screws. The tube T at its 
lower end is closed by a brass disc D. A part of the specimen holder H is tightly 
screwed into D so as to ensure good thermal contact between these two pieces. 
The specimen holder, made of copper, consists of two thick plates H and F, 
both with hemispherical cavities. ‘These are put together with four screws. 
A thick ring with a rod-like extension A, on either side in the plane of the ring, 
is enclosed between the plates and can be rotated about its axis. After H, A 
and F are put together, a short cylinder R can be inserted into the ring from 
outside and, as the cylinder R fits tightly into the ring, good contact is ensured 
between them. The specimen-carrying cylinder R is a piece of thick copper 
tube with a fine wire gauze soldered over the top. On this gauze a transmission 
specimen can be placed. In the main cone there is a small cone M which can 
be rotated about its axis. To its lower end a violin string is attached and also 
fixed to a small wheel W at the end of the axis A. Thus by rotating cone M, A can 


Electron diffraction from anthracene 433 


also be rotated. With these arrangements it is possible to ensure all the necessary 
movements, viz. up and down by three screws, rotational movement of the whole 
cone B, and rotational movement of the specimen-cylinder R. A calibrated 
copper-constantan thermocouple is soldered to H. All the temperature measure- 
ments were made with a potentiometer bridge. 

After the camera was 
evacuated, liquid air was 
introduced into the tube T 
and it generally took 5-10 
minutes to cool the speci- 
men. ‘The amountof liquid 
air necessary to keep the 
specimen at low temper- 
ature is small. In the 
initial experiments at low 
temperature it was found 
that the pattern observed 
on the fluorescent screen 
became blurred after a few 
minutes. It was probably 
due to condensation of 
organic vapour from the 
diffusion pump on cold 
crystal surfaces and it was 
necessary to introduce on 
both sides, H and F of the 
specimen holder, two: 
shields Z in the form of 
tubes about 2-5 cm. long. 
After this modification the 
crystal could be kept for 
a long time without any 
appreciable change in the 
diffraction pattern. 

It should be noted, Figure 1. Specimen holder arrangements for orienting 
however, that as soon as specimen and cooling device. 
the apparatus is cooled, a 
thermal contraction takes place and removes the observed crystal from the beam. 
Although it should be possible with screws to bring the specimen holder back 
to its previous position, in practice it was found impossible to find a very small 
crystal again once it was lost from the beam. A slightly different arrangement 
was therefore tried. ‘T'wo pin holes, S,; and S,, were fitted into the shielding 
tube on the side from which the beam is incident, so that the nearest hole was 
only about 0:5 cm. from the carrier R. The distance between the holes was 
about 3 cm. and their diameter was 0-1 mm. With this arrangement it was possible 
to find the same crystal after lowering its temperature. It also produced a very 
fine beam and reduced the area of the crystal upon which the beam was falling. 
The experiments with this device proved particularly successful in elucidating 
the problem of the fine structure of the diffuse areas. 
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§3. CRYSTALLOGRAPHIC, DATA, PREPARA LION AND MICROSCOPIC 
EXAMINATION OF ANTHRACENE CRYSTALS 

Anthracene (C,,;H,)) has a melting point of 217° c. and sublimes slowly 
in vacuo at room temperature. It belongs to the monoclinic system, the 
elementary cell having the following dimensions :—a=8-58a., b=6-02a., 
c=11-18a., with B=125°. Anthracene crystals occur in the form of hexagonal 
flakes, the plane of the flake being (001); they can easily form twins. ‘The 
unit cell contains two molecules, the atomic coordinates of which were obtained 
by Robertson (1933). As anthracene crystals are often contaminated with. 
naphthacene, very pure material was used in order to eliminate the possibility 
of the presence of foreign molecules which might affect the diffuse areas. It was 
found that this highly purified anthracene formed better crystals than the com- 
mercial material. 

Thin crystals of anthracene suitable for electron diffraction were prepared 
by heating a sample of the pure substance in a test tube. A large number of 
minute crystals showing interference colours was formed from the vapour and 
as they floated in the air they were easily picked up on a fine gauze. Another 
way of preparing suitable crystals was to evaporate xylene solution very slowly 
on a very thin collodion film. Crystals formed by each of these methods were 
examined and photographed under the microscope. ‘Their appearance differed 
according to the method of preparation. 

While crystals grown from xylene solution formed regular flakes lying 
sometimes on top of each other, those grown from vapour often formed clusters, 
twins or intergrowths at different angles. A general view of crystals grown 
from vapour on to a glass slide as they appear under the microscope is shown 
in plate 1. During observation unde+ the microscope the crystals disintegrated 
by developing several holes, although the flake itself was still thick. It seems 
plausible that the electron beam covers portions of the crystal of very different 
thicknesses. 


§4. GENERAL DESCRIPTION OF THE SPOT PATTERN OF 
ELECTRON DIFFRACTION PHOTOGRAPHS 

Most of the electron diffraction photographs obtained by transmission 
through anthracene crystal flakes showed a fairly extensive pattern due to the 
ab cross-grating inclined to the beam at different angles. Sometimes the pattern 
was very extensive. ‘This corresponded to a very thin crystal, possibly only 
a few molecules thick. Other photographs showed only a few intense spots 
and were correspondingly due to much thicker crystals. The pattern on the 
photographs was mainly rectangular, or nearly so. Most photographs, whether 
corresponding to thick or thin crystals showed irregularities in the cross-grating 
array of spots, frequent doubling or tripling of the spots, convergence of some 
rows, and generally untidy and irregular grouping of spots. From a large number 
of photographs taken with crystals grown from vapour, only a few showed a 
perfect cross-grating pattern such as would be expected from a single crystal; 
even on these a very noticeable regular doubling of spots was observed. Crystals 
grown from xylene solution on to collodion did not show these double spots, 
but unfortunately the rings due to the collodion film tended to obscure con- 
siderably the details of the anthracene pattern, especially near the centre. 

It is of interest to note that spots forbidden by the structure factor, of the type 
(2n+1, 0, m), (0, 21 +1, 0), where n and m are integers, appeared on all photo- 
graphs. ‘They were rather faint and quite sharp. 
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§5. PATTERN INDEXING AND ESTIMATION OF CRYSTAL 
AMSOK CARIN SS) 

; The plates were indexed by assuming the values for the cell dimensions 
given by Robertson (1933). Sections of the reciprocal lattice were drawn to 
scale and to allow for the crystal being thin, lines were added to every lattice 
point in a direction perpendicular to the surface of the crystal. This was taken 
to be an (001) face. ‘The sphere of reflection was then orientated to give the 
observed pattern. 

Crystal thickness can be estimated from the length of the extension of the 
reciprocal lattice points. If some spots fail to appear on the photographic 
plate then the reflection sphere misses the corresponding (extended) reciprocal 
lattice point because its extension is too short. Once the length of the extension 
is found the crystal thickness is estimated from the formula: /=1/t, where ft is 
the crystal thickness and / is the length of the extension measured from the centre. 
It may be remarked that the reciprocal lattice points of anthracene are fairly 
close together. Hence, with crystal thickness below 20 a., the reciprocal lattice 
points are so much extended that they form continuous lines. Thus it is not 
possible to estimate crystal thickness below that value. 


NODE SC RUPTON OF VEE Be D THRs 2 PAM ERINIS 
(a) Experimental results at room temperature 


The character of the diffuse background in the diffraction patterns appeared 
at first to depend on the crystal thickness and the plates were accordingly 
classified in this way. ‘lhe deductions below were made from a large number 
of plates of which examples, illustrating the salient points, are reproduced. 

(i) Thick crystals (t>200 a.). Thick crystal patterns are characterized by 
relatively few Laue spots, most of which are intense. ‘hese patterns show, 
on well exposed photographs, an array of diffuse spots and streaks, but nothing 
which can be termed a diffuse area. Plate 2 shows conclusively that the diffuse 
background has a well defined structure, the details of which can be seen in the 
enlarged portion (plate 3). ‘These bear a striking likeness to the diffuse back- 
ground of benzil, in plate 106 of the paper by Lonsdale and Smith (1941). 
To obtain detail, it is particularly important to use plates which show good 
contrast (e.g. photomechanical plates show more detail than process plates). 
It is also necessary to use a fine electron beam to avoid the pattern being a 
composite one of a number of crystals. In these composite patterns the super- 
posed diffuse background may have the appearance of a diffuse area, as the 
detailed structure tends to be blurred. That the diffuse pattern is very faint, 
can be seen by comparing its intensity with that of Laue spots visible in the same 
photograph (plate 3). The dependence of the pattern on temperature will be 
discussed below. 

Another feature which has been noticed on a large number of plates is that 
some of the diffuse spots are elongated and not circular. In special cases these 
are very elongated and are found on both sides of the central beam. Plate 4 
shows a typical example of this. These elongated spots were fairly persistent 
when the orientation of the crystal was changed. 

(ii) Thin crystals (t~100 a.). Thin crystal patterns are characterized by an 
extensive array of cross-grating spots. ‘These spots vary a great deal in intensity 
(roughly in accordance with the structure factor) and forbidden spots are also 
present. A typical example of the thin crystal pattern is shown in plate 5: 
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this shows also the type of background called diffuse areas by Charlesby, Finch 
and Wilman (1939). This is a composite pattern due to a number of crystals 
in more or less the same orientation. Both the print and the plate show some 
signs of fine structure, though it is not so evident as for the thick crystals. 

In order to eliminate one disturbing factor—the composite nature of the 
specimen—a special arrangement described in § 2 was introduced which ensured 
a very fine electron beam. ‘This did not result in perfect single crystal patterns 
but under the best conditions clearer photographs were made possible. Further, 
particular care was taken in the time of exposure. ‘Too short exposure gives 
very little contrast resulting in vague patches of background and a loss of detail. 
Excessive exposure leads to halation and again loss of detail. ‘The choice of the 
right exposure time is thus essential. 

Plates 6 and 8 are patterns of the same crystal and practically the same 
orientation taken at two different times (see below) with a fine beam. In the 
first case, plate 6, the exposure is rather short, and though a fine structure can be 
seen on the original plate, it is not very evident in the reproduction; in plate 8, 
on the other hand, the exposure is stronger and the fine structure of the background 
(for instance the streaking between spots) prominent. 


(b) Experimental results at low temperatures 


The object of the experiments at low temperatures was to find whether the 
diffuse background disappears or weakens at the temperature of liquid oxygen. 

An ideal experiment would be to take a photograph at room temperature 
and then, with unchanged crystal setting and the same beam intensity, to cool 
the crystal down to the temperature of liquid oxygen and take another photograph. 
This was, in fact, done by Lonsdale and Smith with organic crystals using 
x rays. Unfortunately, in electron diffraction with the present experimental 
arrangement, it proved to be impossible, as neither the setting of the crystal 
nor the beam intensity could be strictly controlled. The procedure adopted 
was to take a photograph at room temperature, cool the specimen-holder to the 
temperature of liquid oxygen, find the same crystal and take several photographs 
with different exposures but using the same voltage. It was difficult to find 
two photographs which were exactly comparable; very slight movement of a 
specimen could change the spot pattern, though the new pattern could be identified 
as being due to the same crystal with slightly different orientation or a slightly 
different thickness. 

Plate 6 taken at room temperature shows diffuse patches which, on the original 
plate, have a fine structure consisting of diffuse streaks joining cross-grating 
spots. ‘The pattern indicates that the crystal is thin. Plate 7 shows the same 
crystal cooled to the temperature of liquid oxygen. Although the exposure is 
stronger than that of plate 6 very little diffuse streaking is visible on the original 
plate and the strength of the whole background is considerably reduced. The 
cross-grating spots are also stronger than in plate 6. Plate 8 shows the same 
crystal warmed to room temperature again; the background is more prominent 
than in plate 7, and the fine structure more obvious than in plate 6. These 
three photographs show that the temperature phenomenon is _ reversible. 


A number of sets of photographs of the above type were taken, all of which showed 
the same features. 
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§7. DISCUSSION 


The existence of diffuse spots and streaks, observed with anthracene crystals, 
suggests that the effect is the same as that found by Lonsdale and Smith (1941) 
in their investigation on the x-ray diffuse scattering from some organic crystals. 
As mentioned above, the x-ray diffuse pattern of benzil (Lonsdale and Smith 
1941, plate 10 6) bears a striking resemblance to the pattern shown in plate 3. 
In other photographs, e.g. plate 4, the characteristic elongated streaks on either 
side of the central beam, and their observed persistence with tilt of the crystal 
are similar to the streaks observed in x rays with sorbic acid (Lonsdale and Smith 
1941, plates 12 6 and 21 a, b, c). 

The temperature effect cannot be demonstrated as clearly as with x rays 
because of considerable experimental difficulties, but the examination of several 
sets of photographs, e.g. plates 6, 7, 8, shows conclusively that the diffuse pattern 
intensity diminishes considerably at the temperature of liquid oxygen. This 
decrease is especially well marked in the case of diffuse streaks joining the diffuse 
spots; the streaks disappear almost completely at low temperature. ‘Thus it 
can be concluded that the effect of diffuse scattering as found in x rays has a 
close analogy in electron diffraction. 

The diffuse scattering in x rays was interpreted successfully in terms of the 
crystal vibrations and it is reasonable to conclude that a similar phenomenon 
in electron diffraction is due to the same cause. It is, therefore, unnecessary 
to assume that the diffuse spots and streaks are due to the molecules vibrating 
independently, which has been postulated by Charlesby, Finch and Wilman 
(1939); indeed, the fine structure observed, plates 2, 3, 4 cannot be reconciled with 
their theory. 

Emphasis has been placed by Charlesby, Finch and Wilman (1939) on the 
fact that the diffuse patches observed on their patterns from anthracene crystals 
are due to independent scattering from anthracene molecules, and that, therefore, 
the molecular structure of anthracene can be deduced from the relative positions 
of the diffuse areas on the patterns. While it is possible that some molecules 
at the surfaces and edges of the solid may scatter independently, the molecules 
forming the lattice do not. 

Examination of the patterns obtained in the present work showed that in 
general, the diffuse background was strongest in those positions where Laue spot 
intensities were strongest; this is to be expected on the Faxen-Waller theory. 
It can be shown on this theory (Born 1942-43) that for any crystal the intensity 
of the diffuse scattering is directly proportional to the square of the lattice structure 
factor, in the region near the Laue spot. One would therefore expect that 
an examination of the structure factor, F, would show where the diffuse back- 
ground would be most intense. Contour lines of F’, for anthracene, corre- 
sponding to the orientation in plate 5, have been drawn, and are shown in 
figure 2. ‘The positions of the strong regions coincide reasonably well. 

On the molecular scattering theory, the intense regions of the background 
should be given by each molecule scattering independently. This results in 
a slightly different set of contour lines from the previous set (see figure 3) but they 
are not significantly better as a rough representation of the intensity distribution 
found on the plate. 

The orientation chosen in plate 5 is a fair one from the point of view of the 
molecular scattering theory, since the cross-grating spots are grouped, as regards 
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intensity, into only three regions. Neither set of contour lines (figures 2 and 3) 
shows any detailed structure over a region comparable with the cross-grating 
unit. A careful examination of plate 5, however, shows that there is a streaking 
between spots and a fine structure. The molecular scattering theory fails, 
therefore, to provide a basis for a full description of the experimental results. 
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F Figure 2. Contour lines of intensity scattered Figure 3. Contour lines of intensity scattered 
from a unit cell of anthracene. Beam from two independent molecules oriented 
direction as in plate 5. as in the anthracene cell. Beam direction 


as in plate 5. 
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DISCUSSION 


Dr. K. Lonspate. Mr. Skrebowski’s paper is of great interest, in that it suggests 
‘that experimental data on diffuse scattering, which appeared to have been satisfactorily 
accounted for by the assumption of independently vibrating molecules, can be explained 
equally well by assuming that the diffuse intensity is dependent upon the ordinary structure 
factor, being greatest in the neighbourhood of reciprocal lattice points with the largest 
reflecting power. That this is true in general is quite clear to anyone who has worked in 
this field; it may be seen illustrated in plates 5 and 6 of Lonsdale and Smith, Proc. 
Roy. Soc. A, 1941, 179, 8, for x ray diffuse scattering; although of course the rule cannot 
be applied in detail for all structures. 

It is interesting also that the photograph (author’s plate 7) taken at liquid air tempera- 
tures shows very clearly not only the decrease of intensity of the diffuse pattern but also 
the increase of intensity of the Laue spots (especially of those at large diffraction angles), 
which is a necessary corollary. This fact has misled Dr. Wilman into thinking that plate 7 
must have been taken using a thinner part of the crystal than that which gave plates 6 and 8. 

Another point which comes out more clearly perhaps from the electron diffraction 
photographs than from x rays (although quantitative ionization spectrometer measurements 
given by Laval for KCl, and Geiger counter measurements on urea and on oxalic acid 
dihydrate shortly to be published in the Acta Crystallographica, do prove this fact without 
any doubt) is that the falling off of diffuse intensity with increasing angle of diffraction 
is markedly less than that of the Laue intensities. This also follows from the Faxen-Waller 
theory, as extended by Born and others. 

I see no particular significance, however, so far as this controversy on interpretation 
is concerned, in the presence or absence of diffuse ‘“‘ areas ”’, meaning by these the overall 
blackening of considerable regions of the photographic plate. The practical applications 
of the Faxon-Waller theory so far made have only considered the parts of reciprocal space 
in the neighbourhood of reciprocal lattice points, but in all but the simplest cases (and 
anthracene is not a simple case) it is quite certain that optical waves will have to be taken into 
account, and that the effect of these will be greatest between the reciprocal lattice points. 
In one sense the effect of these optical vibrations is likely to be much more like that of 
independent molecular vibrations than the acoustical vibration effects would be, but I am 
not clear that the present method really provides a crucial experiment by means of which 
a distinction can be made. ‘The presence or absence of diffuse scattering in the neighbour- 
hood of reciprocal lattice points of zero Laue intensity is not such a criterion, for any 
vibration or distortion which disturbs the crystal lattice will necessarily also disturb the 
conditions which cause the Laue reflection to disappear. The surprising thing, perhaps, 
is that such “‘ forbidden ”’ reflections do not appear more often in diffuse form. Some 
unpublished work by Mr. P. G. Owston on the diffuse pattern of ice shows. that such 
forbidden reflections do appear in that case, where the intermolecular forces are much 
stronger than they are in anthracene, and independent molecular vibrations are very 
unlikely. ; 

My opinion is that the electron diffraction photographs of thick anthracene crystals 
are the most useful for comparison with x-ray diffraction effects, because the use of very 
thin crystals will necessarily introduce a possible confusion between surface and internal 
structure, and because with very thin crystals some of the long wavelength standing waves 
will no longer exist. 


Professor Sir GEorGE THOMSON. ‘This is a paper which raises an interesting theoretical 
issue. It is one of those cases in which two widely different theories lead to nearly the same 
answer. Born’s theory has been well confirmed for x rays by the work of Dr. Lonsdale and 
others, and there can be no doubt of its essential truth. Unfortunately it is very complicated, 
and involves a number of elastic constants. For substances as complicated and crystal 
forms as unsymmetrical as anthracene, it is not possible to work it out in detail. One is 
therefore driven to qualitative predictions, as, for example, that the diffuse effects are most 
marked near points of the reciprocal lattice with large structure factor. ‘Tt is inherently 
extremely probable that the Born effects will also occur with electron diffraction, and I 
am convinced that the effects first observed by Finch and Wilman and studied in this paper 
are the electron analogue of those studied for x rays by Dr. Lonsdale. The general agree- 
ment in appearance, the close correlation with structure factor, the character of the fine 
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structure observed by Dr. Skrebowski, all confirm this. But that is not to say that the 
theory put forward by Charlesby, Finch and Wilman, which ascribes the effect to diffraction 
from molecules vibrating independently under thermal agitation, is necessarily wrong. 
The Born theory is a theory, and I believe a complete one, of the effect of thermal agitation 
on the scattering of x rays. It is expressed in terms of certain elastic constants. These 
elastic constants are determined by the structure of the crystal and the forces between 
the atoms in it. If certain groups of atoms are held together very strongly so as to form 
molecules, and if the forces between molecules are relatively weak, these facts will be 
reflected in the values of the elastic constants, and by the time you have gone through all 
the mathematics it is probable you will get roughly the same answer as if you had started 
by assuming the molecules to be independently moving rigid solids as Charlesby, Finch 
and Wilman have done. Only you will not get the fine structure which only the more 
refined theory is capable of predicting. I regard the independent rigid motion of the 
molecules as inherently probable though not certain; a flat molecule like anthracene might 
easily bend, but I do not think that experiments of this kind are capable of proving it, 
since there is so little difference, at least in the cases tried so far, between the distribution 
of intensity deduced from the structure factor and that deduced from the hypothesis of 
independent molecular motion. 


Dr. A. CuarLesBy. The problem with which Mr. Skrebowski deals, i.e. the presence 
of diffuse bands in electron diffraction, is of considerable value in the determination of 
molecular structure. Under certain conditions it enables the arrangement of atoms in 
the molecule to be determined directly, without going through an elaborate Fourier analysis 
of Laue spot intensities. Several papers have shown their occurrence and their correlation 
with molecular structure. 

Mr. Skrebowski attempts to explain away the existence of these diffuse bands by 
ascribing them to the blurring or smearing out of relatively sharp streaks. These sharp 
streaks are, however, quite distinct and have indeed been previously observed in electron 
diffraction and x-ray work. I fail to see how photographic plates can smear out a streak 
perhaps 1 mm. wide into a band 5 mm. in width, while yet reproducing sharp spots less 
than 0-1 mm. across. 

Using Process plates, the author was unable to obtain diffuse bands, but observed sharp. 
streaks. This does not disprove the existence of these diffuse bands. His observation 
arises simply from the sharp contrast and limited intensity range of these plates, which sup- 
press low intensities. It is for this reason that after experimenting, we preferred Special 
Rapid to Process plates. Obviously it is not a sufficient explanation to show that under 
unsuitable conditions of exposure a phenomenon does not appear. It is far more important 
to explain why it does appear under other conditions. 

I would like to refer to another question raised in the paper, namely the comparison 
of diffuse streaks and diffuse bands. On the basis of rather superficial similarities, mainly 
that both phenomena are more diffuse than the spot pattern, Mr. Skrebowski draws the 
conclusion that bands and streaks are identical, whereas a closer examination reveals wide 
differences between them. 

My conclusion is that we are dealing with three distinct sets of phenomena: (a) sharp 
spots corresponding to the crystal structure; (b) diffuse streaks and spots, observed both 
in x-ray and electron diffraction work; these streaks arise from thermal motion and corre- 
spond to elastic waves throughout the crystal; (c) diffuse bands, so far only observed in 
electron diffraction, whose non-appearance in present-day x-ray work is readily explained 
(Charlesby, Proc. Phys. Soc., 1942, 54, 379). ; 

In the case of a few simple structures, the effect on the diffraction pattern of elastic 
vibrations throughout the crystal may be evaluated. These vibrations give rise to streaks 
and spots, the presence of which is confirmed by x rays.. In the case of more complex 
crystals the mathematical difficulties appear insuperable, and one is obliged to proceed 
by qualitative comparisons. Similar streaks and diffuse spots have been observed in 
electron diffraction patterns and reported in Nature. It may be presumed that they arise 
from the same continuous type of thermal vibration. 

Instead of considering these continuous thermal vibrations throughout the crystal 
one can go to the opposite extreme, and consider each molecule as vibrating independently 
of its neighbours. On this alternative assumption, the mathematical calculation can be 
readily completed even for a complex crystal such as anthracene. There is very good 
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Plate. 1. Anthracene crystals grown from vapour. Magnified x 50. 


Plate 4. Thick crystal pattern of anthracene showing fine structure 
and two elongated diffuse streaks. 
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Plate 6. ‘Thin crystal of anthracene at room temperature, 


Plate 7, Same crystal as on plate 6 at temperature of liquid oxygen. 


Plate 8. Same crystal as on plate 7 warmed up again to room temperature, 
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Top :—Diffraction pattern obtained from an 
anthracene crystal, orientated so that the two 
molecules in the unit cell are approximately at 
right angles to the incident beam. 


Opposite right is a sketch of the pattern to be 
expected from independent molecules, orientated 
as in the crystal. The contributions of the two 
separate molecules are distinguished by vertical 
and horizontal shading. ‘This sketch is not the 
same as that of the structure factor. 

Points to note are : 
i. The large dimensions (especially width) 


of the bands, 


C 
! 

ii. the occurrence of strong bands at points Bee! = 
where the spot pattern disappears owing to the ile) 
structure factor, = 


iil. the absence of any marked interference 
effects when the contributions from two separate 
molecules overlap. 


Reproduced from Proc. Phys. Soc., 1939, 51, 479 
‘Charlesby, Finch and Wilman). 


Plate A. 
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Plate B. 


Plate C, 
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agreement between the calculated pattern and the diffuse bands observed. This theory 
was put forward in 1939 by Charlesby, Finch and Wilman. It was realized that this 
explanation was only a first approximation, since it did not account for the presence of diffuse 
streaks and spots, which could often be observed superimposed on the diffuse bands. 

It appears to me that both theories are necessary to account for the full pattern. Of 
the total thermal motion of a molecule one component can be considered to form the 
elastic waves continuous throughout at least fairly large areas of the crystal. This component 
leads to streaks and diffuse spots. The remaining component is independent of that of 
other molecules and gives rise to the diffuse bands. But I must emphasize that an inde- 
pendent component must be accepted to account for the diffuse bands. 

Mr. Skrebowski assumes that for a complex crystal such as anthracene the intensity 
of ‘diffuse streaks follows the structure factor of the unit cell. No such calculation has 
been made, nor does the diffuse pattern observed in simpler crystals such as diamond 
follow this intensity distribution. 

In his remarks the author points out that the structure factor for the unit cell is not 
very different from that calculated from each molecule separately, and attempts to sub- 
stantiate this statement with a direct comparison between the two. Unfortunately the 
single orientation he has shown is very badly chosen to make this comparison. In other 
orientations the difference is much greater, and in such orientations the observed diffuse 
pattern does not agree with that calculated from the structure factor, but rather with that 
calculated for each molecule separately. 

Anthracene provides a good example of this distinction. In certain directions the 
amplitudes scattered by the two molecules are equal but opposite in phase so that the resul- 
tant intensity is zero. Spots should therefore vanish (forbidden reflections), as should 
any diffuse streaks in the immediate neighbourhood. Observations of the area where diffuse 
patterns from different molecules overlap show no such phenomenon of interference. 
Strong diffuse bands occur in areas where the interference between the two molecules in 
the unit cell should eliminate them, and does in fact eliminate the Laue spots. On the 
other hand the theory of independent molecular scattering accounts both for the appearance 
and extent of diffuse bands in these areas (plate A). 

In view of the importance of this argument, it is worth considering it in further detail. 
On page 522 of the original paper by Charlesby, Finch and Wilman (Proc. Phys. Soc., 1939, 
51, 479) the difference was pointed out between the two factors in the expressions for 
the spot and diffuse band intensities. 

The expressions for the intensity distribution of the spot and diffuse band pattern 
respectively in a direction h, k, / were given as 

1, =(A?-- B*)x(sin® M,hr/sin® hz) exp (—167? sin? 6kT/ Af), 

I,=NX(A,?+ By?)(1 —exp (— 162? sin?0kT/ Xf). 

p 

A,+tBy is the amplitude of the wave scattered by molecule p in the unit cell. In the 
expression for the spot pattern (d?+B*) is equal to (SA,)?+(SB,)? summed over all 
molecules in the unit cell, the phase difference between different molecules being taken 
into account. In the corresponding expression for the diffuse bands ¥(4,)*-+- 3(B,)? no 
phase difference between molecules is involved. 

In the particular case of anthracene there are two molecules per unit cell p, g. As 
the molecules are symmetrical B,=B,=0. The structure factor for the spot pattern—and 
according to Mr. Skrebowski for the diffuse streaks—is 

(A, +A, or A,*?+A,?+2A,A, cos m(h--k). 

According to the theory of independent molecular vibration the structure factor for the 
diffuse bands is A,?+ A,’. 

The distinction between these two expressions is best seen in diffractions for which 
Ap=Ay (eg. hOl and ORO). The intensities are then proportional to 4A,” cos® h/2 or 
44,° cos’ zk/2 (spots and diffuse streaks) and 24,” (diffuse bands). In the theory of 
independent molecular scattering the intensity of the diffuse bands is modified by the 
molecular structure factor and by a thermal factor while that of Laue spots and diffuse 
streaks in addition depends on the value of the expression cos? 7h/2 or cos? 7k/2. 

If the author’s thesis is correct, the intensity distribution within a diffuse band which 
extends over an area corresponding to several Laue diffraction spots should show charac- 
teristic fluctuations from 44,” to0. Although certain diffuse bands do show some variations 


PROC. PHYS. SOC. LXI, 5 29 


442 F. Skrebowski 


in intensity, there is no doubt that many bands occur with high intensities at points where 
the factor cos* zh/2 or cos? zk/2 would cause them to disappear. Dyffuse bands are often 
observed at points where the spot pattern is eliminated by the structure factor, e.g. (h02) 
diffraction with h odd, or (ORO) diffraction with k odd. In fact the intensity distribution 
in the diffuse bands is proportional to the structure factor of each molecule, and not to the 
structure factor of the unit cell. 

One is therefore driven to the conclusion that the diffuse bands although dependent 
on the atomic arrangement in the molecule are not dependent on the arrangement of the 
molecules in the unit cell. As has long been suggested—and verified by the author—the 
appearance of diffuse bands is a temperature-dependent phenomenon. It seems difficult 
to escape the conclusion that these bands derive from that component of the thermal 
motion of each molecule which is independent of that of its neighbours. ° 

I have mentioned the existence of a number of characteristic differences between 
diffuse bands and diffuse streaks. The following observations summarize these differences : 

(i) In the patterns shown to us the streaks run from one sharp spot parallel to one of 
the spot arrays. The bands on the other hand bear no obvious relationship to the direction 
of the spot pattern. (ii) The intensity of the streaks is modified in the proximity of a 
diffraction spot. ‘The bands appear to be quite independent of spot intensity except in 
so far as they both vanish in directions in which the structure amplitude from each molecule 
is low. (iii) Each streak is relatively narrow. Bands have a width often equal to several 
times the distance between Laue spots. (iv) The intensity of the streaks, in the photographs 
shown, is less than that of spots in the neighbourhood. Diffuse bands at some distance 
from the central spot are far more intense than Laue spots. (v) Although diffuse streaks 
may be observed superimposed on diffuse bands near the central part of the pattern where 
they can be readily distinguished, only diffuse bands are observed further away from the 
centre. (vi) The diffuse streaks are related to the direction of thermal vibrations in the 
crystal, and hence to the crystalline axes. Their intensity is mainly dependent on the 
direction of these thermal vibrations. The diffuse bands are independent of crystal 
orientation except in so far as this determines the molecular orientation. 

In conclusion, I would like to state that the assumption of independent vibration leads 
to a number of deductions which may be readily verified in the diffraction patterns. These 
include (1) the distribution and dimension of diffuse bands in the diffraction pattern; 
(11) the relative intensity distribution within the bands; (iii) the increased intensity of 
the bands relative to the spot pattern as the angle of diffraction is increased; (iv) the 
non-occurrence of interference between bands from different molecules. 

A simple extension of this idea provides an explanation of the diffuse spots seen by 
Preston and Bragg in inorganic crystal diffraction patterns. 

The theory in the simple form in which it was first presented does not explain the 
occurrence of diffuse streaks or spots. 

The alternative theory of continuous thermal vibrations as propounded by Mr. Skrebowski 
explains the appearance of diffuse streaks and spots, though only qualitatively. It does 
not explain the following observations : 

(i) the non-intensification of bands near Laue spots; 

(ii) the non-interference of diffuse bands; 

(111) the occurrence of strong diffuse bands in areas where the structure factor eliminated 
the spot pattern; : 

(iv) the appearance of diffuse bands at large diffraction angles, when the spot pattern 
has vanished. 

The most likely assumption appears to be that the thermal vibration of each molecule 
can be split into two components : 

a component dependent on elastic vibrations throughout the crystal and giving rise 
to diffuse streaks and spots; 1 
a component which is independent of that of its neighbours, leading to diffuse bands. 

From the practical point of view, the latter component is of greater importance in so 
far as it provides a simple and effective method of molecular analysis. 

In view of the increasing interest shown in these extra diffraction phenomena, might 
I suggest that some less confusing terms be applied to diffuse streaks, diffuse spots and 
diffuse bands, now that their causes are being separated. Much confusion in this paper 
arises from the use of a single adjective to describe such diverse phenomena. 
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Dr. H. Witman. Although it seems reasonable, as we have said before, to expect 
the diffuse maxima to become weaker relative to the Laue spot intensity at low temperature, 
comparison of plate 7 with 6 and 8 is largely vitiated because plate 7 appears to correspond 
to a thinner part of the crystal, and the print of plate 7 seems to have had a longer exposure 
than those of plates 6 and 8. The plate itself still contains considerable intensity of 
diffuse maxima. ‘This may be due to considerable heating of the crystal by the rather low 
energy (28 kv.) electrons used. Mr. Skrebowski has thus not conclusively shown the 
dependence on temperature. I have no doubt, however, but that such a temperature 
effect must exist. 

The apparatus is not new in type. No references to previous work with liquid-air 
cooled specimens are given, though essentially identical apparatus, avoiding impurities 
condensing on the specimen, was described by Gen, Zelmanoff and Schalnikoff in 183% 
Hass in 1937, 1938, Was in 1939, and Charlesby in 1945 (Pros. Phys. Soc., 57, 510). 

The author suggests that he has discovered pattern details which had, in fact, already 
been described by Charlesby and Wilman (1942). In plate 5 Mr. Skrebowski has drawn 
three sets of parallel equidistant straight lines through the spots as the Laue zones—this 
would mean that the three crystal axes were all co-planar and normal to the beam, which 
is impossible and inconsistent with the spot intensities even as an approximation. Farther, 
the crystal appears imperfect and the electron beam too broad. 

In the introduction and discussion, Mr. Skrebowski expresses views which contradict 
facts published by Charlesby, Finch and Wilman (1939), and Wilman in 1942 (Proc. Phys. 
Soc., 54, 350), Charlesby (1942) and Charlesby and Wilman (1942), which we now illustrate 
again with slides. ‘The x-ray patterns cannot confirm our electron diffraction results or 
otherwise, owing to the difference in wavelength, energy, crystal thickness etc. K. Lons- 
dale (Proc. Phys. Soc., 1942, 54, 352) has pointed out that the different crystal thickness 
in the two cases will probably lead to a different distribution of the thermal waves. The 
theory of thermal-elastic waves in the crystal can account for the diffuse spots in the x-ray 
patterns and also for the similar spots in the electron diffraction patterns, but has not been 
able to account for the more diffuse bands and areas of intensity hitherto only observed 
in the latter. Both the thermal-elastic wave theory and our molecular scattering theory 
fail to explain the whole of the electron diffraction phenomena, though a mixture of the 
two does appear to account for the main characteristics (plat2 A). 

Mr. Skrebowski confuses groups of the thermal-elastic spots with the broad diffuse 
maxima on which they often, but by no means always, lie and which he has not effectively 
recorded owing to under-exposure and use of plates with unsuitable characteristics. Similarly, 
the X-ray patterns do not show even the thermal-elastic spots if under-exposed. Our 
experimental proof of the correspondence of the diffuse maxima with the patterns which 
would arise from the separate molecules. orientated as in the crystal is a fact which he 
seeks to interpret as due to the thermal-elastic spots being strongest (and the weaker outer 
parts tending to coalesce) where the normal structure factor is high. He admits his F? 
contour diagram, figure 2, is not significantly different from that for F, and Fs, figure 3, and 
that neither show subsidiary maxima within the broad diffuse maxima. In most orientations 
the F? distribution shou/d, however, in certain regions show maxima and minima separated 
by distances of the order of that between the Laue spots. ‘This has not so far been reliably 
observed—we see no appreciable interference of this sort where the maxima of the two 
differently orientated molecules overlap. Moreover, the diffuse bands and areas are 
still relatively strong, as expected from the independent-molecule theory, in parts of the 
pattern distant from the undeflected spot, where the Laue and Bragg spots are no longer 
strong enough to be visible. Hence, it seems to be the F,+F, diagrams, as we used them, 


which apply. 


Professor G. I. Fincu. In discussing this paper J should like in the first place to outline 
the history of the discovery of these thermal effects in electron diffraction patterns. The 
diffuse areas of intensity were first recorded by Finch, Quarrell and Wilman (Trans. 
Faraday Soc., 1935, 31, 1051). In 1937, Finch and Wilman (Ergebn. exakt. Naturwiss., 
1937, 16, 353) were able to say that these areas were associated with the molecular structure 
rather than with the crystal lattice; and in 1938, Finch and Tyson (15th Ann. Rep. British 
Emp. Cancer Campaign, 1938) accounted in a similar manner for the diffuse area patierns 
they had obtained from a variety of carcinogenic compounds. In 1939, Charlesby, Finch 
and Wilman (Proc. Phys. Soc., 1939, 51, 479), made a comprehensive study of the anthracene 
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pattern and gave detailed proof that the diffuse areas originate from the molecules. In 
1942, Wilman (Proc. Phys. Soc., 1942, 54, 350) published a pattern from an inorganic 
substance, cadmium iodide, showing diffuse lines along the main spot rows similar to those 
found with phthalimide (Finch and Wilman 1937) and later, in 1942, Charlesby and Wilman 
(Nature, Lond., 1942, 149, 411) pointed out that these diffuse lines correspond to the diffuse 
spots or lines which had then recently been found in x-ray patterns, and that they occur 
together with these still more diffuse areas of intensity in electron diffraction patterns 
obtained from a wide variety of molecular aliphatic and aromatic compounds, including 
anthracene. 

The first clear recording of the coherent thermal scattering effects with x rays was 
by Preston in July, 1939, some four or five months after the first full publication of our 
results and four years after we first announced their discovery. It is remarkable that 
since then the rdle played by electron diffraction in recording thermal scattering effects 
seems to have been largely overlooked by x-ray workers—I have seen only two references to 
it, by Mrs. Lonsdale and the late Sir William Bragg. 

Thus in 1942, it was clearly shown by us that electron diffraction patterns may exhibit 
three main types of features, (i) normal Laue spots together with, in thicker crystals, more 
or less well-developed Kikuchi lines of normal type, (ii) more diffuse spots due to coherent 
thermal vibrations, like those now so much studied in x rays and, finally, (111) diffuse areas 
of intensity due to molecules scattering as independent units as a result cf incoherent thermal 
motion. Mr. Skrebowski, as far as I can see, has added little or nothing to this picture 
and indeed seems to me to have closed his eyes to the facts which show the occurrence of 
independent molecular scattering. 

I need only draw your attention to the much larger thickness of the crystals used in 
xX rays to show that the absence of these incoherent thermal-motion effects in the x-ray 
pattern is not surprising. The molecules at or close to the surface of a crystal are obviously 
subject to unsymmetrical forces resulting from atoms (or molecules) being on each side 
of them and below, but not above, them in the crystal, so that their motion must be less 
regular than that of the molecules deeper inside a large crystal; and, moreover, there must 
be an absence of the longer standing waves in directions steeply inclined to the main surface 
of a thin crystal sheet. In large crystals there is a small ratio of surface to mass, and the 
relative importance of the diffraction from the surface layers is small, but with the very 
thin crystals necessarily used in electron diffraction it will be much larger, and may even 
predominate in crystals which are only a few molecules thick. This is indeed precisely 
what is found experimentally, when an adequate technique is employed, not only in pro- 
duction of the pattern, but also in its scrutiny. 

I have hinted elsewhere at the potential value of this independent molecular pattern 
to molecular and crystal-structure analysis. The coherent thermal-motion pattern has 
been in this respect, as far as I am aware, of only limited value in x-ray analysis. It will 
therefore, be of some interest if I tell the meeting that Miss Fisher, working in my laboratory, 
has in proof with the Society *, a paper which describes how an appreciation of the inde- 
pendent molecular scattering pattern has enabled her to make a valuable advance in our 
knowledge of the molecular and crystal structure of rubber, a field in which hitherto we 
have had to rely only on x rays as a tool. 


AuTHOR’s reply. Some misundeistanding has arisen about part of my paper which 
I should like to correct. Contour lines of the intensity scattered from a unit cell (Agure 2) 
are drawn for the purpose of indicating where the cross-grating spots (if they appear) 
are strongest. ‘Thus a cross-grating spot lying between, say, lines 140 and 180 will have 
an intensity between these values. These contour lines tell nothing about the diffuse 
background between cross-grating spots. ‘Their sole purpose is to show where strong 
cross-grating spots will tend to group, as well as the shape of these groupings. Thus 
Dr. Charlesby’s discussion on this part of my paper is inapplicable. 

Owing to the length of the discussion I should like to recapitulate briefly the issues 
involved. Charlesby, Finch and Wilman claim the following : 

1. ‘The diffuse areas are primary phenomenon. 

2. They are without any detailed structure. 


8.) Whey, can be explained by the molecular scattering theory which predicts their 
positions and shapes. 


* Now published (Phys. Soc. Proc., 1948, 60, 99). 
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tho Sometimes diffuse streaks and spots appear, but these have nothing to do with the 
diffuse areas. These streaks and spots can be explained in terms of the Faxen- 
Waller theory and are similar to those found in x rays recently. 


In my paper I attempted to show : 


1. The diffuse areas (in the Charlesby, Finch and Wilman sense) are secondary phe- 
nomena. 

2. They always have a definite fine structure which is more prominent in thick crystals. 

. With thin crystals this fine structure is very often smeared out and the appearance 

of diffuse bands or areas without any fine structure prevails. 

4. The diffuse streaks and spots are grouped into bands or areas and the position and 
shape of these groups is equally well (if not better) predicted by structure factor 
considerations than by the molecular scattering theory. 


W 


The Charlesby, Finch and Wilman view has many conflicting elements. If the diffuse 
streaks and spots exist independently, then why do their contour lines drawn for the diffuse 
areas include these streaks (figures 10 a, 14a) ? 

In the discussion there was no mention of the different appearance in the pattern for 
thick and thin crystals. This is, however, an important point in the present controversy. 
‘The pattern in plate 2 might be taken for the diffuse areas (in the Charlesby, Finch and 
‘Wilman sense) but for plate 3 which shows that these areas havea fine structure. If, however, 
the crystal shown on plate 2 were thinned a large number of cross-grating spots would appear. 
This would tend to obliterate the fine structure and produce something like that observed 
by Charlesby, Finch and Wilman (1939), i.e. diffuse areas. Jn the discussion none of these 
authors attempted to explain the absence of diffuse areas on plates 2 and 4. 

One of the main claims of the molecular scattering theory is that it can predict the 
position of the diffuse areas satisfactorily. In my paper I have shown that these can be 
found equally well from the structure factor consideration. They coincide with the regions 
where the structure factor values are high. ‘This is not a coincidence as Dr. Charlesby 
suggests, for I have calculated this with widely different crystal orientations. Unless he 
produces quantitative contour lines for a particular photographic plate his argument will 
have little value. In fact, the qualitative treatment adopted by Charlesby, Finch and 
Wilman (1939) is open to criticism. Often their diffuse areas are very sma!l (figures 10a, 
14 a) and cover 1—2 reciprocal lattice cells. A quantitative treatment would give a different 
result. The diffuse areas in the latter case are much broader (figure 3) and are without 
sharply defined limits, in contrast to the Charlesby, Finch and Wilman estimations. 

Drs. Charlesby and Wilman in interpreting the diffuse pattern apply either the Faxen- 
Waller theory or the molecular scattering theory or a mixture of both without proving 
first that these can be applied simultaneously and are not mutually exclusive. 

Experiments have been performed on other molecules in order to test the molecular 
scattering theory further, for this theory requires large and highly symmetrical molecules 
to produce the effect (Charlesby 1942). Now phenanthrene is not a symmetrical molecule 
and its unit cell contains four molecules. If, as the molecular scattering theory requires, 
each of these molecules were scattering independently, one would expect the resulting 
pattern to be a maze of diffuse bands and areas. Experimentally, however, the pattern 
is very simple and highly symmetrical (plate B). It appears to follow the structure factor 
and is composed of diffuse streaks (visible on the original plate). Naphthalene was the 
other compound examined. Its crystalline structure is almost a replica of anthracene 
somewhat reduced. According to the molecular scattering theory, the naphthalene 
molecule, being much shorter than anthracene, should give very broad diffuse areas. 
The pattern observed (plate C) has diffuse areas (with fine structure) of the same angular 
spread as anthracene. 

I found Dr. Charlesby’s argument about the use of photographic plates in electron 
diffraction unconvincing. 

In replying to Dr. Wilman, I would point out that he has misunderstood the indexing 
shown on my plate 5. The lines shown are not Laue zones. His comment regarding my 
apparatus is ill-founded in view of the statement made by his collaborator, Dr. Charlesby 
(Proc. Phys. Soc., 1942, 54, 383, § 2), ‘‘ Experimental difficulties have so far prevented 
direct investigation of the effect of temperature on the diffuse band pattern ”’. 
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ABSTRACT. Delayed fracture experiments were made in tension, torsion and com- 
pression, the times to fracture under a series of constant loads being determined. The 
materials tested were brass in ammonia vapour, polymerized methyl methacrylate (p.m.m.) 
in liquid carbon tetrachloride, and glass in air. 

In all cases delayed fracture is attributed to the gradual spread of Sigeks caused by the 
preferential attack by the surrounding medium of the highly stressed material at the ends 
of the cracks. In tension, fractures were transverse and in torsion they were helical. In 
the case of brass and glass, subjected to tension and torsion, fracture of the respective 
materials occurred in not greatly differing times when the materials were subjected to the 
same maximum nominal principal tensile stresses. P-m.m. withstood higher nominal 
principal stress in torsion than in tension for equal times to fracture. "The p.m.m. results 
were complicated by the absorption of carbon tetrachloride, which was greatly accelerated 
by tension or torsion stresses. During the time occupied by the tests, absorption was 
inappreciable when the p.m.m. was stress-free or stressed in compression. 

In no case did delayed fracture occur in compression. 
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ANY materials if left under constant load will ultimately fracture. ‘The 

delayed fracture is usually associated with attack of the substance under 

stress by the surrounding atmosphere. In the absence of stress, the 
attack may be inappreciable, but when the substance is stressed, very localized 
attack may take place, leading te the formation of cracks and ultimately to fracture. 
For metals, the phenomenon is usually called stress corrosion cracking. Under 
certain conditions of heat treatment, many aluminium- and magnesium-rich alloys 
suffer from stress corrosion cracking when exposed under stress to atmospheres 
containing chlorides; copper alloys are attacked by ammonia, and iron alloys by 
many salt solutions (Mears et al. 1944). The very remarkable delayed fracture of 
glass in air, for which a three-fold reduction in stress leads to a proportional increase 
in time to fracture of 10°: 1, has been shown by Preston and his collaborators 
(1946) to be mainly due to the attack of the atmosphere. Organic plastics also give 
rise to delayed fracture, especially if surrounded by organic liquids or vapours. 
A good example is delayed fracture of polymerized methyl methacrylate (p.m.m.) 
when immersed in carbon tetrachloride, which gives inappreciable attack during 
the same period when the p.m.m. is stress free. 

In nearly all the delayed fracture experiments so far made the materials have 
been tested in either tension or bending, whereas in practice materials have to 
withstand other stress systems. In the present experiments delayed fracture of 
materials subjected to tension, torsion and compression, has been investigated. 


§2. MATERIAL 


In the present series of experiments three materials were tested. The first 
was brass to Spec. B.S.251 of the following composition: Cu 62-4, Sn 1-2, Zn 36-4. 
Before machining, the material was heat-treated in air at 500°c. for 45 minutes 
and then cooled in air. After machining, the test pieces were heat-treated in air 
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at 270°c. for one hour and cooled in the furnace. . Control tests were made in 
tension, torsion and compression. The nominal ultimate tensile stress when 
tested in air was 24-9 ton/in? and the 0-1°% proof stresses were 9-5 tons/in? (tension), 
and 9-6tons/in? (compression). The figures quoted are the means of four tests. 
The delayed fracture experiments were done in an atmosphere saturated with the 
vapour of aqueous ammonia solution of 35% w/w. NHs. 

The second material tested was polymerized methyl methacrylate (p.m.m.) 
sheet *. After machining, all test pieces were heated in air at 80° c. for half an hour 
and cooled in the oven. The ultimate tensile strength at a rate of increase of 
stress of about 5000lb/in?/min. was about 10,000lb/in?. Delayed fracture 
experiments were made with the test pieces immersed in carbon tetrachloride. 

The third material was a soda-lime glass in the form of tin. diameter rods, 
the chief ingredients of which were present in the following percentages: silica 69, 
alumina 3, lime 6, magnesia 3,soda 18. It was annealed by the manufacturer, and 
subsequent measurements of internal stress by photoelastic methods indicated 
a mean value of surface compressive stress of about 6001b/in® with a standard 
deviation of about 400lb/in?. The rods were packed without wrapping and 
tested without previous surface treatment. Delayed fracture experiments were 
made in air. 


§3. APPARATUS 


For delayed fracture experiments in tension and compression a spring-loaded 
cage was used. The test piece was screwed into plates at each end. Three 
screwed rods fitted with nuts and springs passed through each end'plate. For 
the tension tests the springs were between the plates, one end of each 
bearing on an end plate and the other bearing on nuts screwed on to the rods. By 
tightening the nuts, the tension on the test piece was increased. For compression 
tests the springs were outside the end plates. In both tension and compression 
tests three Huggenberger tensometers were mounted on each test piece at angles 
of 120° to each other, and the nuts which loaded the springs were adjusted until 
the readings on each tensometer corresponded to the strain corresponding to the 
required stress. This arrangement ensured axial loading. For the torsion tests 
pulleys were screwed on to each end of the test pieces and located by lock nuts. 
Each pulley carried two separate cables attached to its rim at each end of a diameter. 
The cables were in contact with the pulley rims for an arc of about 90°, and then 
they left the pulley in opposite parallel directions. The directions of the cables on 
the pulleys at each end of a test piece were parallel and were arranged so that when 
all four cables were under equal tension the test pieces was subjected to equili- 
brating end torques. The equality of the tensions was arranged by attaching the 
two upward cables to the opposite ends of a beam and by attaching the two 
downward cables to the opposite ends of another beam ; the whole apparatus 
could be put inside the cages used for the tension and compression experiments, 
which were adapted to apply central loading to the two beams. For small loads 
the beams were pulled apart by dead loading. 

The brass test pieces were 2}in. long threaded din. B.S.F. for a length of 
3 in. ateach end, and reduced to 5/16 in. diameter over a central length of 3in. 
They had a central axial hole of 3/16in. diameter running the full length of the test 
piece. ‘The p.m.m. test pieces were 3 in. long threaded } in. B.S.F. for a length 


* Containing 8% dibutyl phthalate as plasticizer. 
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of lin. at the ends and reduced to Lin. diameter over a central length of 3 in. The 
glass test pieces were round parallel rods of }in. diameter. ‘They were tested in 
bending instead of in tension, on account of the difficulty of obtaining axial 
loading. A four-point loading rig was used, the length subjected to uniform 
bending moment being 3in. For torsion, the glass rods were cemented into 3 in. 
long steel sleeves with an organic cement. Pulleys were screwed on to the sleeves, 
and the tests were made in the same rig as that used for brass and p.m.m. 


§4, EXPERIMENTS 


The test pieces were set up in the apparatus already described and the times to 
fracture under a series of constant loads were recorded. ‘The mean temperature 
during the tests was 22°c. This was thermostatically maintained to within + 1°c. 
for the p.m.m. tests, but for the brass and glass tests the temperature was not so 
closely controlled. Figure 1 shows the results for brass, figure 2 those for p.m.m. 
and figure 3 those for glass. For the torsion tests the nominal maximum principal 
stress computed on the basis of a linear stress distribution is plotted against the 
logarithm of the time to fracture. For the tension tests the nominal average 
tensile stress has been plotted. There were five repetitions of each brass experi- 
ment, three repetitions of each p.m.m. experiment and twelve repetitions of each 
glass experiment. ‘The points shown in figures 1, 2 and 3 are median points 
(half the test pieces failed at a lower stress). "The median was chosen rather than 
the mean, as the former can be computed without waiting for all test pieces to 
break. ‘The standard error of the median of a normal population is about 1-25 
times that of the mean, so that the median is not a greatly worse statistic than 
the mean. ‘The standard errors of the median points (on the log (time) scale) 
varied somewhat with stress, the standard error increasing with decreasing stress, 
and becoming very large at the lowest stresses where the curves tend to become 
parallel with the time axis. Except at the lowest stresses, the standard errors did 
not vary greatly with stress, and average standard errors of the median points have 
been computed from the mean of the ranges of results for each stress, omitting 
the results for the lowest stresses. The values for tension and torsion respectively 
were 0-09 and 0-08 for brass, 0-10 and 0-05 for p.m.m., and 0-5 and 0-5 for glass. 

When brass in ammonia and p.m.m. in CCl, were stressed in compression at 
the highest stresses which could be applied without buckling the test pieces, no 
delayed fractures occurred and no cracks were noticed. Even when a short 
cylindrical piece of p.m.m. was immersed in CCl, and compressed so that it was 
deformed into a barrel shape, no cracking was observed. Compression tests on 
glass were not made, but it is known that glassis many times stronger in com- 
pression than in tension, so that delayed fracture in compression would not occur 
at stresses at which it occurred in tension. 

When stress-free brass was left in the ammonia vapour it became covered with 
a blue deposit of more or less uniform thickness. When stress-free p.m.m. was 
left in CCl, no attack was detectable during the time occupied by the delayed 
fracture experiments, but after some months immersion at 22°¢. the surface of 
the p.m.m. was seen to be swollen and covered with a gelatinous layer and the test 
piece had lost weight, showing that solution was taking place. The soda-lime glass 
was not noticeably attacked by the atmosphere when stressed or unstressed, but itis 


known that attack takes place after prolonged exposure. At elevated temperatures 
glass is rapidly attacked by steam. 


JIN? 


TENSILE STRESS LB. 


Delayed fracture under tension, torsion and compression 449 


20,000 


ene | FIG.) BRASS IN AMMONIA — 


10,000 


5000 


ea 


- a FIG. 2 P.M.M. IN CARBON 
| TETRACHLORIDE 


® 


TENSION 


FIG.3 GLASS IN AIR 


©—» INDICATES THAT MEDIAN 
1S FURTHER TO RIGHT 


TORSION 
BENDING 


10! 102 10> 104 105 106 107 
TIME IN SECS. - LOG SCALE 


450 C. Gurney and Z. Borysowski 


§5. APPEARANCE OF FRACTURES 


In the tension tests fracture of the brass test pieces started by the formation of | 
approximately transverse intercrystalline cracks. Many test pieces had more than 
one crack of about the same depth, and final fracture often consisted in bending 
off the material between two cracks that had started from opposite sides of the test 
piece. There was no appreciable reduction in area. In the torsion tests, fracture 
of the brass test pieces started by the formation of cracks at about 45° to the axis 
at the test piece. Completion of the fracture usually involved bending and tearing 
of the test piece. 

The tension fractures in the p.m.m. test pieces were transverse and the torsion 
fractures were helical. Examination of the fractured test pieces showed a very | 
interesting phenomenon. ‘Test pieces which had been subjected to comparatively | 
low stresses showed fractures which appeared to start at the surface, as might be 
expected; but fracture in the test pieces which had been subjected to higher | 
stresses started, not from the surface, but from a point in the interior. Apart from 
examination of the actual fractures, there was other unmistakable evidence that | 
fracture started internally; after tests had been in progress for some time, both | 
tension and torsion test pieces developed clearly visible internal cracks which were | 
transverse in the case of tension and helical in the case of torsion. 

The fractures of the glass rods call for no special comment as their appearance |f 
is the same as that of fractures occurring during quick-loading tests. The bending | 
fractures were usually S-shaped and the torsion tractures were helical. 


§6. DISCUSSION 


In all three materials delayed fracture is regarded as being due to the gradual 
spread of cracks. The cracks were clearly visible in the brass and p.m.m. test 
pieces. ‘They were not visible in the glass test pieces, but there are good reasons |ff 
for attributing the delayed fracture of glass to the growth of very small cracks |] 
(Gurney 1947 a). When material containing cracks is stressed there is a concen- |f 
tration of strain energy at the ends of the cracks and consequently the material | 
at the end of the cracks has a different chemical potential from that of the rest of | 
the material; and under some stress systems it is preferentially attacked by any | 
chemically active surrounding medium. Appropriate thermodynamic relations | 
when the chemical attack takes the form of swelling and solution are given in a | 
recent paper (Gurney 1947b). Other cases could be calculated using similar | 
methods. ‘The thermodynamic relations are somewhat complex and there is no | 
reason to suppose that any simple parameter of the stress system can be used to | 
express results for different stress systems and different chemical substances. 
Nevertheless the fact that delayed fracture did not occur in compression suggests | 
that the maximum principal tensile stress may prove useful in expressing results | 
for simple stress systems such as unidirectional stress and torsion. Accordingly 
this parameter has been plotted as ordinate in figures 1, 2 and 3. 

For brass and glass the curves for tension and torsion are not very different, 
and in the absence of better information the maximum principal tensile stress 
criterion might be used as a rough working rule. There are several factors which 
would tend to cause the tension and torsion curves for brass to differ. The 
assumption of linear stress distribution in torsion would be in error at the higher 
stresses. At the lower stresses cracks occupying a considerable part of the cross- 
section developed before fracture occurred. The nature of the stress distribution 
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would be entirely altered by the presence of these large cracks. In particular, the 
stress at the ends of the cracks, which presumably is an important factor affecting 
the rate of crack growth, might be very different in the torsion tests and the tension 
tests, even though the nominal maximum principal tensile stresses assuming 
uncracked material were equal. In the case of glass these factors are inoperative 
as the assumption of linear stress distribution in torsion is probably accurately 
fulfilled, and fracture occurs before cracks of appreciable length compared with 
the diameter of the test pieces have developed. The experimental points for glass 
in torsion and bending are not significantly different, and a single curve has been 
drawn to represent the results. 

The difference between the tension and torsion curves for p.m.m. is greater 
than that between the corresponding curves for brass or glass. The factors which. 
have been suggested to account for the difference in the case of brass are also 
operative for p.m.m., but it seems doubtful that these effects would be so much. 
greater for the p.m.m. test pieces than for the brass test pieces as to account for the 
large difference. Examinations of the fractured p.m.m. test pieces showed the 
peculiarity already mentioned—that at high stresses fracture started internally. 
This is evidence for an additional factor affecting the p.m.m. results which finds. 
aready explanation. When stress-free p.m.m. is immersed in CCl, for prolonged 
periods, it is found that the surface is swollen by absorption of CCl,, and that some 
of the p.m.m. has actually gone into solution in the liquid. The first effect of 
immersing stress-free p.m.m. in CCl, therefore will be a gradual absorption of the 
CCl, accompanied by swelling, but after prolonged immersion solution in the 
CCl, takes place. The effect of stress of the effects can be deduced from general 
thermodynamic reasoning. Natural processes tend to a reduction of the free energy 
of systems, and part of the free energy of a system under load lies in the poten- 
tial energy of the external forces. If arod under tension absorbs fluid, the external. 
forces do positive work, and therefore reduce their potential energy for two reasons. 
Firstly, there is a direct increase in the length of the rod due to swelling, and 
secondly, there is an increase in the length of the rod due to the reduction in 
Young’s modulus consequent upon swelling. The effect of tension, therefore, is 
to increase the rate of swelling. The effect of torsion is also to increase the rate of 
swelling; the reduction in shear modulus of the swollen test pieces allows the end 
couples to rotate and do work. 

The swelling of the surface material of a rod under tension or torsion causes two 
effects. The swollen material is weakened and therefore more likely to fracture; 
but its size is increased and its elastic moduli (Z and G) are reduced, while the 
elastic constants and dimensions of the unswollen interior material remain the same. 
The swelling of surface material tends to set up a system of internal stress in which 
the surface material is in two-dimensional compression. Both in tension and 
torsion, therefore, the maximum principal tensile stress in the surface is reduced. 
A further reduction is caused by the change in elastic moduli. Absorption of 
liquid by surface layers therefore weakens them, but at the same time causes a 
reduction in the stress they have to withstand. 

These ideas lead to a plausible explanation of the experimental results that at 
low stresses fracture of the p.m.m. test pieces started from the outside, and at high 
stresses fracture started from the inside. It is only necessary to assume that at 
low stresses the weakening of surface material by absorption of CCl, exceeds the 
reduction in the stress associated with the absorption; fracture would then start 
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from the outside. On the other hand, at high stresses the reduction in stress due 
to swelling and change in elastic moduli exceeded the weakening of the material 
due to absorption, so that fracture started in the interior, where the stress was 
already increased by the internal tension stress caused by the swelling of the surface 
layers. 

The broad conclusions to be drawn from the experiments are that delayed 
fracture does not occur in compression. For brass in ammonia and glass in air, 
delayed fracture in tension and torsion occurs in approximately equal times at 
equal principal tensile stresses. P.m.m. in carbon tetrachloride will withstand 
higher principal tensile stresses in torsion than in tension for equal times to fracture. 


ACKNOWLEDGMENT 


The delayed fracture experiments on glass were made by Mr. S. Pearson. 
Measurement of internal stress was done at the Research Laboratories of the 
‘General Electric Company. 


REFERENCES 
Gurney, ©., 1047 a, Proc. Phys. Soc., 59, 169 ; 1947 b, Ibid., 59, 620. 
Mears, R. B., Brown, R. H., and Dix, E. H., 1944, Symposium on stress corrosion cracking 


of metals (Philadelphia : A.S.T.M.). 
PRESTON, F. W., Baxer, T. C., and GLatsHart, J. L., 1946, J. Appl. Phys., 17, 162. 


8-Disintegration and the Sargent Diagram 


By INO FEATTHER ann“ Ho Os WeRICHARD SON 
Department of Natural Philosophy, University of Edinburgh 


IMS. received 6 May 1948; to be read 5 November 1948 


ABSTRACT. Sargent diagrams constructed for negative-electron-active species having 
80< Z< 94 (as also for the restricted group with 80< Z< 84), and for both negative-electron- 
active and positron-active species having Z< 22, are discussed in the light of theory. The 
majority of first-forbidden transitions in the group of heaviest radioelements give a Sargent 
line parallel to the allowed line, corresponding to disintegration constants smaller by a 
factor of roughly 10 only than the disintegration constants of allowed transitions of equal 
energy. ‘Transitions with no parity change between states of zero spin appear to give a 
second-forbidden line also parallel to the allowed line. 

With the lightest radioelements ‘“‘ super-allowed”’ transitions are identified which 
require Wigner selection rules for their interpretation, but there is also a group of transitions 
which appears to be ordinarily allowed (i.e. subject to spin and parity selection rules only). 
‘Compared with these allowed transitions in the lightest elements, allowed transitions in 
the heaviest radioelements are less probable by a factor of about three only, when comparison 
is made for equal energy release and account is taken of the effect of the coulomb field of 
the nucleus. 


On the assumption that all ((y¢i}) nuclei have even parity and zero spin in the ground 


state it is concluded that the f-particles of maximum energy observed with RaB, RaC 
and MsTh, are emitted in transitions to excited states of the final nuclei RaC, RaC’ and 
MsThg, respectively. This conclusion makes it likely that the accepted value of the total 
y-ray emission from RaB (Gray 1937) is too small, that the total disintegration energy 
Ra.C.C’ is 3:78 Mey. and that MsTh, emits a low-energy y-radiation which has not as yet 
been detected. On the same assumption it is concluded that the disintegrations USS 
and UXg. Uy possess features at present unknown or insufficiently understood. 
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REVISED Sargent diagram for the heavy radioelements has recently 
A been published by one of us (Feather 1948) showing strong evidence 

for a hitherto unidentified “line ”, almost parallel to the line for allowed 
transitions, and not very distant from that line. If the suggestion that the new 
line can belong only to first-forbidden transitions be accepted, then obviously 
the degree of forbiddenness corresponding to the other lines, previously identified, 
and lying lower on the diagram, must be revised—the first-forbidden line of the 
old diagram becoming the second-forbidden line of the new diagram, and so on. 
It is one object of this note to point out that this assignment has the support of 
theory: the occurrence of first-forbidden transitions, differing relatively little 
in probability from allowed transitions of the same energy, has been predicted 
by Konopinski and Uhlenbeck (1941), at least for the heavy elements. For the 
lightest radioelements, according to these authors, such transitions should not 
occur : the calculations show that a much larger difference in transition proba- 
bility is to be expected as between allowed and first-forbidden disintegrations of 
equal energy. 

The formal predictions of the theory are briefly as follows. Let R), be the 
ratio of the disintegration constants, A,(W), A,(W), for allowed and first-forbidden 
transitions in which the total disintegration energy (including the rest-energy 
of the electron) is Wmc?. Let Z be the charge number of the f-active nucleus, 
and let the nuclear radius be ph/2z7mec. Then when 


CL LO ee | We ihn ae (1) 
: being the fine-structure constant, 1/137, we have, alternatively, 

IGA PA ae NG ahaa (2) 
or iit UA’ a? a nn Oe e (3) 


depending upon the change of spin in the transition. If /is the total spin quantum 
number of the nucleus, (2) applies to those transitions which involve a change 
of parity of the nuclear state and for which AJ=0 or +1 and (3) applies when 
Al= +2. (Here it might be noted that first-forbidden transitions with AJ = +2 
are possible only on the basis of the axial vector and tensor forms of interaction, 
which are governed by Gamow-Teller selection rules.) The theory also predicts 
that the shape of the energy spectrum should approximate to the allowed shape 
when (2) applies, and differ considerably from this shape, showing a deficit of 
medium-energy f-particles and an excess of high-energy particles, when (3) 
is valid. 

For the heavy elements (80 <Z <94) for which (1) is amply fulfilled, (2) points 
to a first-forbidden line parallel to the line for allowed transitions and distant, 
approximately, by AlogA=1-0 (8-5 <(«Z/2)-? <11-7)—and (3) to a lower-lying 
line, quite steeply inclined to the allowed line at its upper end, but still representing: 
transitions which are only once forbidden. On this analysis of the situation there 
is clearly a good case for our identifying the newly-recognized line with the first- 
forbidden line predicted by (2) : it exhibits the close parallelism with the experi- 
mental line for allowed transitions which is required, and the absolute value of the 
mean separation (0-75 <A logA<1-1) agrees remarkably well with prediction. 

Having elicited this measure of agreement, the next object of this note is to: 
continue the comparison of experiment and theory in respect of transitions of 
higher degrees of forbiddenness. 
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The theoretical result concerning Rog, the ratio of the disintegration constants 

for allowed and second-forbidden transitions of energy Wme?, is that 

Roo~(aZp W/2)-? Mio cholo d (4) 
Or Ros~(p W)* SNOMED 6210 (5) 
depending upon whether, without change of parity, AJ= +2 or See (Again, 
the higher value of the spin-change, as in the case of first-forbidden transitions, 
is permissible only under Gamow-Teller selection rules.) — 

In order to make further comparison with least ambiguity due to variations 
in Z and p, we restrict consideration at this stage to disintegrations of B-active 
nuclei having Z=81, 82 and 83 (204<A<214). Figure 1 contains the 27 
representative points for these disintegrations, now considered separately from 
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the 23 points for 87<Z<93 with which they were combined in the Sargent » |] 
‘diagram recently published. As starting point we take the point for RaE. Here 
Z=83, A=210, W=3-29 and Ry, =616 (AlogA=2-79). Writing (4) for this 
-disintegration in the form 


Rog =Sa(aZ/2)*. (pW) 


on the assumption that a second-forbidden transition, [=2—+I]=0, is involved 
(Greuling 1942), with (ph/2mmc)=1-5 x 10-4? (cm.), giving p-1=43-2, 
we obtain a=0-066.. The factor 5 (the value of 2/AJ/+1, for RaE.F) has been 
introduced into the equation for Ry, because the transition is taken as one in which 
‘the total spin decreases (to zero) rather than increases (Marshak 1942, Greuling 
1942) ; the constant a then refers to the most probable type of second-forbidden 
transition (having AJ= +2). We do not comment further on the absolute 
value of a, which appears as considerably smaller than “ of the order of magnitude — | 
unity ’’—this aspect of the problem has already been treated in detail by others 
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(Marshak 1942, Greuling 1942), who find in it evidence for narrowing the choice 
of possible interaction forms to the polar vector (Fermi rules) or tensor (Gamow- 
Teller rules) forms—but we attempt to follow the assumption that a closely 
similar value of a is appropriate to the other second-forbidden transitions to 
which (4) applies (within our chosen range of Z). Curves Ila and Ila’, therefore, 
have been constructed, the latter according to (6) through the representative 
point for Rak, and the former parallel to it and higher in the diagram by 
A log A=0-699 (=log 5). In constructing these curves, which are of course 
drawn with respect to the allowed curve O as reference curve, the small variations 
in Z and p, possible for the various points in the diagram, have been neglected. 
Curve O, it should be pointed out, is appropriate only when AJ=0 or +1; 
curve O’ belongs to those allowed transitions for which AJ = —1 (Alog\ =0-477). 
Similarly for the curves [a and Ia’ which belong to the first-forbidden transitions 
governed by (2) : Ia is applicable when AJ=0, +1 and Ia’ when AJ= —1 *. 

In considering the whole diagram at this stage it will be seen that one point 
only lies between curves Ia’ and Ila, and that not more than two “ observed ” 
points lie below curve IIa’. Formally, therefore, at least 19 of the 22 ‘‘ observed ”’ 
points would appear to be satisfactorily accounted for—even though four of the 
five ‘‘ unobserved ”’ points, which lie below curve Ia’, require further considera- 
tion. We shall return later to the general question of the ‘“‘ unexplained ”’ points, 
for the present we direct attention to the group of four transitions which has 
emerged as second-forbidden according to (4), through our construction of curves 
IIa and Ila’ as above described. In addition to the ground-to-ground state 
transition RaE.F, we have the ‘“‘ unobserved ”’ ground-to-ground state transition 
RaD.E (for which a disintegration energy of 60+4kev. and an intensity of 
<0-01 has been assumed) and two transitions ThC.C’ (those leading to the 
ground state of ThC’ and to the excited state of 726 kev. energy, respectively). 

In the first place it is extremely satisfactory that the two ground-to-ground 

. State transitions belonging to the successive disintegrations RaD—-RaE—~RaF 
should be found grouped together in this way. The initial and final nuclei, 
2i9RaD and *{{RaF, are both (ve?) nuclei, and in the ground state each would 
be expected to have J=0 and even parity. Thus the f-disintegrations RaD.E, 
RaE.F which involve the same intermediate state (the ground state of RaE) 
must be expected to be characterized by equal and opposite values of AJ and by 
the same parity-change (either each “‘ yes”’ or each “no”’). The two transitions 
should, in short, be forbidden to identical degree. This is exactly what is found. 
(To be precise, the point for the ground-to-ground state transition RaD.E ought 
to lie on curve IIa, if that for RaE.F is on curve IIa’, but it cannot be held that the 
distinction is important having regard to the uncertainty as regards intensity 
and the extent of the extrapolation involved in the comparison just made.) 

We next discuss the two points for ThC.C’. In the successive disintegrations 
ThB+ThC+ThC’ the situation in respect of A and Z is similar in all respects 
to that just considered: 7J3ThB and *s7ThC’ are each of them (¥e") nuclei ft, 


* According to Marshak (1942) curves O’ and Ia’ belong strictly to those transitions characterized 
by I=1>/=0. Other transitions for which 4J=—1 should have representative points between 
buirves O and O%, if allowed, or, if first-forbidden, between Ja and Ia’. See ‘etter by Strachan in 
this]issue, p. 481. 

+ The classification ‘Gkeah Cos etc., of this paper refers throughout to neutron and proton 


‘ 4 ies ; a ad ay . 
numbers, i.e. is a classification in respect of ( Z y not of (7). 
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thus we should expect the ground-to-ground state transitions ThB.C, ThC.C’ 
to be identically forbidden. We have mentioned that the point for ThC.C’ 
falls within the limits of curves Ila and Ila’, but in this case the corresponding 
point for the previous disintegration does not. In fact the point for the ground- 

to-ground state transition ThB.C is the single point which on figure 1 lies mid-way 
between curves Ia’ and IIa. Formally this point might be interpreted as being 
on acurve Ib, determined by (3)—then |A/| would have the same value, 2, for the 

transitions which are thus associated—but we should not in this way solve the 

problem, for the parity changes would then be different. If our assumptions 
regarding spin and parity of the ground states of ({¥gi) nuclei are correct, there 

would, in fact, appear to be only one solution—the ground-to-ground state 

transitions ThB.C, ThC.C’ must be of a special type which we have not hitherto 

considered. This suggestion might be thought to be strengthened by the 

empirical fact that the two representative points concerned are very closely the 

same distance below the line for allowed transitions—and, from the point of view 

of theory, by the fact that transitions /=0— J=0 constitute a naturally special 

class when certain of the five standard interaction forms are employed (this is 

untrue only with the scalar and polar vector interactions, which carry the Fermi 

tules). We adopt this suggestion, then, drawing the line OO, parallel to O and 

distant by AlogA=1-85—and we make the further natural assumption that the 

new line (OO,) is characteristic of a whole class of /=0—0 transitions generally 

(within the chosen range of Z). ‘That this class of 00 transitions is that for 

which there is no change of parity may be deduced from the following con- 

siderations. ‘The representative point for the main ThB.C transition falls on 

the allowed line. ‘Thus there is no parity change in this 8-transition, which leads 

to the 238 kev. state of ThC. ‘The y-ray emitted in transition from this state 

to the ground state is converted internally as a magnetic dipole radiation. The 

excited and ground states, therefore, have like parity. The ground-to-ground 

state 6-transition ThB.C is likewise, therefore, a transition without change of - 
parity. 

Now 0-0 transitions without change of parity are permitted according to 
the scalar and polar vector forms of interaction as ordinary allowed transitions, 
and such transitions are completely forbidden according to the axial vector and 
pseudoscalar forms of interaction. On the tensor form only are they permitted 
merely as (second) forbidden transitions. It appears therefore that our line 
OO, can only be explained, on the basis of a single (unmixed) interaction, as a 
‘“second-forbidden ”’ line for 0-0 transitions obeying Gamow-Teller selection 
rules and describable in detail using the tensor form of interaction. It is satis- 
factory that theory in this case predicts a line which is parallel to the allowed 
line (Konopinski 1943). It is also satisfactory that the distance of OO, from the 
allowed line is approximately twice the distance of Ia from this line, for amongst 
the transitions to which Ia might belong are the 0-0 transitions with change of 
parity which, on the basis of the tensor form of interaction, behave as first- 
forbidden. Using an obvious notation we might in fact say that a line OO, 
coincides roughly with the line labelled Ia in figure 1. 

We have obviously digressed from our discussion of the two points for ThC.C’ ; 
returning now to this specific topic we may note that whilst one of them has been 
taken out of the group to which the curves IIa and Ila’ apply, there is nothing to 
contradict the identification of the other as the representative point of a second- 
forbidden f£-disintegration having |A/| =2. 
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A third case of successive f-disintegrations for which the initial nucleus is 
an (even) nucleus is provided by the sequence RaB-RaC-RaC’. Although 
it cannot be claimed that sufficient is known regarding these much-studied 
disintegrations to provide full information for a discussion, if we make the same 
assumptions as before regarding spin and parity certain conclusions follow without 
question. The most important of these is that the ground-to-ground state tran- 
sition RaB.C is “ unobserved”. This is in line with the opinion originally held 
(cf. Ellis 1934), but it has been rendered doubtful by results of Gray (1937). 
The present argument is as follows. We may assume for the moment that it is 
not certain whether the ground-to-ground state transition RaC.C’ is that of 
3-173 Mey. maximum energy given by the experimentally determined end-point 
of the RaC f-particle spectrum, or an unobserved transition of 0-606 Mev. higher 
energy (Ellis 1934, Bothe and Maier-Leibnitz 1937), but in either case it appears 
that the transition is at least second-forbidden. On the other hand, the experi- 
mentally determined end-point of the RaB f-particle spectrum at 0-72 mev. 
energy (Constantinov and Latyschev 1941) belongs to a transition which is 
obviously allowed. This is enough to show that the two transitions are not 
consecutive ; they cannot have the ground state of RaC as a common intermediate 
state; in other words the ground-to-ground state transition RaB.C is “unobserved”. 

Concerning the disintegration RaC.C’ again, it has just been stated that it 
is not certain, on the direct experimental evidence, whether or not the ground-to- 
ground state transition is observed. We can now use figure 1 to establish the 
conclusion that it is almost certainly not observed. Evidence from the long- 
range s-particles of RaC’.D shows that a state of RaC’ of 2:88 mey. energy is 
excited in the previous f-disintegration, and it is practically certain that de- 
excitation of this state involves the emission of the strong (quadrupole) y-rays 
of 1-12 and 1-76 Mev. energy in cascade (Rutherford, Lewis and Bowden 1933). 
The absolute intensities of these y-radiations have been shown to be closely the 
same (Latyschev 1947) and calculated values are such that an initial excitation 
of the 2-88 Mev. state of at least 20°% must be assumed. If the total energy of 
disintegration is 3-17 Mev., excitation of the state in question involves a B-particle 
partial spectrum of 0-29 Mev. maximum energy only, and the highest excitation 
possible (assuming an allowed transition) is then about 1:5°%. Onthe other hand, 
a partial spectrum of (0-29 +0-606) Mev. maximum energy would in fact give 
22:5°, or 67:5% excitation, depending upon whether curve O’ or O were 
appropriate to the case. Curve O’ would apply if the transition were one for 
which AJ= —1; on this matter all that can be said is that it is very unlikely 
that AJ is positive, though it might be zero. In any case, however, the indication 
is clear enough ; the total disintegration energy must be considerably greater 
than 3-17 Mev. ; thus the ground-to-ground state transition RaC.C’ is not observed. 

It is worth while discussing the unobserved transitions in this disintegration 
somewhat further—for there is clear evidence that there are many of them. 
Thus Rutherford, Lewis and Bowden (1933) detected twelve groups of long- 
range «-particles from RaC’—in this way providing evidence for the excitation, 
direct or indirect, of twelve states other than the ground state in the disintegration 
RaC.C’—and the same authors were forced to postulate the excitation of at least 
four other states in order to account even for the stronger lines in the y-ray spectrum 
of RaC. Now, we have already suggested somewhat more than 20°% excitation 
for the state of highest energy (2-88 Mev.), and the general shape of the f-particle 
spectrum (Constantinov and Latyschev 1941) is inconsistent with the direct 
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excitation of more than a few other states to any comparable extent (say about 
25°, for the state of 0-606 Mev., and roughly 50% for a state or states at about | 
2-2 mev. excess energy). So our statement may be made more precise : probably | 
not more than four of the sixteen excited states which are regarded as definitely | 
established are excited by direct B-transition RaC.C’ to more than 10%. Several | 
of these states cannot be excited to more than a few tenths per cent at the most. _ 
Assuming that in the ground state the nucleus *4;RaC’ has /=0 and even parity, | 
we thus conclude that there are probably other states, of higher energy than that | 
of 1-41 mev., which like the latter state (Fowler 1930) also have J=0. Some | 
of these may not be involved in y-ray emission to any appreciable extent. De- | 
pending upon whether AJ=0 or —1 for the transition exciting the 2-88 Mev. | 
state, bearing in mind the nature of the cascade y-ray emission from this state, | 
we are led to assign J=4 or J=5 to the ground state of RaC. Thus these un- | 
observed transitions must be at least third-forbidden (|A/| >4). 

The general features of the level scheme begin to emerge from this discussion, 
but clearly it would be unprofitable to continue it further until the results of 
experimental studies of B-y and y—y correlation provide more detailed informa- 
tion. Here it should merely be remarked that only three points are plotted 
for RaC’on figure 1 (ZH =3-173 Mev., 23 45% 7 H=15140-02mev4 402107, 
E=3-78 mev., <0:5%) : to have included all the “ unobserved” points would ff 
have been neither practicable nor very helpful at this stage. | 

It has already been stated that not more than two of the “‘ observed’ points | | 
of figure 1 lie below the line IIa’. One of these is the point for the 3-173 Mev. | 
transition RaC.C’ which has just been discussed at length, the other is the point jf 
for the ground-to-ground state transition “*Q/T12°4% Pb. This point almost — 
certainly belongs either to curve IJb’ or to curve I[L1a’—the accents reflecting the 
fact that since the final nucleus is an (<Vh) nucleus AJ is presumably negative 
in this case. On the basis of either attribution AJ = —3 (parity being unchanged 
if IIb’ applies and changing if IIIa’ is in question). Again, since the energy 
dependence of A is the same according to either assumption (compare (5) with jf 
Rox~(aZ/2)-*(p W)~4 which is applicable to curve [Ila’) we may, without further /f 
decision between the possibilities, construct a curve through the point under 
consideration and another curve representing values of A greater by a factor of | 
7 (2|AZ|+1, with |AZ|=3). These curves are shown on figure 1 and appropriately |} 
labelled. It will be noted that they intersect curve Ia’ in the general region of the 
point for RaC (£=3-173 mev.), thus no definite assignment of this point appears _ 
possible on the basis of the diagram alone. A similar remark is obviously true | 
of the unobserved point representing the 2-38 Mev. transition ThC’”.D. As 
regards the ‘“‘unobserved’”’ points for the ground-to-ground state transitions 
RaC’’. D, ThC’’. D and RaC. C’ (see above), it would appear from their positions | 
on the diagram very likely that these belong to transitions which are at least fourth- 
forbidden (see Richardson 1948). 

Figure 1, it will be remembered, was constructed by omitting the 23 points } 
representing f-transitions in elements having 87<Z<93 from the Sargent | 
diagram recently published. We may reconsider these transitions briefly now, | 
making the following general observations : 


(1) 7 of the omitted points would lie, on figure 1, either between curves | 
O and O’ or somewhat above curve O. 
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(2) 8 of them would lie between curves Ia and Ia’ (thus ‘t nappens that 
the evidence that these curves are significant is statistically stronger on the 
basis of the higher-Z group data than on the lower-Z group data) 


(3) 1 point (belonging to a rare mode of the disintegration UX,. Uy) 
would appear to lie on the line OO,,. 


(4) 4 points would lie between curves Ila and Ila’. 


(5) 1 point would be found somewhat below the line IIb’ —II]a’ (this 
is the only point plotted for MsTh,, the disintegration of which is obviously 
almost as complex as that of RaC but much less thoroughly studied). 


(6) The remaining two (“ unobserved ”’) points would appear to belong 
to transitions which are third- or fourth-forbidden (these are the points for 
the ground-to-ground state transitions UZ.U,, and *Npe ona). 


Clearly there is nothing in these general observations to suggest that our 
detailed treatment of the points included in figure 1 would fail if applied to the 
omitted points ; the effect of the change in average Z value between the two groups 
(from Z=82 to Z=91) in fact appears to be even smaller than might have been 
predicted. 

But there is one particular case which merits further discussion, for on the 
evidence at present available it seems that a real difficulty is involved. The 
disintegration sequence UX,;>UX,+U,, between the (°%") nuclei a UX, alid 
*92U, provides the case in point. These disintegrations have been very carefully 
studied by Bradt and others (Bradt, Heine and Scherrer 1943, Bradt and Scherrer 
1945, 1946), and, in spite of the complications due to the UX, —UZ isomerism, 
experimentally it would seem that the evidence is definite and in the main com- 
plete. But the positions of the representative points for the partial spectra of 
highest energy in the disintegrations UX,.UX,, UX,.Uy certainly lie on different 
Sargent curves on the revised diagram—the former between curves Ia and Ia’ 
and the latter on curve O. There is no evidence at present against the obvious 
assumption that the partial spectra in question characterize truly consecutive 
transitions, or against the assumption that the nucleus Uj, is left in the ground 
state after the second of these transitions, yet the Sargent diagram analysis points 
to a change of parity in the first transition (UX,.UX,) and no change of parity 
in the second (UX,.Uj;)._ The apparent overall change of parity as between 
the initial and final (¢y%") nuclei disclosed by this analysis is the difficulty already 
mentioned. 

There is another more subtle difficulty which concerns the disintegration 
UX,.Uj,; considered in isolation from the preceding transformation. If the 
two component radiations of about 0-95 Mev. and 1:5 Mev. energy identified by 
Bradt and Scherrer (1945, 1946) were the only y-radiations emitted by UX, 
the difficulty would not arise : observed intensities fit well with the assumption 
that for these radiations excitation of the y-emitting states (of U,,) is by first- 
forbidden transitions, “ cxcitation”’ of the ground state being by an allowed 
transition as above described. But there are two or three other radiations (of 
782, 822 and possibly 806 kev. energy, respectively) to be taken count of, and 
there is strong circumstantial evidence that some or all of these radiations are 
emitted in 0-0 transitions leading to the ground state. ‘The serious difficulty 
in this interpretation is that, if the emitting states and the ground state have the 

30-2 
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same spin ([=0) and parity, excitation of the emitting states by f-disintegration 
should be to the extent of about 20°% whereas in fact it is almost certainly less 
than 0:5°%. Even the assumption of different parities would not entirely remove 
the discrepancy. Empirically the representative points for the -transitions 
exciting the emitting states in question (only one such point was plotted in the 
Sargent diagram recently published) appear to lie on the curve OO, (0-0 tran- 
sition, no parity change), as mentioned in our general observation (3) above. 
This, of course, is consistent with the assumption that the subsequent radiations 
involve 0-+0 transitions to the ground state of the (<v) nucleus 73;Uy,: it is 
rather the main f-transition UX,. Uy, which does not appear to fit into the general 
interpretative scheme—it cannot similarly be described as a 0-0 transition 
between states of even parity in view of the position of the representative point 
on the allowed curve O of the revised diagram. It would seem, therefore, that 
unless our general assumption regarding spin and parity of ({yeh) nuclei in the 
ground state is at fault, some new feature of the UX,.UX,, UX,.U;, disintegra- 
tions remains to be discovered. More detailed studies of the known y-radiations 
in emission, and a search for a possibly highly converted low-energy y-radiation 
of high intensity, would appear to be called for. It might be mentioned here 
that no difficulty arises in assigning degrees of forbiddenness to the various 
6-transitions in the associated disintegration UZ. Uy, : energies and intensities 
of excitation, and the nature of the y-radiations, can be satisfactorily explained 
on the assumption that the /-transition of energy 0:45 Mev. is allowed, that of 
energy 1-2 Mev. second-forbidden and the unobserved ground-to-ground state 
transition of about 1:95 Mev. energy third- or, more probably, fourth-forbidden. 

There remains finally the case of the successive disintegrations MsTh,— 
MsTh,—RdTh, the last case amongst the “‘ classical”? radioelements of successive 
B-disintegrations between initial and final nuclei which are both (¥@}) nuclei. 
Here, as already mentioned, information is far from complete, but it is complete 
enough to sustain one clear prediction, so long as our basic assumptions remain 
valid. According to Lee and Libby (1939) the maximum energy of the f-particles 
of MsTh, is54+4 kev. On the basis of a single disintegration mode this implies 
a first-forbidden transition according to our revised diagram. On the other 
hand, the ground-to-ground state transition MsTh,.RdTh is at least second and 
probably more highly forbidden (see our general observation (5) above). ‘Thus 
we must anticipate that the ground-to-ground state transition MsTh,.MsTh, is 
unobserved, that is, that there is an intense (soft) y-radiation from MsTh, still to 
be discovered. 

Our discussion of the Sargent diagram for the heaviest radioelements has 
been based on the tacit assumption that the only selection rules involved are those 
concerned with the parity and total angular momentum of the nucleus—and |] 
empirically there is no indication that these rules are inadequate. With the |] 
lightest elements, however, there is abundant evidence that further rules are 
required, and Wigner (1939) and others have formulated such additional rules 
employing a “ super-multiplet ” classification of nuclear energy states. The 
empirical evidence for the necessity of some such extension of theory is obvious 
from figures 2 and 3. These diagrams refer, respectively, to electron-active and 
positron-active species having Z<22*. In the absence of additional selection 

* No point is plotted on figure 2 for *Li, the B-spectrum of which cannot be resolved into a 


finite number of discrete partial spectra because the disintegration energy varies continuously 
from zero to about 12 Mev. For a full discussion see Wheeler (1941). 
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tules we should expect the allowed curves to be well-defined, and that they 
should occupy identical positions on the two diagrams, the effect of the coulomb 
field of the nucleus being all but inappreciable for elements in this range of Z, 
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Figure 2. Ain sec-!; E in Mev. 
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Comparison of the diagrams shows, however, that these predictions are far from 
the truth. There is the appearance of an allowed line (A) in figure 2—if we leave 
on one side for a moment the points for the B-active species °H, SHe and “B— 
though it is not so well-defined as the allowed line for the heavy radioelements | 
in figure 1, but if this line is transferred from figure 2 to figure 3 (A’) it is seen to 
fall well helow the majority of points on that diagram. By contrast, in figure 3 
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Figure 3. insect; E in Mev. 
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there is a perfectly well-defined line (S) parallel to A’, which line we can only | 
classify empirically as “super-allowed”’. In this connection the treatment of 

Wigner is obviously significant: of the 15 points defining our super-allowed | 
line, 12 belong to the ground-to-ground state transitions of 12 of the 16 (477) | 
species “3C to 3}Sc (and there is fairly clear indication from the known periods that 


the corresponding points for three at least of the remaining four species of this 
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Sequence, viz. i1Na, {3Al and 3;Ca—only {jK is at present unknown—also fall 
on this pp ead the other three represent the ground-to-ground state tran- 
5 ( A Nett Bb a 26 / / 5 7S arte j 
sitions of *¢C, SF and 78Al. According to Wigner, only transitions belonging 

to the classes 
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with Z odd are not necessarily forbidden by any rules. It would appear that 
for Z <22 * all positron-active species falling into the first of these classes undergo 
super-allowed transitions (from ground to ground state) and that some of the 
species of the second class behave similarly, though others do not. No super- 
allowed transitions, however, occur for positron-active species which belong 
to neither class. Amongst the partially forbidden transitions of species of 
ee second class we might mention the ground-to-ground transitions of {jNa, 
isP, Cl and {§K—though apart from the Wigner rules the case of }{2P might be 
regarded as that of an ordinarily allowed transition (see figure 3). | 

Because practically all species of the type (*Z_{) with Z<22 are stable species, 
examples of super-allowed transitions are much less frequent amongst negative- 
electron-active bodies than with positron-emitters. In fact {H is the only 
representative of the first Wigner class, and $He, '{Be and '(C are the sole repre- 
sentatives of the second class, amongst the radioelements for which figure 2 has 
been drawn. As with the positron-emitters, 7H, as a class I body, undergoes a 
transition which empirically appears super-allowed {, whilst pile Of clascel lus, 
but '{Be and '4C are not, in the super-allowed group according to the diagram. 
In this connection it may be mentioned that although the very varied behaviour 
of the different class II species has yet to receive a satisfactory explanation on any 
theory, the fact that super-allowed transitions occur with some species of this 
class is unambiguous proof that for these species the Fermi-type interactions are 
inadmissible (Konopinski 1943, Wigner and Feenberg 1942). A final obvious 
remark is that the position of the point for '2B (figure 2) is entirely anomalous 
according to the Wigner rules. 

Summing up, we may say that the Wigner rules claim evident success in 
describing the otherwise inexplicably short lifetimes of the (°4,') positron 
emitters, but to be set against this is the very considerable concentration of points 
about what appears to be the ordinary allowed line A for negative-electron-active 
species in figure 2. There would appear to be some consideration which restricts 
the effect of these rules even in the class of the lightest elements to which they 
should most clearly apply. In this connection the exact lie of the line A of figure 2 
is of some importance. This line has been drawn purely empirically as the best 
straight line having regard to the experimental points. As transferred to figure 3 
(A’) it is seen to represent equally closely the slope of the super-allowed line S. 
Furthermore it may be shown to be the straight line which best approximates 
to the theoretical predictions of the Fermi theory. The dotted line F in figure 2 
gives the Fermi function f (0, Wo) (Konopinski 1943) fitted to A in the region of 
E=1mey. The theoretical line obviously appears as in many ways an even 
better representation of the experimental results than is our empirical line A; 


* The choice of this particular range of Z for figures 2 and 3 was made precisely because no 
species belonging to the first Wigner class are at present known for Z>21. 

+ The f-disintegration of the neutron should also be super-allowed: in this case interpolation 
from figure 2 would indicate a half-value period of 741 minutes. 
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thus the well-determined point for 35S is no longer significantly out of place, and 
the super-allowed points for 3H and ®He are more nearly the same distance above 
the line. Also, the anomaly regarding the point for 1B is clearly reduced. 
Having established the significance of the allowed line A in this way by com- 
parison with theory, it is of interest to discuss the exact lie of the empirically- 
drawn allowed line O of figure 1 ina similar manner. We can do this most easily 
by reference to line F, figure 2, which provides the normalization for the theoretical 
results. In the table we give the difference of the ordinates (A=logA(O) — 
logA(F), for E constant) of corresponding points on the two lines as a function 
of the disintegration energy E. he theoretical differences, A’(81), A’(93), 
defined by A’(Z) =log[f(Z, W,)] —log[f(0, Wo)] (see Konopinski 1943) are given 
for comparison *. Inthe last line of the table the quantity 6 is the second difference 


Table 
E(mev.) 0-:0158 0-0511 0-158 0:307 0-511 irO22 1-533 2:044 
log E 2-20 2708 1:20 1-487 1-708 0-009 0-185 0-310 
A 1-47 1:47 1-32 1:16 1:13 0:97 0:84 0-69 
A’(81) Dee, 223 1-82 ibeSyil 1-45 1-27 1:20 1:10 
A’(93) 2-91 2°58 2:14 1°81 i OUP 1-51 1-42 1:30 
5 125 0:94 0-66 0:50 0:46 0-42 0:47 0-51 


3[A’(81) + A’(93)] —A. The constancy of 6 over the range 0-2 Mev. <E <2-0 Mev. 
indicates that the straight line O, figure 1, has the full support of theory as 
representing the variation of A with £ over the range in question. ‘This conclusion 
may be stated conversely and expanded as follows : the Fermi theory of allowed 
transitions appears to be entirely adequate so far as predictions of the dependence 
of disintegration constant on energy and charge number are concerned. ‘The 
latter point is of interest, and the numerical value of 5 over the range of EF for 
which it is constant indicates that allowed transitions of the heaviest elements 
are less probable than the transitions which we have now with some confidence 
identified as ‘‘ ordinarily allowed” for the lightest elements (see lines A and F, 
figure 2) only by a factor of roughly 3 (for the same energy release). This is of 
importance in relation to the problem of the extent to which the wave-functions 
of initial and final nuclear states “‘ overlap’ for nuclei of such different degrees 
of complexity. 

The increase of 6 for E<0-2 Mev. which is evident from the table implies 
that the theoretical allowed line for the heaviest elements lies above the empirical 
straight line O of figure 1 at the low-energy end, the deviation being somewhat 
similar to the deviation of the theoretical line F from the straight line A of figure 2. 
‘This conclusion is not surprising and does not invalidate any of our previous 
discussions, for the purpose of which the empirical allowed line has been taken 
as reference line. Only with the heaviest elements it should be pointed out that 
for the extreme low-energy end of the line (H <0-02 Mev.) the predictions of the 
simple theory are obviously inadequate, since they do not take count of the energy 
derived from the “ adiabatic.”’ reorganization of the extranuclear electrons 
which for these elements is of the order of 0-01 Mev. 


* In this calculation we have falsified theory only to the extent of using Do, the maximum 


momentum, instead of p, an undefined mean momentum which enters as a slowly-varying term 
in the expression for f(Z, W,). 
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; There is just one feature of figures 2 and 3 which remains to be considered. 
We have so far concentrated attention on the identification of a number of transi- 
tions which appear to be allowed according to the spin and parity rules and for 
which there is no empirical evidence of the intervention of the Wigner rules. 
It will be obvious, however, that there are many points, particularly on figure 2, 
which do not come within the scope of this identification. Some of these may 
refer to transitions on which the Wigner rules impose increased forbiddenness— 
some of them, therefore, may be allowed points so far as spin and parity changes 
are concerned—others may represent transitions which are “ ordinarily’ for- 
bidden to first or higher degree. In this connection, since the condition expressed 
by (1) is no longer fulfilled, we should not expect to find an ordinary first-forbidden 
line parallel to the allowed line A. Also, since fairly large relative changes in p 
are involved, as between the lightest and the heaviest radioelements in the range 
Z<22, there is further reason why the higher-order Sargent lines should not 
stand out so obviously in the diagram. The question of order of forbiddenness 
must therefore be settled on other evidence for a number of transitions before 
much progress can be made in empirical classification. 

Partial evidence of one kind is already provided by the experimental results 
considered in this paper. Figure 2 includes the representative points for 13 
radioelements of which the f-disintegration is known to be complex, and with 
these species some transitions occur which would appear to be ordinarily allowed. 
In two cases of two-component spectra (j?Na and #4P) each of the observed com- 
ponents appears to result from such an allowed transition; in six other cases 
CEN, *9F, 73Mg, 73Al, Cl and {JA) one or more allowed components are present 
together with another component or components which are to some degree 
forbidden. It seems reasonable to assume that in these cases the other B-tran- 
sitions are ordinarily forbidden, that is, that the Wigner rules are of no more 
effect for them than they are for the allowed transitions with which they compete. 
If we make this assumption, we can at once identify the representative points 
numbered 6, 12, 19, 21, 27, 28 and 30 on figure 2 as ordinarily forbidden. On 
the further likely assumption that some of these points, at least, represent tran- 
sitions which on spin and parity rules alone are first-forbidden, we then have a 
general impression of the lie of the second Sargent curve on our diagram *. 
With this approximate curve in mind we can recognize a further group of points 
(e.g. 14, 23, 35, 38, 41) which lie between the ordinarily allowed and first-forbidden 
lines—and we are led to conclude that these belong to species (7{Na, ?{Si, 33Ca, 
3°Sc, 2°Sc) for the B-transitions of which the Wigner forbiddenness is greater, 
and more varied, than it is for the other transitions which we have been con- 
sidering. It is somewhat surprising to find the heavier species predominant 
in this group, but this is clearly the division to which our tentative considerations 
lead. 
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ABSTRACT. A detailed description is given of a semicircular focusing f-particle 
spectrograph in which a thin source is centrally placed and sheaves of oppositely moving 
electrons are investigated using two thin-windowed Geiger counters. ‘These counters can 
be used in a coincidence arrangement, and each is capable of movement so as to accept 
electrons which have been focused with any radius of curvature between 2:5 cm. and 6 cm. 

The arrangement was tested using sources of thorium active deposit. During this work 
the (coincidence) $-particle spectrum for the ThB.C transition leading to the 238 kev. 
state of ThC was obtained separately for the first time and its end-point determined as. 


355 kev. A feeble higher-energy B-component due to ThB was also detected and the 
end-point estimated at 589 kev. 


Redeterminations of the end-points of the B-particle spectra of ThC and ThC” (by, B-y- 
coincidences) gave values of 2:256 Mev. and 1-805 Mev., respectively. The intensity of the 
F-line of ThB was found to be 0-288 photoelectrons per disintegration. 


SS LN ERO DUCm TON 


HE coincidence counting method of investigating time-correlations in the 

emission of particles and quanta by unstable nuclei has already become one 

of the standard techniques in the study of radioactivity (Dunworth 1940, 
Feather 1940a). A development of this technique, first proposed by Feather: 
(1940 b), has now been realized in practice : in the present paper the method will be: 
described in detail and the results which have been obtained in its initial appli- 
cation to the active deposit of thorium will be presented and discussed. 

The essential feature of the new method is the study of coincidences between 
the various components of the particle radiations emitted by a radioactive source, 
e.g. between the internal conversion electrons (belonging to the B-ray “lines 
and the disintegration electrons (of the continuous spectra). For this purpose a 
special form of semicircular focusing spectrometer is employed, designed so that. 
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(for example) a f-line may be focused on one Geiger counter while another 
counter is used to scan the continuous spectrum. The spectrum of coincidences 
will then be the continuous spectrum of those disintegration electrons time- 
correlated with the photoelectrons of the f-ray line, and thus, by inference, with 
the y-ray quanta of corresponding energy. 

It is well known that magnetic analysis by semicircular focusing is inconsistent 
with the use of large solid angles for collection—in our particular arrangement (as 
shown in § 2) these are of the order of 10-?-10-2 depending upon position; on the 
other hand, the sensitivity of any coincidence experiment increases rapidly with 
the collection angles which can be employed. A natural development will 
therefore be to use the present method with other types of 6-particle spectrograph 
ci greater collecting power (see e.g. Roberts, Elliott and Deutsch 1944). Nothing 
new in principle would necessarily be involved in such a development. 


Se Eien sby Cal RO\G REAP En 


A detailed description of the spectrograph has already been given (leather 
1940b). The essential features of construction are illustrated in figure 1. Slits 
5, and S, select two sheaves of oppositely moving f-particles coming from the 
source A situated midway between the slits. [he two sheaves of electrons are 
focused by a magnetic field, perpendicular to the plane of the diagram, in the 
horizontal planes containing respectively the windows of the counters C, and Cy, 
and the defining slits. 

The counters can be moved independently in the planes of focusing thus 
allowing the radius of curvature of either sheaf of electrons to be varied between 
2:5cm.and6cm. Electrical connections to each counter are made by two spring 
strips attached to the counter carriage. One of the strips bears against the side 
wall of the box and serves to earth the cathode. The other strip makes a sliding 
contact with a copper wire insulated from the box and supplies the anode voltage 
(see also figure 8). 

A shutter B can be moved down to cover the slit S, and so to cut off the electrons 
from C,. Movement of the counters and the shutter is obtained by means of 
screws which can be rotated from outside the vacuum box by keys turning in 
ground joints. ‘The inside of the box is covered with aluminium and the slits are 
bevelled to reduce scattering as much as possible. The box is evacuated to a 
pressure of 10-*mm. Hg by a three-stage oil diffusion pump and is locked in 
position in the 5-5cm. gap between the poles of a large permanent magnet (Cock- 
croft, Ellis and Kershaw 1932). The use of a permanent magnet allows the 
magnetic field to be kept constant over the long periods necessary for coincidence 
counting and the provision of field coils enables the field strength to be set at any 
value up to 2000 gauss in a few minutes. ‘T’o avoid opening the box in order to 
measure the magnetic field a search coil is mounted on the top of the spectrograph 
box and the change of flux produced when the coil is rotated through 180° in the 
stray field of the magnet is measured with a fluxmeter. For standardization this 
change of flux can be compared with that measured with a search coil of known 
area-turns inside the box. An accuracy of better than 1 °% can be obtained in this 
way. still more accurate value for the effective field strength can be obtained by 
finding the position of one or other of the prominent B-ray lines in the spectrum 
of the thorium active deposit. ‘The Hp values of these lines have been measured 
with great care (Ellis 1932, Arnoult 1939, Siegbahn 1944) and since the position of 
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the counters can be determined to 0-01cm., the effective field strength can be 
obtained to an accuracy of 0:25°%% by dividing the accepted value of Hp by the 
observed value of p. Since several lines can usually be obtained for one field 
setting, a check on the accuracy of the measurements is generally possible. 

The profile of a B-ray line as observed with a constant field setting (figure 2) 
extends over an appreciable range of p. It is necessary to determine which point 
on the counting-rate—position curve is to be taken as fixing the position of the line, 
and what correction is necessary in order to obtain the true value of p from the 
value corresponding to the position so fixed. If in figure 3 (which is drawn in 
perspective) we consider the source Q as emitting electrons of one momentum 
only, P, then the trajectories in a plane perpendicular to the magnetic field will 
be circles of radius pp» = P/He. An electron from the centre of the source passing 
through O, the centre of the slit, will arrive at R on the centre line of the plane of 
focusing at a distance ay from O, where 4p)2=a,)?+6?. Taking axes X, Y at R 
as shown, we may calculate the intensity at a point (x, y) by a method similar to that 
employed by Li (1937), for a source of dimensions given below. 

The variation of intensity with y for fixed values of x is shown for a typical case 
in figure 4. If we integrate one of these curves with respect to y between the 
limits +/and —/, we obtain the intensity over a narrow strip of width dx and length 
21 for a fixed value of x, 2/ being the length of the counter window. ‘This is in 
effect the solid angle, expressed as a fraction of 47, subtended by the narrow strip 
2ldx at the source, for electrons of momentum P. ‘The variation of this solid 
angle with x is shown in figure 5. If the width of the counter window is w, we 
can then obtain the solid angle subtended by the entire window at the source (the 
transmission factor T(x)), when the centre of the window is at x, by integrating 
the curve in figure 5 between the limits x +w/2 and x —w/2. Figure 6 shows the 
variation of 7(«) with x, or, in other words, the variation of the counting rate for 
electrons of momentum P as the counter slit is moved in the plane of focusing in a 
‘constant magnetic field, i.e. the line-shape of a 6-ray line of momentum P= Hepp. 
Having regard to figure 6 we have chosen as the point fixing the position a of a line, 
the point A at half the maximum height on the steep right-hand edge of the 
counting-rate—position curve, since this point can be located with greater accuracy 
than the position of the peak. The calculated line-shape then enables us to 
determine the small correction term (less than 0-25 °%) which must be subtracted 
from a in order to obtain a) from which the value of py is calculated. 

The resolving power of the spectrograph, defined as Hp/8(Hp), where 8(Hp) is 
the change in the Hp value corresponding to the width of a B-ray line at half the 
maximum height, is obviously p/dp in the case of a constant-field spectrograph. 
‘This ratio can be determined from the calculated line-shape. It is found to change 
with the radius of curvature, from about 33 for the smallest value of p to about 60 
for the largest p employed. For the setting of the F-line in figure 2 the calculated 
value of 53-8 may be compared with the experimentally determined value of 56:5, 
in good agreement with the former. 

To map out a continuous spectrum the counting rate in counter (1) is deter- 
mined for a series of positions of the counter, first with the shutter B open, then 
with it closed. Differences give the counting rates due to B-rays alone. These 
differences are corrected for counting-rate losses (usually negligible for the 
continuous spectrum), and for decay of the source. In order to give a true picture 
of the continuous (momentum) spectrum, the corrected B-particle rate must 
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further be adjusted to refer to the same solid angle of collection and the same 
interval of Hp. Although this problem might have been treated analytically it 
was in fact solved empirically.* A value of Hp in the continuous spectrum was 
chosen and observations were made with various field strengths H,, H,.... at 
radii of curvature p,, py.... such that H,p,=constant. Let the corrected B- 
particle rates under these conditions be N,, N,.... Now, for a constant radius 
.of curvature and a variable field, it is known that the normalizing factor is 1/H 
(Gurney 1925), thus under the conditions of our standardizing experiment the 
counting rates at p,, py.... for a constant field and a continuous (momentum) 
spectrum of uniform height would be as N,/H, : N,/H,.... This means that 
H,/N,, H,/N2.... are respectively proportional to the normalizing factors corre- 
‘sponding to p,,p2..-.. Plotting these quantities as a function of p, a normalizing 
‘curve is obtained for the particular counter used. This curve could be checked 
by the method described in § 6. 

Normalization of the continuous spectrum having been effected, it is still 
necessary for many purposes to know the solid angle for collection to which the 
spectrum has been normalized. ‘This may be calculated from curves similar to 
that shown in figure 5. For a radius of curvature of 5-8cm., for example, and a 
counter window of area 1-9 x 0-10 cm?, we obtain 3:4 x 10-3 of 47. To be precise, 
this fraction represents, in the conditions specified, the probability that an electron 
of the continuous spectrum will enter the counter, the probability being averaged 
over the momentum range Hp to Hp +A(Hp), where Ap is the difference between 
the radii of curvature corresponding to entry by the extreme edges of the counter 
window. During the experiments described below entrance windows of breadth 
from 0-10cm. to 0:15cm. were used at various times. 

The last problem in dealing with the continuous spectrum is that of determining 
the ‘‘representative”’ p value for any setting of the counter. This was solved by a 
method similar to that used by Lawson and Tyler (1940) for the same purpose. 


$3. THE SOURCE 


Sources of thorium active deposit were obtained by the exposure of aluminium 
foils of 0-28 mg/cm? thickness to the thoron given off by a “dry”? emanating 
preparation of radiothorium, of strength approximately 40 millicuries. For 
activation a foil is bent over a suitably formed brass block in an arrangement which 
allows an area of 19 x 2mm?* of foil only to be exposed. Both short- and long- 
exposure sources were prepared in this way. With long-exposure sources, 
transient equilibrium having already been attained, the ratio of the absolute 
activities of ‘hC and ‘ThB remains constant at A,,/(A; —Ag) during the subsequent 
‘decay ; with short-exposure sources this ratio varies with time. The numerical 
value of the constant ratio is 1-11. Important for our experiment is the fact that 
with a source in transient equilibrium 1-11 is very closely the ratio of the numbers 
of electrons in the continuous spectra of Th(C +C”) and ThB, respectively, the 
disintegration constant of ThC” being so much greater than the disintegration 
constants (Ap, A;) of ThB and ThC. 


* Attempts to calculate the normalizing factor in our arrangement are necessarily difficult on 
-account of the form of the counter window (§ 4). In general the electrons do not strike the plane of 
focusing normally, but arrive at angles to the normal spread over a wider range, and more inclined 
the smaller the value of p. Thus a fraction of the electrons entering the counter slit (which fe 
finite depth) strikes the far wall of the slit and this fraction increases with decreasing p. Some of 
these eJectrons, after scattering, may produce a discharge in the counter but many of them will not 
‘Obviously the effect 1s not susceptible of exact evaluation except by experiment. 
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With short-exposure sources information concerning the relative contributions 
of the B-particles of ThB and Th(C +C’) to any given interval of the experimental 
continuous spectrum can be obtained by the observation of the growth (or decay) 
of the B-particle counting rate when a source is used which has just been removed 
from the exposure pot and the Geiger counter is adjusted to accept the appropriate 
momentum range in the spectrograph. If we consider a source which has been 
exposed for one hour only, then, whilst the activity due to ThB should decay as 
shown by curve B, figure 7, the activity due to Th(C + C”) should follow curve C 
(calculated using standard equations). If an observed activity curve, such as A, 
lies between the curves C and B then the fraction of the activity due to ThB is given 
by the ratio XZ/YZ. Clearly, for a particular setting of the counter, this ratio 
should be constant within the limits of experimental error for any vertical secant to 
the curves. In order to achieve this result in practice it has been found necessary 
to wait for about half an hour after the end of an exposure before starting activity 
measurements of this type. It would appear that during activation atoms of 
ThC", given off by recoil from the active deposit which invariably accumulates on 
the walls of the pot, are collected together with ThB on the foil. If this additional 
ThC" activity is not allowed to die out before measurements are begun, the ratio 
XZ/YZ will be found to vary along the curves. 

In general, in order to prevent the escape of recoil atoms from a source in the 
spectrograph, the activated foil is covered with a similar aluminium foil (0-28 
mg/cm? thickness) before being clamped in the source holder. In position, the 
active deposit, sandwiched between the two thicknesses of foil, is inclined at an 
angle of 45° to the plane containing the central axes of the two defining slits and is 
symmetrically placed with respect to this plane. In this connection, the shape of 
the F-line of ThB was determined first with a bare source and then with a “ sand- 
wich” source, but no detectable difference was observed either in the height or 
shape of the line. It therefore appears that the effect of absorption and straggling 
produced by the 0-28mg/cm? covering of the source is negligible—at least for 
6-particles of energy as great as 147 kev. 


§4. COUNTERS 


‘The Geiger counters are necessarily of small dimensions since they have to be 
capable of movement in the restricted space between the magnet pole pieces. 
While this means that the advantages of a relatively low background counting rate 
(5-10 min“) and a short recovery time (30 usec.) are enjoyed (Muehlhause and 
Friedman 1946), at the same time it implies that at high counting rates the life 
of a counter is short (counting rates as high as 60000 min“ had to be sustained for 
several hours when a counter was registering a strong photoelectron line, or was 
placed close to the source to count «-particles). 

In figures 1 and 8 are shown the details of counter construction. ‘The cathode 
is made of brass suitably finished to provide a slit 1-0-1-5 mm. broad for the 
entrance of the f-particles. End pieces are of perspex and are threaded with 
a central wire of platinum of 0-1mm. diameter as anode. ‘The thin mica film 
(~1 mg/cm?) which is used to cover the slit is mounted on a concave surface in 
order to withstand better the atmospheric pressure to which it is subject during the 
filling process. A small amount of lead screening, lcm. thick, is mounted to- 
gether with the counter on a carriage of ebonite. 

A mixture of 10°% ethyl alcohol and 90° argon at a pressure of 10cm. Hg was 
used to fill the counter after the cathode surface had been washed in alcohol and 
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polished with silk. The filling line was simply constructed, use being made of 
rubber connections and vacuum stop-cocks, but the counters were evacuated to 
10-4mm. Hg with a diffusion pump. Otherwise the technique employed was 
similar to that described recently by Regener (1947). Satisfactory “‘ plateaux”, 
200-300 v. long, were obtained with threshold potentials ranging from 900-1100 v. 
and positive plateau slopes of about 0-10-0-15% per volt, though this latter 
quantity was found to increase somewhat when the counting rate was high. 
Another troublesome effect observed during periods of rapid counting was a 
gradual increase of threshold voltage (up to 20 or 30 v.) which occasionally made it 
necessary to use a radium standard to check the counter sensitivity throughout the 
course of an experiment. 

The importance of obtaining an accurate correction for absorption and scatter- 
ing in the mica window cannot be over-estimated since it becomes very large for 
energies below 70kev. (Hp~900 gauss-cm.). Additional thin sheets of mica 
covering the counter slit were therefore used to plot absorption curves for electrons 
of a number of energies in the thorium spectrum and in this way, by extrapolation 
of the results to zero thickness, the correction curve given in figure 9 was obtained... 
This curve shows reasonable agreement with that of Martin and ‘Townsend (1939), 
though it should be emphasized that the magnitude of the correction factor for 
electrons of any energy depends on the precise arrangement of the counter slit 
and of the mica covering the slit. 


$5.) BLE CARI CATE ERGO laes 


In any coincidence experiment it is necessary that provision should be made 
both for determining the individual rates of entry of particles into the two counters 
and for counting the coincident pulses which may occur as a consequence of the 
association in time of these particles. 

Figure 10 illustrates the basic circuit used for single channel counting. All 
H.T. voltage supplies are stabilized. It will be noted that a load resistance of 
1 MQ is used in the counter circuit, but even so, with the small output pulse 
obtained with a counter of small volume, it is necessary to use a conventional 
cathode-follower unit mounted close to the counter. The cathode-follower 
operates through a screened cable into an amplifier designed for the equalization 
of the pulses (the higher the counting rate the larger the proportion of pulses which 
are smaller than normal). ‘The particular equalizing circuit used is a modification 
of the cathode-coupled multivibrator (West 1944). This has proved to be 
exceedingly stable and reliable in operation. 

_ The 50kQ rheostat in the screen grid lead of valve V, is used to obtain suitable 
conditions for operation from triggering pulses, and can in fact be varied over a 
large range. ‘The output pulses, of 5 usec. duration, from the screen grid of Vs are 
used to operate a scaling circuit having a factor either of 100 or, for individual 
rates greater than 50000 min-’, of 3200. The Dynatron Model X200 “ Scaler’’, 
with its decimal counting system, has proved in every way satisfactory and 
extremely convenient in use. ‘The input resolving-time of this scaler is 5 ysec., 
thus it is evident that the 30sec. recovery time of the Geiger counter alone 
determines the reso]ving-time of the counting system as a whole. 

It was found that, even at the highest counting rates used (100 000 min), 
both the size and shape of the output pulses from the amplifier remained constant. 
Thus it can be taken as established that the resolving-time for coincidence counting 
was independent of the single counting rates throughout the entire investigation. 
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The output to the coincidence mixing set is taken from the anode of V3 in each 
amplifier. A conventional Rossi coincidence circuit is used, with a fixed grid bias 
on the output valve but with a choice of grid leaks (2000 Q and 20000 Q) for the 
Rossi valve. In this way two values of the resolving-time, 2:5 x 10-7sec. and 
1:8 x 10-8sec., are available. The requirement of the highest possible resolving 
power (smallest +) has already been stressed (Dunworth 1940). However, 
attempts to reduce t much below 2 x 10-7 sec. do not succeed, because of the statis- 
tical time-lag between the entry of the ionizing particle and the appearance of the 
voltage pulse across the counter (Dunworth 1940, Bradt and Scherrer 1943 a); 
instead genuine coincidences are lost. The larger of the two resolving-times 
available is used for purposes of checking that there is no loss at the smaller + and 
also when delayed coincidences are in question—as when the f-particle of ThC is 
followed by the «-particle of ThC’ after a mean time interval of (3:740-6) x 10-? sec. 
—the average life of ThC’ (Bradt and Scherrer 1943 b). 

A separate recorder was not provided for determining the coincidence rate, 
although this was only of the order of a few counts per minute ; instead, registration 
on the same scaler as was used for the single counting rates was arranged without 
inconvenience. 


§6. EXPERIMENTS WITH THE ACTIVE DEPOSIT OF THORIUM 


Thorium active deposit sources were chosen for the initial tests of the new 
method, and of our spectrograph generally, for two reasons : first because there is 
present in the complete £-particle spectrum an intense photoelectron line (F-line) 
eminently suitable for coincidence studies, and secondly because so much inform- 
ation is already available * regarding the radiations in general. In the course of 
our tests we have, however, been able to add to this information in certain parti- 
culars, details of which will be apparent from the descriptions which follow. 


(i) Investigations using a single counter 


A study was first made of the complete continuous f-ray spectrum of the 
thorium deposit in transient equilibrium. Because this spectrum extends as far 
as Hp = 9000 gauss-cm., it was necessary to make measurements at four different 
values of the magnetic field (250, 500, 1000 and 2000 gauss). From these the 
distribution curve shown in figure 11 was obtained after the appropriate normali- 
zations (see § 2), corrections for decay and absorption corrections (see § 4) had been 
made. ‘The fact that a smooth curve resulted from such a series of overlapping 
junctions provides a good test of the accuracy of the method of normalization. 

A further test is provided by a comparison of the areas of the ThB and 
Th(C+C’”) spectra deduced from the curve (the justification for using the dotted 
projection of the (C + C”) curve to separate the two spectra is given below). These 
areas, which should be proportional to the numbers of ThB and Th(C +C”) atoms 
disintegrating per second, were found to be in the ratio 1 : 1-12 in good agreement 
with the ratio 1 : 1:11 quoted in §3. 

The continuous 6-spectrum of the ThB having been obtained, the absolute 
intensity of the strong F-line of ThB may be calculated. The area under the 
F-line profile (§ 2, figure 2) represents the number of photoelectrons emitted per 


* Gurney 1926, Ellis 1932, Sargent 1932, Sze Shin Yuan 1933, Henderson 1934, Wang 1934 
Surugue 1937, Arnoult 1939, Flammersfeld 1939, Siegbahn 1944, Siegbahn and Johansson 1946. 
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second through a known solid angle, hence, adjusting the solid angle for collection 
to be equal to that assumed in normalizing the continuous spectrum, we have 


area under F-line profile _ | 599 


number of F-electrons per disintegration = 
area of b-spectrum 


This result agrees more closely with that of Flammersfeld (0-284) * than with the 
results of Arnoult (0-222) or Gurney (0-250). 

For the end-point of the ThC spectrum a Fermi plot of the single counting 
rates for Hp >7500 is given in figure 12, where (N/F)? is plotted with F calculated 
for Z=83. From the point of intersection with the energy axis we obtain 
Evnax = 2°256 Mev. (Hp =9070 gauss-cm.), in good agreement with the result of 
Henderson (1934), Hp =9040. 

The end-point of the ThB spectrum cannot similarly be determined by 
simple extrapolation. It will be shown in (i1) below how the main component of 
this spectrum (associated with the yp transition) can be separated out. However, 
it appeared from figure 11 that a component of higher energy might also be present, 
and experiments with a single counter were made to test this assumption. A series 
of measurements was performed of the growth and decay of the f-activity at 
different points of the spectrum, using sources exposed for one hour only. We 
have already discussed (§3) the way in which the activity of such a ‘‘non-equili- 
brium”’ source should vary after a “‘short”’ exposure (see figure 7). 

With the counter at a setting corresponding to Hp =3500 gauss-cm., experi- 
mental points were obtained showing the presence of a pure Th(C+C’) activity 
(figure 7, C). At Hp=2500 gauss-cm., on the other hand, the growth of activity 
was less marked (figure 7, A), indicating the presence after 50 minutes of an effect 
due to ThB amounting to 12+3% of the total. It can be shown that this effect 
would correspond to 7 +2% of ThB activity once transient equilibrium had been 
reached (after 50 minutes the ThC in the source is still growing, while the ThB 
is of course decaying in amount). ‘Thus in figure 11 we can identify the point X 
with certainty as belonging to the (C+C”) spectrum. To complete the low 
energy end of this spectrum (Hp <3000 gauss-cm.) the distribution given by 
Flammersfeld (1939) has been used. 

We have just seen that the ThB spectrum must have a component extending 
beyond Hp = 2300 gauss-cm., the currently accepted end-point. It is not difficult 
to imagine that this weak high energy component, with apparent end-point about 
589 kev., corresponds to a ground state-ground state transition between ThB and 
ThC. The maximum energy of the partial component associated with the y, 
transition is shown in (i1) to be 355 kev., and the energy of the F y-ray (emission of 
which is believed to lead to complete de-excitation of the nucleus) is 238-2 kev. 
‘The numerical check (355 +238~589) is fairly convincing.+ 


* Flammersfeld’s value 0-284 was deduced on the assumption that the areas of the continuous 
spectra of ThB and Th(C+C”) should be equal. He apparently did not regard the difference of 
35% in area which he found as significant. If his source was in transient equilibrium there should 
have been a difference in area of 11% as stated above, and in that caes values of 0-289 and 0:310 
would have been obtained if he had referred his F-line intensity to the observed areas which he 
assigned to ThB and Th(C+C”) respectively. 

} We have to thank Dr. H. O. W. Richardson of this department for drawing our attention to his 
work on the ThB spectrum. His results, as yet unpublished, and those of Flammersfeld (1939) 


both clearly show the existence of the high-energy component, though it is remarkable that 
Flammersfeld does not comment upon this feature of his curve. 


2-6 || 
al 
g 
aes 
as rm 
: a * 
oa LS 
33) oo 
eo 2 
3 ogi 
eile oe 
1S) 18 3 & 
L & 
. ) 
1-0 1 Ripe S) 0 a Say 
500 1000 1500 2000 2500 600 1000 1500 2000 
Hp (gauss-cm.). Hp (gauss cm.). 
Figure 9. Correction curve for mica Figure 13. Spectrum of coincidences 
window absorption (1:3 mg/cm?). of disintegration electrons with 
electrons of the F-line. The 
vertical lines represent probable 
errors. 
| VWV— 
At Pec ia 
+300v 5 | 
e ee ey 
cale 0} D ' 
call "3200 ne ae 
Rossi Stage oa ' 
re H 
8 ot \ 
iE 
EDS 
>) : 
: \ 
ith ae 
C6 5 : ome 
Ss 
zs) 
(a4 


oe | ane 
2000 4000 6000 8000 10000 


FH[p (gauss-cm.). 
Figure 11. The momentum distribution of 


= Vs 120r-- 
Figure 10. Single channel circuit. 


COMPONENTS the disintegration electrons of the active 
Resistances deposit of thorium in transient equilibrium. 
R1 11082 =R67:5x1083@ R11 10 10°Q Ordinates normalized to equal momentum 
ROO 5 G10 mR et OOK 105 SeRI2 a2 102 intervals and corrected for counter window 
R3 10 x 108 R8 320103 R13 220 103 absorption. 
R4 50x 108 ROTOO< 102 R14 ex 1102 
R5 100 ROBO 102 RISES 0a 02 
Condensors Valves 
C1 100 ver C5 1000 upr V1 
C2 150 upyr C6 8 uF V2 >EF 50 
C3 25puuF C7 8uF V3) 


C4 50 uur C8 8uF 


8 
5 
3 40 
z 
ad oak 
= 2-0 
J S ; 
—— i ’ 0 
4546 48 50 52 5-4 UO) ie aie sl 2:5 3-0 3-5 4:0 45 
(1+n*)3. (1-+n%)s (14-7?) 
re 12. Fermi plot giving ThC Figure 14. Fermi plot of Figure 15. Fermi plot of B-y coincidences 
end point. points given in figure CNN” spectrum), 
13, F calculated for Probable errors are shown. 


Z=82. 
31-2 


476 N. Feather, 7. Kyles and R. W. Pringle 


(ii) Investigations using the coincidence technique 


-photoelectron coincidences. ‘The number of genuine coincidences to be 
expected between photoelectrons entering counter (1), when it is centred on a 
B-line, and disintegration electrons entering counter (2), as it is moved through 
the continuous spectrum, may be calculated as to order of magnitude as follows. 

Let photoelectrons enter counter (1) at the rate of N, min. Then of the 
corresponding disintegration electrons a fraction w,F(P)AP will enter counter (2) 
where 47, is the solid angle subtended by the counter window at the source, 
F(P) is the momentum distribution function describing the partial spectrum for 
the disintegration mode with which the f-line is associated, and AP is the range in 
momentum, about the mean momentum P, of the particles which enter the counter. 
The value of ws, as we have seen in § 2, may be assumed to be of the order of 4 x 1073, 


and since F(P) must be normalized so that ie i(P)dP=1, PP, being the upper, 
limit of the (partial) momentum spectrum, F(P)AP may be of the order of 1/50. 
The coincidence rate to be expected is then Ng=N,w.F(P)AP ~ 4 min", if we 
take N,, as 50000min“'. This value of N, could conveniently be attained in our 
arrangement for the F-line of ThB. 

When counter (1) was fixed on this line, coincidences were observed with 
counter (2) at various positions in the continuous spectrum. A coincidence 
resolving-time of 2:5 x 10-* sec. was used, and a typical set of results obtained near 
the momentum peak of the main ThB spectrum reads : 


N, (continuous spectrum +y background) = 7000 min“! 
N, (F-line + continuous spectrum +y background) =40000 min 
N42 (coincidences observed) =5-1 min‘? 
N, (accidental coincidences = 2N,N,7) =2°3. min? 
Ng (genuine coincidences) =2-8 min 


It may be noted that with a counting period of one hour the accuracy of 
determination of Ng is only about 8 % in this case. 

An extended series of observations was carried out with the object of deducing 
the complete momentum spectrum of coincidences. The usual corrections were 
made—for decay, for absorption in the counter window, and for the variation 
of the effective solid angle of collection with the position of the counter (see § 2)— 
as also corrections for (yy) and (6-y) coincidences (<5 % of the total even in the 
most unfavourable case). ‘There was no necessity to correct for cosmic-ray coinci- 
dences, which were negligible throughout. After all corrections had been made, 
the final spectrum of coincidences was as shown in figure 13. Two field strengths 
(H = 280 and 440 gauss) were used to obtain this spectrum. 

No genuine coincidences with the F-line were detected beyond Hp =2350 
gauss-cm., and the observed points were found to fit a Fermi plot (Z =82) very well 
(figure 14). ‘The Fermi end-point appears at (1 +7?)'=1-695 or E,yjay =355 kev., 
the value used above. 

An increase of resolving-time to 1-8 x 10-8 sec. was found to be without effect: 
on the number of genuine concidences observed, thus we may conclude that the 
lifetime of the 238 kev. excited state of ThC cannot be greater than 10-7 sec. at the 
most. 

While these results are to be regarded as satisfactory, the method obviously 
allows of much improvement in the direction of shortening the periods of counting 
which are required under conditions of great stability of operation of counters. 
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and circuits. A limiting factor at present is undoubtedly in the Geiger counter 
itself, the essential characteristics of which determine the smallest resolving-time 
which can profitably be employed. If the electron multiplier could be used in 
place of the counter, the scope of the method would clearly be greatly extended ; 
this change would probably enable 7 to be reduced by a factor of 10? and the loss of 
accuracy involved in the background of accidental coincidences to be largely 
avoided. 

B-y coincidences. In an experiment similar to that of Bothe and Maier- 
Leibnitz (1937) one of the counters, shielded so as to record only y-rays, was placed 
as near as possible to the source (dotted position in figure 1) and another counter 
used to traverse the continuous spectrum as before. Coincidences were observed 
for electron energies above 650 kev. (Hp >3500 gauss-cm.). . 

Since the disintegration electrons of ThC are followed by only a small amount 
of medium-energy y-radiation, the spectrum of coincidences which is obtained in 
this arrangement will follow very closely the B-spectrum of ThC” (which emits a 
2:62 Mev. y-ray in every disintegration). The Fermi plot (figure 15) of the 
experimentally determined coincidence spectrum is a straight line cutting the axis 
at E,.4x =1:805+0-02mev.* Henderson (1934) obtained the value 1-795 Mev. 
using a pure ThC” source by a single counting method in a magnetic spectrograph ; 
the only other published determination with such a source is that of Chalmers 
(1929) who obtained 1-75 mev. by an absorption method. ‘The coincidence 
resolving-time used in our experiments was 2:5 x 10~’sec. 

a—B coincidences. "The «-particles from ThC’ (8-6cm. range) follow the dis- 
integration electrons from ThC at a mean interval of (3-7 +0-6) x 10~“sec. (see 
above). Thus if an «-particle counter is set up near the source, with an 8 mg/cm? 
sheet of mica (to cut out the unwanted «-particles from ThC of 4:7cm. range) 
replacing the lead screen used in the last arrangement, it will register coincidences 
with the disintegration electrons of ThC if the longer resolving-time (1-8 x 10~® sec.) 
is used. In this way the total 6-spectrum of ThC may be mapped out once the 
contribution of the 6-8 { and y—f coincidences has been assessed (<10°%, of the 
total). Because of the increased resolving time the results obtained with the 
present arrangement were not good enough to reproduce accurately the shape of 
this spectrum, but they established the fact that genuine «8 coincidences could 
be recorded up to the known end-point (9070 gauss-cm.). 

Photoelectron-photoelectron coincidences. ‘The B-spectrum of the thorium 
deposit contains no suitable pair of lines for the study of a cascade process, but the 
work is being continued with the active deposit of radium in the B-spectrum of 
which there are a number of f-lines (e.g. A, B, F, G and H) all of sufficient intensity 
to allow of such an investigation. This work will be reported on in a later com- 
munication. 
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Note on the Impedance Variations of an Electro-Acoustic 
Transducer in a Reflecting Field 


By -S. BYARD 
Royal Naval Scientific Service, Department of Physical Research, Admiralty 


Communicated by A. B. Wood ; MS. received 2 March 1948 


HEN a reflecting surface is introduced into the field of an electro-acoustic 
transducer so that a certain proportion of the radiated energy is returned 
to the source, the measured electrical impedance of the transducer will 

vary according to the relative phasing of the returned energy. Consequently if the 
reflecting surface be moved towards or away from the transducer, or alternatively, 
if the frequency of the electrical input be varied, periodic fluctuations may be 
observed in the input voltage or current. The Pierce interferometer uses this 
effect in determining acoustic wavelength in a liquid or gas. The purpose of this 
note is to relate the impedance changes responsible for the voltage and current 
fluctuations to the efficiency of the electro-acoustic device by making use of the 
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reciprocity principle. It is also indicated how, in theory at least, the absolute 
calibration of a projector or a receiver may be obtained by impedance measure- 
ments on a single instrument, and how this is related to Carstensen’s more practical 
“self-reciprocity’’ method of determining transducer. constants. 

Let Ey be the open-circuit voltage response of an electro-acoustic transducer 
when placed in a sound field of pressure Pe~*”. If M is the receiver sensitivity 
expressed in open-circuit volts/dyne/cm?, then 


(RE tela cae ee he (1) 


Let Z, be the nominal electrical impedance of the transducer. Suppose now that 
a voltage E, of the same frequency as E, be applied to the transducer terminals ; 
the equivalent circuit is shown in figure 1, and the resulting current J is given by the 
expression 


Bel Zig gn ae AO haar (2) 
The measured impedance Z under these conditions would be 
Vice WOW A wed Oa Vs SO) Tu i OR ce (3) 
or from (1), expressing E in terms of P, 
Za diy PeyLe eg Meg ee Fe (4) 
| xX 
Tot | | 
| : | - 
ee E <—d |——+ Fi 
be | 
i : 
Figure 1. Figure 2. 


The field P may be assumed to be produced by an identical transducer excited 
by a current J at the same frequency, and situated at a distance d from the first 
transducer as in figure 2. If S, is the pressure produced at this distance d per 
ampere excitation, then 


P= Sie Sahl ee a are (5) 
and substituting this value for P in (4), | 
; a Ly eM Sie C9 eae eee eee (6) 


The factor MS, in (6) above is the change in impedance due to the sound field 
of the second transducer and is added vectorially to Z) in order to obtain Z.  Z will 
have a maximum whenever MS, is in phase with Z), and a minimum when MS, 
is 180° out of phase, and the impedance will show periodic fluctuations of diminish- 
ing amplitude as the transducers are moved away from one another. 

Let two impedance measurements Z, and Z, be made at an adjacent maximum 
and minimum when the transducers are at distances d, and d, apart respectively 
(d,=d,+A/2, where A is the wavelength). Then if S, and S, be the appropriate 
transmitting constants obtaining at the two positions, 


[Z,|=|Zo/+|M@Si13  [Z,|=|25|-| 4482], era) 
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and the difference between these two moduli is 
|AZ|=|2Z,|—|2Z2|=| (Si + S3)]. Bee. (8) 


The reciprocity parameter J which relates the transmitting and receiving responses 
of a transducer is given ae 1926, Maclean 1940) by 


J =| M]S,|=(2Ad/R) x 10-7, 


where A is the w avelength of the sound and R the characteristic acoustic impedance 
of the medium in c.c.s. units. If this is assumed to be valid for the system, it 


follows that 
uaa 
die td. 


Now suppose the second transducer eee by a perfect reflector, XX in figure 2, 
le 
which can be moved to satisfy the required impedance conditions. If d’ and d 
are the two distances of the reflector from the source, satisfying the two impedance 
conditions, then equation (10) becomes 
M?R 


ue ed 
EE IE 2rd 


where d is the mean of d’ andd”. In this equation the distances between source 
and reflector are doubled to give the total lengths (d, and d, in (10)), of the trans- 
mission paths. The approximation is permissible when d >A, and this will 
obtain at those distances where the reciprocity parameter defined in (9) is valid, 
and where the inverse relationship between pressure and distance is established. 


From (9) and (11), the transducer constants may be expressed in terms of 
AZ, Ry Aand ds 


| M|=[(2Ad/R)|AZ|x10-7}#; | S,|=[(R/2Ad)|AZ|x 107}... (12) 


An examination of these expressions with reference to the performance of 
practical transducers shows that their application to the absolute calibration of a 
transducer would be limited to efficient projectors functioning at resonance, and 
that in any case the accuracy will be low since the result depends upon the measure- 
ment of a small change in impedance. Carstensen (1947) has, however, recently 
described a pulse technique in which the transmitter current J and the’ received 
voltage Ey are measured separately. A train of waves is used short enough to 
ensure that the transmission does not overlap the reception of the echo. Refer- 
ence should be made to the original paper for practical details and for a discussion 


of the limitations. The expressions for the transducer constants | | and | S,| 
are, in terms of E, and J, 


pve 


107 eee (10) 


|AZ|= x 107 & 


10h. 8 eae (11) 


| M|=[(4Ad/R)(Ey/Z) x 10-7}; | S,|=[(R/Ad)(Ep/Z) x 107]... (13) 


and it will be observed that the ratio Z)/J in these expressions is the added imped- 
ance of equation (3) and is equal to half the differential impedance AZ. 
Thanks are due to the Admiralty for permission to publish this note. 
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Dependence of @-Decay on the Sign of the 
Change in Nuclear Spin 


The calculation of the probability of B-decay for given values of the magnitudes of the 
initial and final angular momenta of the nucleus may be abbreviated, as has been done by 
Marshak (1942), by application of a sum rule and of conservation of angular momentum. 
Professor Feather and Dr. Richardson have pointed out to me a discrepancy between 
Marshak’s results and a statement of Greuling (1942). Here Marshak’s argument is 
amplified and Greuling’s statement shown to be wrong. 

Under certain approximations the interaction between the nucleons and the electron- 
neutrino field contains terms of the type Vy, Vp, where Vy, Vy are respectively functions of 
the nuclear variables only, and of the field variables only. Let P(j, m:j’, m’) denote the 
probability of transition of the system from a state where the angular momentum of the 
nucleus corresponds to the usual designation 7, m to the state labelled by 7’, m’, the pro- 
bability being summed for all consistent values of energy, linear momenta, angular momenta 
and spin momenta of the electron-neutrino field. ‘Total angular momentum is conserved to 
the usual approximation. Since /y is Hermitian, and since the part of the transition 
probability arising from Vp is independent of the sense of the change of angular momentum, 


TEX GO 34) pO NSIAG a Sp Ht SS (1) 


Define O(j, 7’) to be the probability of the transition from a single m belonging to j to any of 
the 27-+-1 values of m’ belonging to j’: thus 


j' 
Oni = > PG, HTS) ae eee (2) 
m=—j 
‘That O(j,j’) is, as the notation suggests, independent of m, may be seen as follows, ‘The 
usual expansion of the light particle wave functions (evaluated at the position of a nucleon) 
in ascending powers of the radius (see for example Konopinski and Uhlenbeck 1941) means 
that eventually the operator Vy contains, among other things, coordinate vectors of the 
nucleons: if %, %, are wave functions of electron, neutrino the term ¥*V py, contains 
operators (in Vy) not involving spatial coordinates, together with spherical harmonics, etc. 
When an integration over the azimuth ¢ is performed, the integrand contains the azimuth in a 
factor of the form exp (:M¢) and in terms involving coordinate vectors of the nucleons. 
This is now entirely similar to the theory of light emission, and a similar sum rule (see for 
example Condon and Shortley 1935) is valid. 
Summation of (2) for all values of m belonging to j gives, on account of (1), the usual 


result : 


Zi WOG TF) SCI WOU CI). ee ae (3) 


Marshak’s simplification. Let k, » specify the angular momentum of the electron- 
neutrino field so that, using suffixes i, f to denote initial, final values, 


ie ays, a eect (4) 
RS Ss HO BEC (5) 
Suppose first that j;—jp>0. We may choose mj=jj in evaluating @), Aarars 


ji —je=RE>psmj—me=ji— ME. - 


Hence j¢< me and, since also jp> mr, we must have jg—m,. ‘Thus only one term occurs in 


the summation (2), viz. 
Oi. J) =PUis Ji Its Je)» Beas (6) 


‘Then, for ji<jr, (3), (6) give OC, j) =[2ie+1)/2i+1))POe, St: Ji Ji)- 
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The discrepancy. Greuling’s equations (1), (3) contain a summation over m which 
corresponds to our summation over » or mg. He does not explain in detail his derivation of. 
his equations (11), but the general implication of his statement in the second paragraph of 
p. 575 of his paper seems to be wrong. For j; >j¢ Greuling’s summation over m should be 
equivalent to Marshak’s result. Marshak’s summation reduces to one term when mj=jj). 
but the result would be the same whatever value of m; were chosen and however many were 
the terms in the summation. Thus there is no need to multiply half-lives by 2(ji—Jp) +1 
if Greuling’s equations (1), (3) are correctly applied. 

The number of terms in the summation may be easily calculated. By equation (5),. 
for given mi, ji, jf, # may take the values mi—jr, mi—jet+1, ...-, mi+je, and, in order 
that |u| <k, we have by (4) — |ji—je| <mi—set+A<|si—Je|, where A is any integer 
satisfying 0< A< 2)¢. 

Thus there are 2|jj—je|-+1 or |ji—je| —mitJe+1 terms according respectively as 

|ji--je| —mitje20 or <0. Thus, if ji<je, there are 2(je—Ji) +1 terms, while if jj >7¢ 
there are 2(j;—j¢) +1 or jj mj+1 terms according as 2j¢ 2 or <jjt+mj. Thus, if jj >J¢ and 
m;=ji, there is just one term. 


Department of Natural Philosophy, C. STRACHAN.. 
University of Aberdeen. 
ist September 1948. 
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The Attenuation of Ultra-High Frequency Electromagnetic 
Radiation by Rocks 


We were interested by Mr. Cooper’s paper (1948) on the attenuation of very short: 
radio waves in rocks, and would like to draw attention to some very similar work which we 
carried out at the National Physical Laboratory in 1942, the results of which were published 
later (McPetrie and Saxton 1945). 

We made measurements on wavelengths of 2 and 5 metres of the field strength produced 
in two specially excavated tunnels by a horizontal aerial fed from a small power transmitter 
which was erected over each tunnel in turn. The two tunnels were made as nearly identical 
in structure as possible, except that one was at a depth of 1-15 metres and the other at 
1-85 metres. It was hoped that any disturbing “‘ tunnel effect” could be eliminated by 
making comparative measurements in the two tunnels when subsidiary tests made above 
ground showed that the same field existed at the surface of the earth immediately above 
the tunnels. This seems to us a more reliable procedure than that followed by Mr. Cooper, 
in which an attempt is made to evaluate this tunnel effect theoretically. Any such attempt 
-must necessarily be uncertain. In Mr. Cooper’s experiments, fortunately, the transmission 
path through the earth is sufficiently long for the absorption in the earth to be very much 
greater than such losses as may be expected, on any reasonable hypothesis, to arise from the 
effects of the tunnel. 

Our two tunnels were close together on the same site, which was of sandy loam overlying 
gravel, and the observations thus gave directly the attenuation in a path of 0-7 metres 
through this medium. We also made measurements of the. reflection coefficient of the 
earth’s surface on the same site, and these, together with the absorption measurements, 
enabled us to determine both the dielectric constant and the conductivity of the ground, 
assuming uniform properties with depth for the latter. Since this assumption is not 
strictly true, the values so obtained must be regarded only as reasonable mean values. 
Our results indicated that the attenuation of both 2-metre and 5-metre waves was about 


12 db/m., which is not very different from the value of 10-5 db/m. for Pennant sandstone 
quoted by Mr. Cooper. 
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The actual figure for the attenuation obtained, as the author states, depends largely 
upon the amount of water present. We should like to point out, however, that Mr. Cooper’s 
statement that at a frequency of 200 Mc/s. the dielectric losses in water are still negligible 
in comparison with the conductivity loss is not strictly accurate. In fact, at this frequency 
for the usual amount of dissolved salts occurring in fresh water, which leads to an ionic 


conductivity of about 10% E.s.u. (~107!8 E.M.u.), the dielectric and conductivity losses 
are just about equal. 


Signals Research and Development Establishment, J. S. McPETRIE. 
Somerford, Christchurch, Hants., J. A. SAXTON. 
and 


Radio Division, National Physical Laboratory, 
Teddington, Middlesex. 
21st July, 1948. 
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I was very grateful to Drs. McPetrie and Saxton, with whose work I am, of course, 
familiar, for their comments on my paper. I must apologize for not having referred to 
them in it, but I should like to point out that my work was concerned with the electrical 
properties of large thicknesses of rock im sitw, and not with the properties of the surface 
soil-layer. JI had hoped that the consolidated Pennant Sandstone which I used would 
have held less water and therefore show a smaller attenuation than the sandy loam mentioned 
in their work. I agree that their procedure was more reliable, but the experimental diffi- 
culties with which I had to contend prevented any more satisfactory arrangements from 
being made. In fact I find some consolation in the agreement obtained with their results. 

As a matter of interest I also observed a rainfall effect similar to that mentioned in their 
paper during my work ; on the last day of the experiments in Pontypridd tunnel rain fell, 
but it was not until the following day that any diminution was observed in the echo-strength 
inside the tunnel. This effect was not investigated sufficiently well for inclusion in the 
paper, however. 

Finally, I am indebted to Drs. McPetrie and Saxton for pointing out the error mentioned 
in the last paragraph of their letter. 


Department of Geodesy and Geophysics, R. I. B. Cooper. 
Downing Place, Cambridge. 
3rd August 1948. 


The Performance of a New Radio-Frequency Ion Source 


The observations of Hertzberg (1927) on the low intensity of the molecular spectrum 
of an electrodeless ring discharge in hydrogen suggested that this method of excitation 
should prove very effective in achieving a high proton percentage in an ion source. ‘This 
mode of excitation has been used in a recently described ion source (Bayley and Ward 1948). 
The characteristics of such an ion source are of great interest to users of accelerated ion 
beams. The radio-frequency ion source which has been developed in this laboratory 
is reliable and gives high ion currents, a high atomic content and ions of uniform energy. 

A high concentration of atomic hydrogen is maintained in a pyrex discharge tube by 
reducing the exposed metal surfaces to a minimum and producing a high degree of ionization 
by a coil round the tube excited by a 20.Mc/s. 150 watt oscillator. ‘The discharge is self- 
starting provided the probe voltage is off. A tungsten wire anode is isolated from the main 
discharge by a constriction. ‘This reduces the flow into the discharge space of molecular 
hydrogen formed by recombination on the anode. The duralumin probe electrode enclosed 
in a glass cylinder projects into the discharge space as it has been found (Hertzberg 1927) 
that the concentration of atomic gas and hence the proton percentage diminishes towards 
the walls. Because of the high concentration of positive ions the exposed probe surface 
can be reduced to small dimensions and recombination reduced to a minimum. 
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The form of the electric field resulting from the presence of the glass cylinder around 
the probe focuses the ions into the canal. This mode of focusing has been demonstrated 
before (Thonemann 1946). 

Figure 1 shows the mechanical details of the source used for making measurements. 
Approximately 50 watts of high-frequency energy are dissipated in the tube. A small 
air-blower is used to keep the wax joint cool. 
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‘Figure 1. Assembly of the radio-frequency Figure 2. Comparison of the theoretical contour for 
ion source. a homogeneous beam with the D,* ion contour. 
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Table 1 shows the mass spectrum of a deuterium ion beam with about 7% hydrogen 


present. ‘The various masses and break-up products are easily identifiable. The width 
of the D,* peak is approximately $ ma. 


Table 1. The mass spectrum of a deuterium ion beam 


Magnet : Percentage 
current (ma.) Species abundance 

29 D,* from dissociation of D3+ 0-9 

32 H,*; D,* from dissociation of D,* 6:0 

54 Die and skies 79-0 

60 D,* from dissociation of D,* 1-0 

62 EDiwandi Ess 0:9 

80 D,* and H,D* 8-8 

90 (80D) 1-4 

99 IDs 2-1 


The variation of proton and deuteron percentage with pressure in the source is shown 
in tables 2 and 3. The discharge extinguishes at the low pressure limit. Hydrogen was 
supplied through a heated palladium tube. The differences in atomic percentage between 


deuterium and hydrogen may be due to the difference between the radio-frequency powers 
absorbed in deuterium and in hydrogen. 
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When a new source is first operated the discharge may change suddenly from its normal 
brilliant red to a dull white and the proton current may fall by a factor of 10. This is 
probably due to the evolution of hydrocarbons and other gases from electrodes under 
positive ion bombardment in the accelerator. After the system has been outgassed the 
proton current remains stable and shows no deterioration after months of intermittent use. 

The homogeneity of the deuterons after a 10 kv. acceleration was measured on a mass 
spectrometer. A Faraday screen was introduced to eliminate any possible modulation 
of the ion energy by the electrostatic r.F. field. The discharge had then to be initiated 
by bringing a Tesla coil close to the upper section of the source. Figure 2 shows the 
contour for the D,+ ion. The theoretical contour from a homogeneous source, taking 
into account the aberrations in the spectrometer, is given for comparison. It is concluded 
that 90 % of the ions have an energy spread of less than +20 electron volts. The magnitude 


Pressure (microns) 14 16 18 20 23; 2S 26 28 
H,t 84 85 88 89 90 90 92 91 
H,t 14 11 10 8 7 8 5 6 
H;* 2 4 2 8 3 2 3 3 
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of the tail on the low energy side of the experimental curve increases with pressure in the 
accelerator and spectrometer and is attributed to scattering and energy loss by collision 
in the high vacuum side of the apparatus. 

In adapting this ion source for use on our high voltage accelerators we have found 
it more convenient to connect the gas leak and Pirani gauge through the base of the ion 
source. ‘To increase the focusing efficiency the outside diameter of the probe has been 
reduced from 1/4 in. to 3/16in. The same high atomic percentage has been measured 
with this arrangement. With 3 kv. extraction volts and 10 ma. probe current no difficulty 
has been found in focusing a }ma. beam of ions in a 6 mm. diameter spot after passing 
through 3 metres of the accelerating tube. 

P. C. TTHONEMANN. 
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The Absorption Coefficient of the 8 Rays from 7 
in Aluminium 


During the course of some measurements using radio sulphur 32S it became necessary 
to make an estimate of the correction due to self absorption of the Brays. ‘The end point 
of the B-ray spectrum is known (Libby and Lee) but as this lies outside the range of F eather S 
empirical formula it is hoped that the following experimental determination of the absorption 
coefficient in aluminium may be of use to other workers using counters with cylindrical 
See caoements were made with a demountable counter ; it differs only in details 
from that used by Collie and Morgan in 1939. The active barium sulphate was distributed 
as evenly as possible as a film whose mean weight was less than 1/10 mg/cm?* on the inside 
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of a movable glass cylinder. ‘The counter was filled with a carefully controlled air alcohol 
mixture using a bellows manometer of the type described by East and Kuhn (1946). A 
long series of measurements showed that variations in the observed activity of a single 
specimen due to refilling the counter and replacing the specimen did not amount to more 
than 3% (including the statistical error). The aluminium foil used as an absorber was 
wrapped round the perforated cathode of the counter and unwrapped layer by layer. On 
examining the pulses with a cathode-ray tube it was observed that the aluminium foil 
caused the production of a small proportion of double counts. ‘These were insufficient 
to interfere with the stability of the counter at the highest rate of counting used. 

The absorption curve was calculated on the assumption that the number of double 
counts was a constant proportion of the true activity. The correctness of this assumption _ 
was verified experimentally by varying the resolving time of the counting circuit and com- 
paring the results with those obtained with a counter free from doubles. ‘The absorption 
was found to be exponential with an absorption coefficient w=780+25 cm™. The 
greatest uncertainty of the determination is due to the difficulty of measuring the exact 
weight of foil through which the 8 particles have passed. Using this figure for p, the 
activity of a bare source will be halved by 2:4 mg/cm? of aluminium. This result may be 
compared with the observation of Hendricks and his collaborators (1943) that the halving 
mass in mica is 3 mg/cm?. The small difference between these results is not more than 
would be expected when the very different geometrical arrangements are taken into account. 

We have to thank Professor Peters for putting the radio sulphur at our disposal, the 
Board of Education for a maintenance grant and Lord Cherwell for extending to us the 
facilities of his laboratory. 


The Clarendon Laboratory, Oxford. C. H. Commit: 
7th September 1948. M. A. GRACE. 
(First communicated 3rd August 1946). 
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The Effect of Room Temperature Radiation on the Infra-Red 
Response of Lead Telluride Photoconductors 


One of the principal objects in studying the photoconductivity of semiconductors is 
to search for substances which show photosensitivity to radiation as far into the infra-red 
as possible. The practical importance of these substances is in their ready application 
to infra-red detection, and in particular to detection for the purposes of infra-red spectro- 
scopy. It has been generally recognized that as the long wavelength limit of photoconduc- 
tivity is extended, the temperature at which a detector may profitably be used must be 
decreased in order to cut down thermal excitation in the material. However, if the thermal 
radiation from the surroundings of a photoconductor contains appreciable energy at wave- 
lengths to which the detector is photosensitive, there will be an effective background of 
illumination which will increase the conductivity even when no source of illumination is 
deliberately applied. This will be true at whatever temperature the photoconductor is 
maintained, unless it is surrounded by an enclosure at the same, or lower temperature 

Lead telluride, in the form of thin films, has been found to show photocondusuyny 
at least as far into the infra-red as 5-5 », and it has been shown that at — 190° c. its perform- 
ance as an infra-red detector surpasses that of a good Hilger-Schwartz thermocouple 
for wavelengths shorter than 4:74 (Simpson, Sutherland and Blackwell 1948). It has 
recently been discovered that the conductivities of all photoconductive lead telluride films 
are much less when they are cooled to —190° c. and surrounded by an enclosure at that 
temperature than when they are cooled to —190° c. but surrounded by an enclosure at 
‘room temperature. ‘That this effect is due to the thermal radiation which in the latter 
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‘case is emitted by the enclosure at 15° c. is demonstrated in the figure where the logarithm 
of the conductivity of a film is plotted against the reciprocal of the temperature. Curve A 
‘shows the conductivity that is observed when the film is allowed to warm up from —190° c. 
to room temperature, the photoconductor being exposed all the time to an enclosure at 
15°c. It will be observed that below —170°c. (1/T>9-7x 10~*), the conductivity 
decreases with increasing temperature; this in itself is anomalous for a semiconductor. 
The anomaly is removed in curve B where the conductivity of the same film is shown 
under conditions in which the thermal radiation was excluded by surrounding the film 
iby an enclosure at —190°c. The conductivity at —190°c. was more than six times 
greater when the film was exposed to the room temperature radiation; this effect decreased 
rapidly as the temperature was increased, and was negligible at temperatures above — 140° c. 
The effect of the thermal radiation on the photoconductive response may best be 
illustrated by considering the following equation which is found experimentally : 
(o;/oq)?=1+-81, 
where oq is the conductivity at —190° c. with the thermal radiation excluded, and 9; the 
‘conductivity under an intensity of illumination J watts/em?. The constant b depends 
upon the wavelength of the incident light 
and was 2:2 10° watt-! cm? at 2:2» for 10“ 
the film shown in the figure. Using this 
equation, the voltage change across a film 
caused by the absorption of radiation while 
the film is passing a fixed current may be 
readily calculated. It is found that for an 
exciting current of 1a. the voltage output 
is 9-2 v/uw/cm?, A=2-2 », when the film is 10°F eek 
entirely shielded from thermal radiation. 
The output for the same exciting current 
when the thermal radiation is not excluded o B 
is only 0-053 v/uw/cm?. 4 p ; r — 

In order to assess the effect of the ther- Reciprocal Temperature 
mal radiation on the ultimate limit of sensi- 
tivity of lead telluride further data are required concerning the noise characteristics of the 
films. Clark Jones (1947) has shown that fluctuations in the radiation field set the ultimate 
limit for any detector operating under prescribed conditions. It has been shown here that 
by excluding the background of thermal radiation the signal output may be increased; 
the signal-to-noise ratio will be further increased by this means if the noise output of the 
film is predominantly due to statistical fluctuations in the background radiation. ‘There 
is at present insufficient evidence for it to be certain that this is in fact the case. 

The spectral response curve for lead telluride (Simpson, Sutherland and Blackwell 
1948) has so far only been obtained with films exposed to room temperature radiation ; 
there is a possibility that under these conditions the photosensitivity at long wavelengths 
is at least partly saturated. If this were so the photoconductive response would be increased 
even more at long wavelengths than at the wavelength 2:2 » which has been used as standard 
in this discussion. It seems possible that the changes in the spectral response of PbS, 
PbSe and PbTe obtained by Moss (1948) by operating at very low temperatures are not 
unconnected with the above results. ; ; . 

I am indebted to Dr. G. B. B. M. Sutherland and Dr. E. F. Daly for helpful discussions. 

O. Srvpson. 
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Cavendish Laboratory, Cambridge. 
12th August 1948. 
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REVIEWS OF BOOKS 


Radar Beacons, edited by ARTHUR Roperts. Pp. xx+489 with 246 figures. 
(New York : McGraw-Hill Book Co. Inc., 1947.) $6. 


After world war II the Radiation Laboratory of the Massachusetts Institute of Tech- 
nology undertook the heavy task of preparing for publication in book form the war-time 
developments of Radar and Radar Techniques. Several volumes of this series have already 
appeared and are proving of real value to the specialist and student alike. ‘The present 
volume is Volume III of the series and deals with the important subject of Radar Beacons. 

A radar beacon is a small receiver transmitter apparatus which is used at the target to be 
located. It is actuated (or triggered, to use the correct technical term) by the pulsed signals. 
sent out by the radar set, and in turn the beacon sends out pulsed signals which may be 
received at the radar set. The time delays in the beacon system can usually be made very 
small and constant so that range finding by timing the interval between the radar transmitted. 
pulse and the pulse echo received back from the beacon is possible as for “‘ primary radar ”’. 
This latter depends on the target itself reflecting back some of the incident energy. ‘The 
beacon system (called by Watson-Watt “‘ a secondary radar system ”’) has, however, several 
advantages when circumstances allowits use. Thus it allows target location in circumstances. 
where ground clutter would prevent the operation of a primary radar system, greater range 
and, if the beacon transmissions are coded, positive identification of the target echoes. For 
these, and for other reasons, the beacon system was much used during the last war for general 
navigation of aircraft and surface vessels, for approach systems at airfields, for paratroop 
dropping, for bombing aids and for other purposes, and it is now finding application in 
civil aviation and in maritime service. ‘The present book attempts to give a survey of these 
developments and, it can be said at once, that the 24 authors who have collaborated in this 
work have succeeded in providing a very valuable treatise which should be of special interest 
to all those who want to obtain a thorough understanding of the present state of the beacon 
art. 

The book is divided into four parts : Basic Considerations, Beacon Design, Interrogator 
and System Design, and Beacons in the Field. ‘The first part is especially valuable and gives 
a thorough treatment of the principles of the subject. After a brief introductory section 
dealing with the types and functions of beacons it gives a detailed discussion of the require- 
ments of systems using beacons, considering in turn range requirements, propagation and 
coverage, choice of frequency, use of coding and communications, and questions of traffic 
and capacity. All this is explained in a very lucid manner and should be easily assimilated 
by the would-be beacon designer. 

Part If deals with beacon design and includes sections on the design of all the component 
parts of the beacon. There are chapters on the design of the R.F. components, beacon 
receivers, interrogation coding, response coding, beacons as pulse communication systems, 
the design of modulators and transmitters and the associated problems of power supplies and 
performance testing. This is all very well explained, particularly so far as the components 
of microwave systems are concerned, but, in some cases, the treatment of the component 
parts of the longer wave systems is dealt with far too briefly. 

Part III gives details of the Interrogator and system design and includes chapters on 
radar design for beacon operation, the design of interrogator-responsors and_ the 
design of typical beacon systems. This latter chapter is well worth careful study as an 
outline is given of the application of the principles discussed in Part I to meet requirements 
of different types of beacons and will therefore be of particular value to the designer. 

Part IV discusses beacons in the field and includes chapters on installation, operation 
and maintenance. Questions of this sort tend to be forgotten in many books and it is 
pleasing to a reviewer to find these subjects so well presented. 

The book can be recommended with every confidence. DENIS TAYLOR. 
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